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This  thesis  examined  the  feasibility  and  a  method  of 
extending  continuous  domain  Quantitative  Feedback  Theory  (QFT) 
flight  control  system  designs  to  the  discrete  domain.  Analog 
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Abstract 

This  inTveatigation  examines  the  feasibility  and  a  method 
of  extending  continuous  domain  Quantitative  Feedback  Theory 
(QFT)  flight  control  system  designs  to  the  discrete  domain. 

The  results  of  two  previous  OFT  analog  design  efforts  are 
modified  for  a  digital  implementation.  The  first  design 
effort  is  for  the  KC-135  transport  aircraft.  Robust  analog 
fixed  compensators  are  designed  for  three  different  flight 
conditions.  The  second  design  effort  is  for  the  AFTl/F-16 
aircraft.  In  this  design,  parameter  variation  is  due  to  both 
varying  flight  conditions  and  dif^ferent  aircraft  configura- 
tions  caused  by  failed  surfaces.  _ 

This  study  is  accomplished  in  two  phases.  In  the  first 
phase  only  the  effects  of  sampling  are  included.  The  compen¬ 
sators  are  discretized  using  two  transformation  techniques 
after  they  are  modified  to  counteract  the  effects  of  the 
sampler  and  zero  order  hold  (ZOH).  Simulation  results  show 
that  the  Tustin  transformation  using  a  sampling  rate  of  40  Hz 
provides  good  results,  as  compared  to  the  analog  design,  for 
the  KC-135  aircraft.  For  the  AFTI/F-16,  acceptable  results 
are  obtained  using  the  Tustin  transformation  and  a  sampling 
rate  of  60  HZ,  but  only  for  two  out  of  four  flight  conditions. 

In  the  second  phase  both  sampling  and  finite  wordlength 


effects  are  considered.  A  computer  program  is  developed  to 
simulate  the  effects  of  finite  wordlength.  It  is  suffi¬ 
ciently  general  to  allow  the  specification  of  rounding  or 
truncation,  the  controller  wordlength  and  the  number  of 
desired  binary  digits  to  the  right  of  the  radix  point. 
Simulations  are  conducted  using  various  wordlengths  and 
numbers  of  digits  to  the  right  of  the  radix  point  as  well  as 
appropriately  quantized  controller  coefficients. 

The  conclusion  drawn  from  this  study  is  that  the  Tustin 
transformation  provides  good  discrete  system  performance  as 
compared  to  the  continuous  OFT  design,  if  the  compensator  real 
zeros  and  poles  are  located  at  S  -T/2.  To  ensure  the  same 
system  performance  guarantees  as  provided  by  QFT  in  the 
continuous  domain,  it  is  recommended  that  the  OFT  design  be 
performed  directly  in  the  discrete  domain,  as  opposed  to 
extending  analog  QFT  designs. 


STUDY  OF  THE  EFFECTS  OF  DISCRETIZING 
QUANTITATIVE  FEEDBACK  THEORY 
ANALOG  CONTROL  SYSTEM  DESIGNS 

I .  Introduction 

Quantitative  Feedback  Theory  developed  by  Dr.  Isaac 
Horowitz  has  been  used  successfully  to  design  "robust*  control 
systems  (1;  2;  3).  A  robust  control  system  is  generally 
considered  to  be  one  which  has  reduced  sensitivity  to  plant 
parameter  variations  and  disturbances,  i.e.  less  than  in  an 
open  loop  system.  These  parameter  variations  result  from 
uncertain  plant  parameters  or  uncertain  external  distur¬ 
bances.  Plant  parameter  variations  occur  when  an  aircraft  is 
operated  at  a  flight  condition  other  than  the  flight  condition 
used  to  generate  the  linear  plant  model  incorporated  in  the 
basic  control  system  design.  Loss  or  failure  of  a  flight 
control  surface  can  also  cause  large  plant  parameter  varia¬ 
tions  . 

There  are  many  advantages  to  implementing  control 
compensators  using  a  digital  controller.  Digital  control¬ 
lers  have  improved  reliability  and  maintainability  over 
analog  compensators.  A  great  degree  of  flexibility  can 
also  be  realized  with  a  digital  controller  implementation. 

A  completely  new  control  law  can  be  implemented  simply  by 
making  a  software  change  instead  of  a  hardware  change. 

Several  control  laws  can  be  stored  simultaneously.  Depending 


on  mission  requirements,  a  particular  control  law  can  be 
"switched  on"  by  use  of  a  software  algorithm. 

Previous  applications  of  Quantitative  Feedback  Theory  (QFT) 
have  been  limited  to  the  design  of  control  systems  in  the 
continuous  domain.  In  order  to  realize  the  inherent  advant¬ 
ages  of  a  digital  controller,  such  as  those  cited  above, 
either  QFT  or  QFT  designs  must  be  extended  to  the  discrete 
domain.  A  digital  implementaion  of  a  QFT  analog  design  was 
done  for  the  lateral  mode  of  the  YF-16  CCV  aircraft  (4:174), 
but  finite  wordlength  effects  were  not  considered. 

In  this  thesis  investigation,  the  results  of  two 
previous  QFT  design  efforts  (1;  2)  are  extended  to  the 
discrete  domain.  An  analysis  of  these  discrete  designs  are 
made  to  determine  if  they  meet  the  performance  objectives 
of  he  original  designs. 

This  chapter  presents  the  pertinent  background  informa¬ 
tion,  statement  and  scope  of  the  problem,  assumptions,  the 
approach  and  the  sequence  of  thesis  presentation. 

I-l  Background 

Several  difficult  control  problems  involving  highly 
uncertain  parameters  have  been  solved  by  Dr.  Isaac  Horowitz 
using  the  QFT  technique  (5:287-309;  6:81-106).  In  addition  to 
Dr.  Horowitz's  work,  several  AFIT  Master's  theses  have  applied 
the  QFT  technique  to  flight  control  problems  involving 
uncertain  parameters,  two  of  which  are  used  to  provide  the 
starting  point  for  this  thesis  effort  (1;  2).  These  two 


theses  are  discussed  next: 


The  objective  of  the  first  thesis  effort,  by  Captain  j 

Russell,  was  to  design  two  robust  controllers  for  the  KC-135 
transport  aircraft  (1:2).  The  first  design  involved  two 
inputs  and  two  outputs,  and  controlled  the  lateral  mode  of  the 
aircraft  during  two  separate  maneuvers.  The  second  design 
involved  three  inputs  and  three  outputs,  and  controlled  the 
longitudinal  mode  during  one  maneuver.  In  Russell's  thesis, 
the  linearized  lateral  and  longitudinal  equations  of  motion 
are  (assumed)  decoupled  and  treated  separately  through¬ 
out  (see  Chapter  II).  Both  controllers  were  designed  to 
operate  over  the  entire  flight  envelope  of  the  aircraft.  All 
performance  specifications  were  met  for  all  three  flight 
conditions  in  both  the  lateral  and  longitudinal  designs  using 
only  fixed  compensation  (1:99). 

The  objective  of  the  second  thesis,  by  Major  Arnold,  was 
to  investigate  QFT  "inherent  reconfiguration"  design  for  the 
AFTI/F-16  aircraft  (2:2).  In  this  context,  inherent  reconfig¬ 
uration  means  the  ability  of  the  control  system  to  maintain 
aircraft  control  even  under  one  or  more  simultaneous  losses  or 
failures  of  control  surfaces,  inherent  reconfiguration  is 


accomplished  by  a  fixed  controller  using  individual  control  of 
all  (including  redundant)  flight  control  surfaces  of  the 
AFTI/F-16  aircraft.  Thus  in  the  above  work,  the  two  flaperons 
and  the  two  horizontal  tails  were  individually  controlled. 

The  plant  model  explicity  displayed  the  effect  of  each  control 
surface  on  the  two  chosen  output  variables,  pitch  rate  (q) 
and  roll  rate  (p)  (see  Chapter  II).  Two  command  input  sources 


and  were  assumed,  and  quantitative  basically  noninter¬ 
acting  performance  was  specified.  Not  only  was  the  control¬ 
ler  designed  to  operate  over  the  entire  flight  envelope,  as  in 
Russell's  thesis,  but  it  was  also  designed  to  operate  over  6 
different  failure  combinations.  Not  all  performance  specifi¬ 
cations  were  met  in  Arnold's  design  as  the  settling  time 
tolerance  was  exceeded  for  several  aircraft  configuration/ 
flight  condition  combinations  (2:58).  QFT  guarantees  that  the 
original  performance  tolerances  are  satisfied  despite  uncer¬ 
tainty  (5).  No  exact  reason  is  given  for  this  inconsistency, 
but  Arnold  does  speculate  that  one  possible  reason  is  that 
approximations  were  made  for  the  compensators  g^  and  g2  so 
their  order  could  be  reduced  (2:60).  The  order  reduction  was 
required  due  to  limitations  of  the  simulation  program 
utilized. 

The  number  of  compensators  for  each  of  the  designs 
discussed  above  is  determined  by  the  dimensionality  of  the 
system  and  the  number  of  "non-zero"  inputs.  For  instance, 
for  a  three  input-three  output  system  (3  loops  or  chan¬ 
nels),  as  in  Russell's  longitudinal  design,  three  feed¬ 
forward  path  (inner  loop)  compensators  and  one  prefilter 
are  required  to  complete  the  control  design.  In  Russell's 
longitudinal  design  the  only  non-zero  input  was  pitch  angle. 
Note  that  in  Russell's  design,  the  term  "robust"  refers  to  the 
plant  parameter  variation  due  only  to  different  flight 
conditions.  In  Arnold's  design,  "robust"  means  parameter 
variation  due  to  both  varying  flight  conditions  and  the 


different  configurations  of  the  aircraft  caused  by  failed 
surfaces . 

1-2  Problem 

The  main  purpose  of  this  effort  is  to  investigate 
the  suitability  of  continuous  domain  QFT  control  system 
designs  for  digital  implementation.  This  general  problem  is 
separated  into  several  smaller  but  related  problems.  The 
results  of  Russell's  and  Arnold's  designs  are  extended  to  the 
discrete  domain.  This  is  accomplished  by  discretizing  the 
compensators  usina  various  transformation  techniques.  The 
performance  of  the  digital  controllers  are  examined  by  use  of 
a  very  large  wordlength  digital  computer  simulation  program 
(7).  Performance  is  judged  via  comparisons  to  the  original 
time  domain  specifications  and  performance  achieved  in 
references  1  and  2.  Implementation  issues,  such  as  the 
effects  of  finite  word  length  and  sample  time,  are  discussed. 

Recommendations  are  made,  based  upon  the  above  analysis,  to 
aid  the  design  of  s-domain  compensators  using  QFT  so  that 
performance  objectives  are  met  when  the  designs  are  implement¬ 
ed  digitally. 

1-3  Scope 

This  study  extends  continuous  domain  QFT  flight  control 
designs  (1;  2)  to  the  digital  domain.  The  digital  design 
starting  points  are  the  continuous  designs  of  Captain  Russell 
(KC-135)  and  Major  Arnold  (AFTI/F-16).  Only  the  rigid  body 
aircraft  is  considered  for  the  AFTl/F-16.  For  the  KC-135,  the 


first  two  body  bending  nodes  are  also  included.  This  study 
involves  only  constant  rate  sampling.  The  following  informa¬ 
tion  is  provided  concerning  the  above  two  design  efforts. 

For  the  KC-135,  the  digital  designs  are  done  for 
the  decoupled  longitudinal  and  lateral  modes.  The  performance 
of  the  lateral  digital  control  design  is  examined  using  two 
maneuvers:  a  coordinated  turn  and  a  sideslip.  The  longitu¬ 
dinal  digital  control  design  involves  a  pitch  pointing 
maneuver.  The  KC-135  designs  are  for  the  following  flight 
conditions  (F.C.)  (1:3); 

F.C.  #1:  High  altitude  cruise 

(45,000  ft  at  mach  0.77). 

F.C.  #2:  Medium  altitude  cruise 

(28,500  ft  at  mach  0.77). 

F.C.  #3;  Approach 

(sea  level  at  mach  0.21). 

The  AFTI/F-16  design  is  for  the  following  four  F.C.s  (2:97): 


F.C. 

#1: 

Approach 

(30  ft  at  mach  0.2). 

F.C. 

#2: 

High  altitude/medium  speed  cruise 
(30,000  ft  at  mach  0.6). 

F.C. 

#3; 

Medium  altitude/medium  speed  cruise 
(20,000  ft  at  mach  0.9). 

F.C. 

#4: 

High  altitude/high  speed  cruise 
(30,000  ft  at  mach  1.6). 

To  incorporate  the  reconfiguration  aspect  into  the  design, 
the  following  control  surface  configuration  (CSC)  modes  are 
utilized  (2:4): 

Mode  1  -  All  four  surfaces  are  operating  normally. 
Mode  2  -  One  horizontal  tail  is  failed 
Mode  3  -  One  flaperon  is  failed. 


Mode  4  -  One  horizontal  tail  and  one  flaperon  are 
failed  {same  side). 

Mode  5  -  One  horizontal  tail  and  one  flaperon  are 
failed  (opposite  sides). 

Mode  6  -  Both  flaperons  are  failed. 

A  failed  surface  is  defined  to  be  in  the  neutral  reference 
position . 

Although  a  detailed  discussion  of  the  OFT  design 
technique  is  not  a  purpose  of  this  investigation,  two  appendi 
ces  are  included  which  do  present  this  material.  They  are 
included  as  a  convenience  to  the  reader  and  because  of  the 
relatively  small  amount  of  literature  available  on 
Quantitative  Feedback  Theory.  The  single  input-single  output 
(SISO)  system  design  philosophy  is  presented  in  detail  in 
Appendix  A  and  the  multiple  input-multiple  output  (MIHO) 
design  technique  is  presented  in  Appendix  B.  The  emphasis  in 
Appendix  B  is  on  the  (2x2)  and  (3x3)  MIMO  system  design. 

1-4  Assumptions 

Several  assumptions  are  made  in  the  course  of  this 
investigation.  Seven  of  these  assumptions  are  the  same  as 
those  made  by  Russell  and  Arnold  (1:3;  2:4).  They  gener¬ 
ally  cover  various  aspects  of  generating  the  linear,  time- 
invariant,  decoupled  aircraft  equations  of  motion.  Russell 
and  Arnold  also  assume  that  ’continuous  time’  performance 
simulations  executed  on  a  large  wordlength  digital  computer 
provide  realistic  responses  of  the  aircraft  motion.  It  is 
assumed  that  the  design  work  done  by  Russell  and  Arnold  in 
their  respective  QFT  design  efforts  is  correct  as  presented. 


although  some  cursory  checks  of  their  work  are  made. 

This  study  involves  various  computer-aided-design 
(CAD)  and  computer-aided  simulation  tools  to  examine  the 
performance  of  the  digital  control  system  designs.  These 
simulations  are  carried  out  on  large  wordlength  digital 
computers.  It  is  assumed  that  these  simulations  provide 
realistic  responses  when  an  appropriate  time  relationship 
has  been  established  between  the  digital  controller  and  the 
"analog  plant." 

It  is  assumed  that  the  reader  of  this  thesis  is  at 
least  familiar  with  the  work  done  by  Russell  and  Arnold  in 
references  1  and  2,  respectively.  This  is  important,  since 
the  digital  controller  performance  is  compared  directly  to 
their  specifications  and  results,  in  addition,  a  general 
understanding  of  the  QFT  design  technique  is  assumed  on  the 
part  of  the  reader. 

1-5  Approach 

This  digital  design  study  involves  two  general  phases. 
Phase  1  discusses  the  theoretical  aspects  of  the  digital 
design  and  examines  the  performance  of  the  theoretical 
design.  In  this  context,  "theoretical"  means  that  finite 
wordlength  effects  are  not  considered.  Phase  2  discusses 
various  implementation  issues  associated  with  a  finite  word 
length  digital  controller.  A  performance  simulation  is  done 
which  takes  into  account  the  effects  of  finite  word  length. 
With  the  above  design  phases  in  mind,  the  digital  design 
approach  incorporates  the  following  steps; 


The  first  step  of  the  theoretical  phase  developes 
a  suitable  system  model  in  the  form  of  a  system  block  diaqram 
including  samplers  and  zero  order  holds  (ZOH),  The  second 
step  developes  a  means  co  counteract  the  lag  effects  of  the 
ZOH.  Next,  the  sample  rate  is  determined  from  a  theoretical 
standpoint  and  this  sample  rate  is  then  compared  to  sample 
rates  that  can  be  practically  realized.  The  fourth  step  is  to 
transform  the  modified  compensators  to  the  discrete  domain 
using  the  Tustin  and  standard  Z-transf orms .  The  effects  of 
■folding"  and  "warping"  of  poles  and  zeros  are  then  discus¬ 
sed.  The  last  step  of  the  theoretical  phase  involves  a 
simulation  of  the  theoretical  digital  design,  using  various 
sample  rates,  as  necessary,  to  determine  the  theoretical 
performance.  The  results  are  then  compared  to  the  specifi¬ 
cations  and  continuous  domain  performance. 

The  first  step  of  the  implementation  phase  developes  a 
"hybrid"  simulation  program  using  existing  simulation  tools. 
The  second  step  performs  a  hybrid  simulation  using  various 
wordlengths.  The  results  of  this  simulation  are  then  compared 
to  the  original  specifications  and  continuous  domain  perfor¬ 
mance  results. 

Throughout  both  the  theoretical  and  implementation 
phases,  discussion  is  included  with  suggestions  for  additions 
or  changes  to  the  QFT  design  technique,  to  improve  the  results 
obtained  in  an  eventual  digital  implementation.  These 
guidelines  are  summarized  at  the  end  of  this  thesis. 


1-6  Presentation 

This  thesis  consists  of  seven  chapters.  Chapter 
II  contains  a  detailed  discussion  of  the  system  models  to  be 
used.  Chapter  III  presents  the  detailed  theoretical  digital 
design.  The  theoretical  simulation  ("infinite  wordlength")  of 
the  digital  controllers  is  contained  in  Chapter  IV.  Chapter  V 
discusses  the  hybrid  simulation  program  and  digital  imple¬ 
mentation  issues.  The  results  of  the  hybrid  simulation  are 
presented  in  Chapter  VI,  Chapter  VII  contains  the  conclusions 
and  recommendations  for  both  the  digital  implementation  of  the 
QFT  design  technique  and  future  study  in  this  area. 


I I .  System  Models 


k 


I 


I I-l  Introduction 

In  this  chapter  the  various  models  used  in  this  thesis 
are  presented.  First  a  generic  pseudo-continuous-time  (PCT) 
single  input-single  output  (SISO)  system  model  is  developed 
from  a  continuous  time  model.  The  development  of  the  2x2 
(lateral)  KC-135  PCT  model  is  shown  next.  The  form  of  the 
aircraft  equations  given  in  reference  1  are  incompatible  with 
the  PCT  model  configuration  developed,  therefore  the  form  of 
the  KC-135  equations  are  changed  accordingly.  The  rigid  and 
elastic  body  3x3  (longitudinal)  KC-135  ..  JT  model  is  then 
developed.  As  with  the  lateral  aircraft,  the  longitudinal 
aircraft  equations  must  be  massaged  into  the  proper  form. 
Finally  the  F-16  continuous  time  model  is  discussed  and 
extended  to  a  PCT  model.  The  F-16  plant  equations,  already  in 
transfer  function  form,  are  outlined.  Finally,  the  response 
models  for  both  the  KC-135  and  the  F-16  are  discussed. 


II-2  Generic  SISO  Model 

A  general  continuous  time  SISO  system  is  shown  in 
Figure  2-1  where  P  represents  the  plant,  D  is  a  disturbance 
to  the  plant,  G  is  the  fixed  compensator  and  F  is  a  prefilter 
to  the  system.  Since  it  is  desired  to  implement  G  as  a 

I 

digital  controller,  the  continuous  signal  E(s)  must  be  sampled 
to  generate  a  suitable  digital  signal  to  present  to  the 
•  controller.  The  output  of  a  digital  controller  is  a  discrete 

-V'  signal  and  therefore  must  be  transformed  into  a  continuous 
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d(t) 


r  (t)  F 

o  » 


Fig.  2-1.  Continuous  Time  System  Model 


signal  to  drive  the  plant  P.  This  is  because  in  general,  it 
is  undesirable  to  apply  a  signal  in  sampled  form  to  a  plant, 
because  of  the  high  frequency  components  inherently  present 
in  the  sampled  signal  (9:103).  This  process  is  performed 
by  a  hold  device.  The  simplest  data  hold  device  (requires 
no  memory)  and  by  far  the  most  common,  is  the  zero-order 
hold  (ZOH).  This  type  of  hold  device  is  assumed  throughout 
this  thesis  investigation.  Figure  2-2  shows  a  representation 
of  the  sampling  and  data  hold  operations  where  the  ZOH 
transfer  function  is  defined  by 


(2-1) 


It  should  be  noted  that  individually  the  sampler  in  Figure 
2-2  does  not  model  a  physical  sampler  nor  does  the  block 
model  a  physical  hold  device.  However,  it  can  be  shown 
mathematically  (8:71)  that  the  combination  does  accurately 
model  a  physical  sampler/data-hold  device. 

The  sampler  and  ZOH  are  included  in  the  block  diagram 
shown  in  Figure  2-3  to  generate  the  sampled-data  system  model. 


12 


S 


L 


vv 


Note  that  now  a  digital  controller,  G(z),  is  used  as  the 
compensator.  The  output  of  the  digital  controller  block  G(z) 
is  a  discrete  signal.  This  is  emphasized  in  some  digital 
control  system  texts  by  "inserting"  a  fictitious  sampler 
after  the  digital  controller.  This  fictitious  sampler  is  not 
shown  in  the  diagrams  contained  in  this  study.  The  prefilter 
F  is  implemented  as  an  analog  compensator.  The  reason  for 
this  is  that  generally,  QFT  design  technique  generates  a  very 
simple  pref liter.  The  prefilter  is  usually  either  a  first  or 
second  order  compensator.  There  is  also  a  distinct  advantage 


Fig.  2-3.  Sampled-Data  System  Model 

to  implementing  the  prefilter  with  an  analog  circuit.  The 
prefilters  are  essentially  low  pass  devices  and  serve  effec¬ 
tively  as  an  antialiasing  device  where  input  noise  above  the 
prefilter  breakpoint  is  attenuated.  The  prefilters  for  both 
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aircraft  are  shown  in  Appendix  E. 


An  approximation  to  the  sampled-data  system  of  Figure 
2-3  can  be  made  using  the  PCT  model  (9:248-251).  In  the  PCT 
system  model  the  ideal  sampler  of  Figure  2-3  is  approximated 
by  the  first  block  of  Figure  2-4.  This  approximation  is 
appropriate  because  the  sampling  process  modifies  the  funda¬ 
mental  frequency  and  all  its  harmonics  by  a  factor  of 
(9:84)  . 


Swpler  ZOH 

Approxijnation  Approximation 
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Fig.  2-4.  Sampler  and  ZOH  Approximations 


The  ZOH  is  approximated  by  the  first-order  Pade  approxi¬ 
mation  as  follows: 


2T 

Ts  +  2 


G  ( s ) 
pa 


(2-2) 


The  Pade  approximation  is  satisfactory  only  for  small 
values  of  T  (9:251).  Specifically,  the  approximation  is 
good  for  <j_  <  a»„/10.  The  PCT  system  model  is  shown  in 

C  a 

Figure  2-5  with  the  sampler  and  ZOH  approximations  included. 
Therefore  the  sampler  and  ZOH  devices  contained  in  the 
sampled  data  system  are  approximated  in  the  PCT  system  by 
tne  transfer  function 


Fig.  2-5.  PCT  System  Model 


G^(S) 


Ts  +  2 


(2-3) 


The  usefulness  of  the  PCT  system  model  and  Equation  (2-3) 
are  discussed  in  Chapter  3. 


II-3  KC-135  Lateral  Mode  Models 

The  lateral  mode  equations  of  the  KC-135  aircraft 
are  reduced  in  reference  1  to  a  linear  two  input-two  output 
model.  Those  equations  were  used  to  design  the  compensator 
for  the  lateral  aircraft.  A  general  block  diagram  correspond 
ing  to  the  lateral  aircraft  is  shown  in  Figure  2-6  where  F,  G 
and  P  are  the  prefilter,  fixed  compensator  and  plant  matrices 
respectively. 


KC-135  General  Lateral  Mode  System  Mode’ 


A  more  detailed  continuous  time  siqnal  flow  diagram  is  shown 
in  Figure  2-7.  Note  that  both  the  prefilter  and  fixed  compen 
sator  matrices  are  diagonal.  QFT  is  not  restricted  to 
diagonal  matrices,  but  the  designs  in  both  references  1  and  2 
used  diagonal  F  and  G.  The  inputs  to  the  2x2  model  are;  r^^ 
(bank  angle  command)  =30;  Z2  (sideslip  angle  command)  =  5. 
It  is  desired  that  the  output  track  the  input  r^^  and  that 
y2  be  ideally  zero  (with  r2  =  0).  When  r2  is  commanded  it  is 
desired  that  y2  track  r2  while  y^^  should  be  ideally  zero 
(with  r^  =  0 )  . 

The  sampled-data  lateral  system  model  is  shown  in 
Figure  2-8.  As  described  in  Section  II-2,  discrete  rather 
than  continuous  controllers  are  used.  The  PCT  approxima¬ 
tion  to  the  sampled-data  model  is  shown  in  Figure  2-9. 

-1 


Fig.  2-7.  kc-135  Lateral  Mode  Continuous  Time 

System  Model 


The  lateral  aircraft  equations  given  in  reference  1  are 
in  input-output  matrix  form,  where  an  input  matrix  N  premul¬ 
tiplies  the  input  vector  and  an  output  matrix  M  premultiplies 


The  input  vector  variables  are  control  wheel  deflection 
and  rudder  deflection  respectively.  The  output  vector 
variables  are  roll  angle  and  sideslip  angle  respectively. 
For  flight  condition  (F.C.)  #1  the  lateral  mode  equations 
in  input-output  matrix  form  are  (1:36): 


(S^+.74028S-. 010714) 


( .24789S+4. 46897) 


(  .039665s2+.0614  2S+.006  5)  -( S^+ . 227437S+1 . 43226 ) 


'.364455 

.708183 

(  .00074615S+. 020812) 

-( .026647S+1. 01677) 

w 


(2-4 


The  corresponding  input-output  matrix  equations  for  F.C.'s 
#2  and  #3  are  contained  in  Appendix  C. 

As  can  be  seen  in  Figure  2-8,  the  individual  plant 
transfer  functions,  relating  specific  inputs  to  outputs, 
are  needed  for  the  model.  Equation  (2-4)  can  be  written 
in  general  terms  as: 


flY  =  Nd 


(2-5 


Y  =  M 


(2-6 


Using  Equation  (2-6),  the  plant  transfer  function  matrix 
can  be  defined  as: 


P  = 


(2-7 


The  matrix  P  contains  the  individual  transfer  functions 
needed  for  the  system  model  used  in  this  study.  To  obtain 
the  plant  transfer  function  matrix,  a  matrix  inverse  and 


multiplication  must  be  done.  Since  the  individual  elements 
of  the  input  and  output  matrix  are  polynomials  in  S  it  is 
apparent  that  a  great  deal  of  effort  would  be  required  to 
obtain  the  transfer  function  matrix  P,  especially  if  the 
order  of  the  plant  matrix  is  greater  than  2x2.  In  addition 
to  being  a  very  tedious  task,  the  large  amount  of  algebraic 
manipulations  involved  would  invariably  result  in  many  errors 
The  continuous  domain  design  effort  for  the  KC-135  required 
similar  manipulations  but  instead  of  obtaining  an  exact 
solution,  a  computer  program  was  utilized  to  numerically 
generate  the  magnitude  and  phase  angle,  as  a  function  of 
frequency,  of  the  individual  functions.  This  method  is 
undesirable  since  the  transfer  functions  are  not  known 
explicitly.  Using  this  method  one  must  interpret  the  numeri¬ 
cal  magnitude  and  phase  data.  Although  not  necessary,  it  is 
usually  convenient  to  then  synthesize  an  appropriate  transfer 
function.  If  non-minimum  phase  terms  are  present  in  the 
plant,  it  is  very  easy  to  misinterpret  the  data  (10:272). 

To  avoid  this  problem  altogether,  a  computer  program 
called  MACSYMA  (11)  that  performs  symbolic  mathematical 
manipulations  is  used  to  obtain  the  transfer  function  matri¬ 
ces.  This  program  is  chosen  over  other  similar  programs 
because  the  computer  it  is  hosted  on  (Digital  Equipment 
Corp.  PDP-10)  has  a  large  amount  of  core  memory  (as  compared 
to  various  microcomputer  hosted  programs).  Since  the  algebra 
ic  manipulations  are  symbolic  rather  than  numerical,  the 
results  are  in  terms  of  the  variable  S  and  therefore  explicit 


transfer  functions  are  obtained.  Within  the  program  there  are 
functions  that  allow  the  inversion  and  multiplication  of 
matrices  whose  elements  are  polynomials.  It  should  also  be 
noted  that  the  symbolic  algebraic  manipulations  are  very 
accurate  since  intermediate  operations  are  done  with  rational 
numbers  instead  of  real  numbers.  The  resulting  plant  transfer 
functions  corresponding  to  Equation  (2-4),  F.C.  #1,  are: 


.364635(3+. 125 308+1. 2926 3) 
(S+.0104168) (S+. 8021 30) (S+. 0825017+ j 1.27795) 

_ .701577(3+2.39768X3-2.09826) _ 

( 3+. 0104168 ) (3+. 802130 ) (3+. 0825017+ jl. 27795) 

-.00074615 (3+. 290 305) (3+10.1471) (3-1.17926) 
(3+.0104168) (3+. 802130 ) (3+. 0825017+ jl. 27795) 


.026647(3+. 770138) (3+39.189122) (3-. 007821 11) 
(3+.0104168) (3+. 8021 30) (S+. 0825017+ j 1.27795) 


The  plant  transfer  functions  for  F.C.'s  #2  and  #3  ate  contain¬ 
ed  in  Appendix  C. 


II-4  KC-135  Longitudinal  Mode  Models 

Rigid  Body 

The  longitudinal  mode  equations  were  developed  in 
reference  1  as  a  three  input-three  output  model.  Only  one 
maneuver  is  done  in  the  longitudinal  direction  and  it  is  a 
four  degree  pitch  pointing  maneuver.  The  pitch  angle  output, 
'{ 2  should  track  the  pitch  angle  command  input,  1.2.  The  pitch 


pointing  maneuver  is  to  be  done  with  minimum  response  in 


commanded  maneuver  for  the  longitudinal  mode,  £  has  only  one 
element,  (22-  This  implies  that  the  inputs  to  loop  one  and 
three  are  zero. 

The  sampled-data  longitudinal  system  model  is  shown  in 
Figure  2-12  in  block  diagram  form.  Again,  as  in  the  lateral 
sampled-data  model,  the  compensation  is  performed  by  digital 
controllers  driven  by  sampled-signals .  The  ZOH's  provide  a 
piecewise  continuous  signal  to  the  plants,  Pj^j.  The  diagram 
of  the  PCT  approximation  to  the  sampled-data  system  is  the 
same  as  Figure  2-12  except  that  the  sampler  and  ZOH  are 
replaced  by  their  continuous  domain  approximations,  as 
discussed  in  Section  II-2,  and  the  continuous  domain  compensa 
tor  transfer  function  is  used  in  the  PCT  model. 

The  longitudinal  mode  aircraft  equations  given  in 
reference  1  are  in  input-output  matrix  form,  as  shown  in 
Equation  (2-5).  The  input  and  output  vectors  are  respec¬ 
tively: 


The  input  variables  are  defined  respectively  as  elevator 
deflection,  speed  brake  deflection  and  percent  engine  thrust 
(change  from  trim).  The  input-output  matrices  for  F.C.  #1  ar 
shown  in  Equation  (2-12)  (1:69). 


.  0  4 1 1 .9 

( .534779-1 . 096 ) 

(9+.003) 
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(-.07695^-. 04179) 

(  .0079+ . 5441  ) 

.0049 
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(9^+. 75379+2. 797) 
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.011619  -.0827  .0495 

-  -.021319  .032  0 

-1.649  .1733  0 

The  input-output  matrices  for  F.C.'s  2  and  3  are  given  in 
Appendix  C. 

The  symbolic  mathematical  manipulation  program  was  again 
used  to  perform  the  inversion  and  multiplication  operations 
indicated  in  Equation  (2-7)  to  yield  the  transfer  function 
matrix  P.  The  individual  transfer  functions  for  all  three 
F.C.'s  are  listed  in  Appendix  C. 

Plastic  Body 

A  continuous  domain  controller  design  was  also  done  for 
a  non-rigid  KC-135  model  where  the  first  and  second  body 
bending  modes  were  included  (1:89).  Due  the  small  magnitudes 
of  these  two  modes,  the  continuous  domain  QFT  design  process 
generated  the  same  compensators  as  for  the  rigid  body  (1:97). 
To  enable  hybrid  simulations  to  be  performed  on  the  elastic 
body,  the  elastic  body  transfer  functions  must  be  obtained. 

The  following  discussion  explains  this  process. 

The  form  of  the  rigid  aircraft  equation  is 


e 

<5sb 

6^ 


(2-12) 


=  Eul 


(2-13) 


output  matrices  respectively.  The  elastic  body  input-output 
matrices  are  included  in  Appendix  C.  Equations  (2-13)  and  (2- 
14)  can  be  written  in  the  following  form 

1  =  (2-15) 

and  Y  =  j  (2-16) 

where  the  input  and  output  vectors  are  defined  as  indi¬ 
cated  for  the  rigid  body  case  above.  Combining  the  rigid 
and  bending  mode  equations  yields 

Therefore  =  Er  +  Ze  (2-17) 

As  before,  the  manipulation  needed  in  Equation  (2-16)  is 

performed  using  the  symbolic  math  program.  The  results  of 

this  manipulation  showed  that  all  elements  of  P  are  zero 

— e 

except  for  and  Pg2i*  Therefore  the  only  difference 

between  P„  and  P^  are  the  elemental  transfer  functions  P„, , 
and  Pj,2i'  These  two  elemental  transfer  functions  are  obtained 
using  Equation  (2-17).  The  transfer  functions  ^c21' 

^ell  ^e21  in  Appendix  C. 

The  elastic  system  model  diagrams  are  the  same  as 
those  for  the  longitudinal  rigid  body  except  the  appropri¬ 
ate  transfer  functions  are  substituted  as  discussed  above. 

II-5  F-16  Models 

The  F-16  model  is  significantly  more  complicated 


than  the  KC-135  models  due  to  the  reconfiguration  nature  of 
the  original  design.  The  need  to  maintain  stability  and 
control  of  the  aircraft  in  the  event  of  single  or  simultan¬ 
eous  control  surface  failures  requires  cross  coupling  of 
the  compensation  signals.  Because  of  the  added  complexity, 
a  more  detailed  discussion  of  the  model  is  included  for  the 
F-16.  The  following  two  paragraphs  are  taken  from  reference 
2  with  some  minor  modifications. 

The  structure  for  the  continuous  time  domain  system 
model  for  the  F-16  is  shown  in  Figure  2-13  (2:16)  in  signal 
flow  graph  form.  A  two-elevator,  two-flaperon  system  is 
represented  with  each  individual  surface  separately  control¬ 
lable.  The  elevators  and  flaperons  are  capable  of  moving 
individually,  and  symmetrically  or  differentially  in  pairs. 

The  right  and  left  elevators  are  designated  by  and 
respectively  and  the  flaperons  are  designated  similarly  by 
and  ^2"  The  longitudinally  controlled  output  variable  is 
pitch  rate,  q,  while  the  laterally  controlled  output  variable 
is  roll  rate,  p.  The  plant  transfer  function  designator, 
relates  the  i  output  to  the  j  input  (control  variable).  A 
superscript  "a"  or  "b"  denotes  the  right  or  left  control.  For 
example,  denotes  the  effect  of  the  right  elevator  (6®) 
on  roll  rate  (output  y2 ) .  For  a  normally  operating  system, 

=  P^j.  Minus  signs  are  associated  with  P2j^  and  P22  so 
that  the  two  pairs  of  surfaces  move  differentially  to  generate 
roll  rate.  In  the  no  fail  case,  surfaces  move  an  equal 
amount  in  either  direction  to  generate  the  commanded  rates. 


parameter  changes  resulting  from  single  or  multiple  control 
system  component  failure.  In  Figure  2-13,  ^^2  and  are 
constants  which  determine  the  division  of  control  effort 
between  the  flaperons  and  elevators.  In  the  preliminary 
reconfiguration  study  done  in  reference  2,  optimization  of 
these  constants  was  not  addressed,  and  both  n' s  were  set  equal 
to  0.25.  Henceforth,  the  subscripts  are  dropped  on  these  two 
constants . 

Figure  2-14  shows  the  sampled-data  system  model  for  the 
F-16  in  block  diagram  form.  As  for  the  KC-135,  the  sampled- 


data  system  model  is  generated  by  replacing  the  continuous 
domain  compensator  with  the  sampler/digital  controller/ZOH 
combination.  The  PCT  system  diagram  is  the  same  as  Figure 
2-14  except  that  the  sampler  and  ZOH  are  replaced  by  their 
continuous  time  domain  approximations,  as  discussed  in  Section 
II-2,  and  the  continuous  domain  compensator  is  used  in  place 
of  the  digital  controller. 

A  set  of  equivalent  plant  transfer  functions  were 
developed  in  reference  2  to  populate  the  2x2  multiple  input- 
multiple  output  (MIMO)  system  model  discussed  above  for  four 
different  F.C.'s.  Each  set  of  transfer  functions  model  the 
healthy  (no  failures)  aircraft.  The  Pj^j's  for  all  four  F.C.'s 
are  given  in  Appendix  C  (2:13). 

II-6  Response  Models 

As  discussed  in  Appendix  A,  the  QFT  design  technique 
involves  quantifying  the  range  of  desired  responses  for  a 
given  command  input.  The  response  range  for  tracking  is 
expressed  as  an  upper  and  lower  limit  while  the  disturbance 
rejection  response  often  is  just  an  upper  bound.  The  desired 
responses  are  usually  given  as  time  domain  specifications, 
such  as  the  figures  of  merit  t„,  t„,  and  K„  and  are  based 
on  a  step  or  ramped  step  input.  These  desired  responses  are 
used  to  determine  bounded  regions  in  the  frequency  domain  as 
part  of  the  design  process  and  to  check  the  system  perfor¬ 
mance.  The  figures  of  merit  for  the  KC-135  and  the  F-16  QFT 
designs  were  established  as  part  of  the  continuous  domain 
design  process  of  reference  1  and  2.  The  performance  obtained 


by  the  sampled-data  system  designed  in  this  thesis  is  compared 
against  the  performance  of  the  continuous  system  and  its 
associated  response  models.  The  response  models  for  both 
aircraft  are  contained  in  Appendix  D  (1:211;  2:152).  Only  the 
figures  of  merit  are  given.  The  response  model  transfer 
functions  can  be  found  in  reference  1  and  2.  A  detailed 
description  of  how  the  response  models  are  specified  and  used 
is  contained  in  Appendix  A. 

II-7  Summary 

The  system  models  for  the  KC-135  lateral  and  longitu¬ 
dinal  modes  and  the  F-16  are  presented  in  this  chapter.  The 
relationship  between  the  continuous  domain,  sampled-data  and 
PCT  systems  models  are  discussed.  Plant  transfer  functions 
for  both  aircraft  are  given.  Response  models  for  the  KC-135 
and  F-16  are  presented. 


III.  Theoretical  Digital  Design 


III-l  Introduction 

This  chapter  is  devoted  to  the  design  of  robust  discrete 
controllers  for  the  KC-135  and  P-16  aircraft.  Since  the 
intent  of  this  thesis  investigation  is  to  determine  the 
feasibility  of  extending  continuous  domain  QFT  designs  to 
discrete  designs  while  still  maintaining  acceptable  perfor¬ 
mance  and  robustness,  the  original  continuous  domain  designs 
provide  the  starting  point  for  this  study.  The  first  aspect 
of  the  digital  design  is  to  examine  a  technique  for  compensat¬ 
ing  for  one  of  the  sources  of  lag  in  a  sample-data  system. 
Sampling  rates  for  both  aircraft  are  discussed  and  selected. 
Once  the  sampling  rates  (i.e.  frequency)  have  been  estab¬ 
lished,  various  transformations  are  performed  on  the  modified 
continuous  domain  QFT  controllers  to  generate  the  discrete 
controllers.  Both  the  standard  Z  and  Tustin  transformations 
are  examined.  The  last  portion  of  this  chapter  discusses  the 
problems  encountered  in  extending  the  F-16  continuous  domain 
design  to  the  discrete  domain. 

III-2  Counteracting  ZOH  Lag 

QFT  design  technique  is  usually  applied  in  the  contin¬ 
uous  frequency  domain  (although  QFT  can  be  applied  in  the 
discrete  frequency  domain).  As  discussed  in  Appendices  A  and 
B,  desired  response  models  are  used  to  develop  bounds  in  the 
continuous  frequency  domain  which  are  constraints  on  the  loop 
transmission  function,  L,  for  each  channel.  These  bounds  are 


usually  placed  on  a  Nichols  chart  to  show  by  graphical  means 
the  constraints  or  requirements  on  the  loop  transmission. 
Figure  3-1  shows  an  example  of  a  loop  transmission  shaped  to 
meet  bounds  at  various  frequencies.  In  fact  this  is  the 
continuous  domain  nominal  loop  transmission  designed  for  loop 
two  (tracking  loop)  of  the  KC-135  longitudinal  mode  (1:81). 

In  QFT,  an  "optimum*  loop  transmission  design  is  one  for 
which  Ii..,(jtJ^)  lies  just  on  the  boundary  B(j(j.  ).  In  this 
case,  optimum  means  that  the  desired  performance  and  robust¬ 
ness  are  guaranteed  while  obtaining  the  smallest  bandwidth 
possible.  As  can  be  seen  from  Figure  3-1,  the  loop  transmis¬ 
sion  tracks  down  the  universal  high  frequency  boundary  (UHFB) 
with  only  a  few  degrees  of  separation  over  a  frequency  range 
of  about  75  rad/sec.  Remembering  that 

L  =  PG  (3-1) 

it  is  apparent  that  any  change  in  the  compensation,  G, 
which  results  in  additional  lag  could  take  the  loop  transmis¬ 
sion  into  the  UHFB.  If  this  occurs,  and  depending  on  the 
severity  of  the  additional  lag,  some  plants  in  the  parameter 
space  can  exhibit  less  than  desired  performance  or  they  can  be 
unstable  (5:291). 

As  discussed  in  Section  II-2,  a  ZOH  is  included  in  the 
sampled-data  system  models  to  construct  a  piecewise  continuous 
signal  from  the  discrete  output  of  the  digital  controllers. 
Figure  3-2  shows  the  frequency  response  characteristics  of  the 


Deg 

Fig.  3-1,  Example  of  Bounds  and  Nominal  Loop  Transmission 

ZOH  where  the  ZOH  transfer  function  is  shown  in  Equation 
(2-1).  It  is  clear  from  the  magnitude  curve  of  Figure  3-2 
that  the  ZOH  approximates  the  behavior  of  a  low-pass  filter. 
This  is  desirable,  as  previously  mentioned,  because  the  low- 
pass  characteristic  tends  to  smooth  the  high  frequency  noise 
in  the  system.  The  phase  curve  of  Figure  3-2  shows  that  the 
ZOH  has  a  considerable  lag  characteristic.  Since  the  intent 
in  the  sampled-data  system  is  to  maintain  the  frequency 


Fig.  3-2.  Gain  and  Phase  Characteristics  of  a  ZOH 

response  characteristics  of  the  original  continuous  domain 
loop  transmission  (and  therefore  compensators),  it  is  clear 
then  that  the  ZOH  by  itself  can  have  an  undesirable  effect. 
Also,  by  comparing  the  frequency  spectra  of  a  continuous 
signal,  E(j<J),  and  the  associated  output  of  an  ideal  sampler, 

It 

E  { j(J) ,  it  can  be  seen  that  the  output  magnitude  has  been 
changed  by  a  factor  of  (9;85). 

Thus  it  is  desirable  to  modify  the  original  compensa¬ 
tors  in  a  fairly  straightforward  manner  whereby  the  undesir¬ 
able  effects  of  the  sample  and  hold  processes  are  diminished. 
In  Section  II-2,  the  PCT  system  model  is  developed  as  a  means 
of  approximating  the  effects,  in  the  S  domain,  of  the  sample 
and  hold  processes.  Equation  (2-3)  shows  the  transfer 


function,  G^{s),  which  approx inates  these  two  processes.  The 

I 

original  S  domain  QFT  compensators,  G  (s),  can  be  modified  by 
[G^( s ) ]  ^  to  obtain 

G(s)  =  [G^(r) ]'^g' (s)  (3-2) 

The  net  affect  is  that  the  gain  is  modified  by  a  factor  of 
T  and  phase  lead  is  added  to  the  compensation.  The  overall 
effect  on  the  loop  transmission  can  be  shown  mathematically  in 
the  following  way.  Assuming  the  PCT  approximation  to  the 
sampled-data  system  model,  the  loop  transmission  can  be 
written  as 

-I 

L  =  [g^(s3g'  (s)P(s)  (3-3) 

I 

If  the  original  S  domain  compensators  G  (s)  are  replaced 
by  the  modified  compensators  G(s)  using  Equation  (3-2) 
it  can  be  seen  that  the  original  loop  transmission  of  Equation 
(3-1)  results.  All  of  the  S  domain  QFT  compensators  are 
modified  as  shown  in  Equation  (3-2)  before  being  transformed 
to  their  discrete  form.  The  transformations  are  discussed 
in  Section  III-4. 

III-3  Sample  Rate  Selection 

There  are  a  multitude  of  techniques  used  to  select  the 
sample  rate  of  a  sampled-data  system,  many  of  which  often 
depend  on  the  original  design  technique  used.  A  review  of  the 
literature  (12:975)  shows  that  in  general,  control  system 
designs  done  first  in  continuous  domain  and  then  extended  to 
the  discrete  domain,  tend  to  require  a  higher  sample  rate  than 


those  systems  designed  originally  in  the  discrete  domain 
(direct  method).  The  major  considerations  which  influence  the 
selection  of  sample  rate  can  be  divided  into  three  broad 
categories;  closed-loop  bandwidth  or  time  response  require¬ 
ments,  rejection  of  unwanted  disturbances  and  sensitivity  to 
parameter  variations  (12:975).  These  categories  tend  to 
provide  bounds  only  on  the  lower  limit  of  the  sample  rate. 
Generally,  for  an  aircraft  digital  controller,  the  selection 
of  a  sample  rate  represents  a  trade  off  between  cost/accuracy 
and  desirable  characteristics  in  one  of  the  three  areas  listed 
above.  By  lowering  the  sampling  frequency,  the  A/D  and 

D/A  equipment  cost  is  reduced  directly.  Also,  the  increased 
accuracy  obtained  by  slower  sampling  (increased  T)  is  well 
documented  (12:975)  and  can  be  reflected  into  reduced  costs  by 
decreasing  the  microprocessor  and  A/D  wordl.ngth .  Since  cost 
is  not  a  direct  concern  of  this  study,  an  upper  bound  on  the 
sampling  rate  is  established  by  what  sample  rates  are  current¬ 
ly  considered  practical  and  realizable  in  avionics  systems. 
Sample  rates  of  40  to  60  Hz  are  used  as  a  realistic  upper 
bound.  Only  the  first  category  of  the  three  mentioned  above 
is  considered  in  the  selection  of  the  value  of  T  in  this 
thesis  since  the  last  two  categories  are  explicitly  included 
in  the  QFT  design  process.  These  two  categories  are  usually 
reflected  in  the  bandwidth  of  the  system  anyway  (QFT  designs). 

An  absolute  lower  bound  on  the  sample  rate  is  estab¬ 
lished  by  a  desire  to  track  a  specified  command  input  signal. 
This  bound  derives  its  theoretical  basis  from  Shannon's 


Sampling  Theorem,  The  theorem  states  that  in  order  to 
reconstruct  a  band-limited  continuous  signal  from  samples  of 
that  signal,  one  must  use  a  sample  rate  which  is  at  least 
twice  as  high  as  the  highest  frequency  contained  in  the 
signal.  It  should  be  added  that  only  signals  of  "significant 
amplitudes*  need  be  considered  (8:79).  The  implication  here 
is  that  the  signals  are  from  the  closed-loop  system.  This 
theorem,  however,  only  provides  a  fundamental  lower  limit  and 
in  practice  it  tends  to  be  significantly  higher.  Franklin  and 
Powell  provide  a  "rule  of  thumb"  that  the  sample  rate  must  be 
about  4  to  20  times  the  closed  loop  bandwidth  of  the  system  to 
obtain  the  desired  transient  characteristics  (13:276).  This 
rule  is  used  as  a  "guide-post". 

The  closed-loop  bandwidth  of  both  aircraft  can  be 
estimated  from  the  Nichols  charts  used  in  the  designs  of 
references  1  and  2.  For  the  KC-135,  the  largest  closed-loop 
bandwidth  was  approximately  55  rad/sec  which  was  associated 
with  the  longitudinal  mode  tracking  channel  (loop  2).  This 
bandwidth  is  quite  high.  The  other  loops  were  generally  well 
below  30  rad/sec.  The  design  of  the  loop  transmission 
corresponding  to  this  maximum  bandwidth  (1:86)  shows  that 
there  is  some  overdesign,  which  translates  into  a  higher 
than  necessary  bandwidth.  Since  the  initial  selection  of 
sample  rate  is  based  upon  bandwidth  it  is  important  to  stress 
that  the  designer  of  QFT  continuous  domain  compensators 
insure  that  the  bandwidth  be  as  small  as  possible.  As  a  first 
trial  for  both  the  lateral  and  longitudinal  modes  of  the  KC- 


135,  a  sample  rate  of  40  Hz  is  used  which  is  about  five  times 
the  maximum  closed-loop  bandwidth.  As  is  discussed  in  the 
next  chapter,  this  sample  rate  proves  to  be  quite  satis¬ 
factory  . 

The  maximum  closed-loop  bandwidth  for  the  F-16  contin¬ 
uous  domain  design  was  34  rad/sec  (2:37).  A  sampling  rate  of 
40  Hz,  which  represents  a  multiplier  of  about  7,  is  also 
chosen  for  the  F-16,  Intuitively,  one  can  expect  the  F-16  to 
require  a  higher  sampling  rate  due  to  the  naturally  faster 
response  of  the  aircraft,  its  relative  size,  and  the  reconfig¬ 
uration  requirements  on  the  QFT  design.  In  part,  this  can  be 
explained  by  the  fact  that  the  designs  were  done  by  two 
different  designers  who  attained  different  levels  of  profi¬ 
ciency  in  the  QFT  design  technique.  Although  not  accomplished 
in  this  study,  it  is  quite  possible  that  a  sampling  rate  of 

less  than  40  Hz  would  have  provided  good  performance  for  the 

KC-135. 

The  F-16  proved  to  be  unstable  for  all  four  F.C.'s  at 

the  40  Hz  sample  rate  initially  chosen.  The  sample  rate  was 

then  increased  to  60  Hz  which  generated  satisfactory  perfor¬ 
mance  for  one  of  four  F.C.'s  (#1).  Some  modifications  to  the 
Z  domain  compensators  were  made  which  provided  some  additional 
improvement  in  performance.  The  modifications  made  to  the 
controllers  and  explanation  of  the  poor  performance  are 
discussed  in  the  last  section  of  this  chapter. 

The  continuous  domain  QFT  compensators  for  the  lateral 
mode  of  the  KC-135  are  shown  in  Equations  (3-5)  and  (3-6) 
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after  having  been  modified  as  discussed  in  Section  III-2 
unmodified  compensators  can  be  found  in  reference  1.  As 
discussed  previously,  the  sample  rate  chosen  was  40  Hz. 


1081075(5+4) (S+4) (S+17) (S+50) (S+80) 


(S+10) (S+10) (S+83) (S+100) (S+60+j80) 
295125(3+7) (S+80) 


(S+50 ) (S+96+ jl28 ) 


At  the  beginning  of  this  thesis  effort  a  gross  overdesign 
was  discovered  in  the  loop  transmission  (loop  1)  of  the  KC- 
135  lateral  mode.  This  was  due  to  an  error  in  plotting 
the  UHFB  on  the  Nichols  chart.  After  reference  1  was  finish¬ 
ed,  the  author  (Capt  Russell)  removed  some  of  the  overdesign 
by  adjusting  the  location  of  the  complex  pole  pair  of  Gj^(s). 
The  previous  location  was  at  s  =  -900+jl200.  As  can  be  seen, 
this  was  a  significant  modification  and  explains  the  differ¬ 
ence  between  Equation  (3-5)  and  Gj^(s)  in  reference  1.  The 
modified  S  domain  QFT  compensators  for  the  KC-135  longitudinal 
mode  and  the  F-16  are  given  in  Appendix  E. 


III-4  Transformation  of  Compensators 

There  are  several  ways  of  obtaining  G(z)  from  G(s),  but 
only  two  are  considered  in  this  thesis.  As  mentioned  pre¬ 
viously,  these  two  methods  are  the  standard  Z  transform  and 
the  Tustin  transformation  (also  called  the  bilinear  transfor¬ 
mation  or  trapezoidal  integration) .  Some  of  the  other  methods 
are  forward  rectangular  integration  (forward  difference), 
backward  rectangular  integration  (backward  difference). 


matched  Z  transform  (pole-zero  mapping)  (8:338;  13:61)  and 
hold  equivalence  (13:62).  Of  the  other  methods,  the  matched  Z 
transform  and  the  hold  equivalence  seem  to  provide  the  best 
results  while  the  rectangular  integration  methods  generally 
give  poor  results  (13:66). 

III-4-1  Standard  ^  Transform.  The  standard  Z 

transform  is  chosen  primarily  because  it  encompasses  the 
definition  of  the  Z  transform.  The  Z  transform  can  be 
stated  in  its  most  simple  form  as  a  mapping  from  the  Laplace 
domain  to  the  Z  domain: 

z  =  e"^®  (3-6) 

To  transform  G(s)  to  G(z)  using  the  standard  Z  transform 
in  a  closed  form  one  must  first  obtain  a  partial  fraction 
expansion  of  G(s).  Generally,  the  poles  of  each  term  in 
the  resulting  summation  are  then  mapped  using  Equation 
(3-6)  and  the  residues  of  each  term  are  calculated  (9:113). 

The  standard  Z  transform  is  best  applied  using  a  computer- 
aided  design  program  such  as  TOTAL  (14),  especially  if  G(s)  is 
of  higher  order  than  two  and/or  has  repeated  poles.  The 
CAD  package  of  reference  14  is  used  to  obtain  G(z)  from  G(s) 
using  the  standard  Z  transform.  Simulations  using  these 
compensators  showed  that  the  system  is  unstable.  Closer 
examination  of  the  discrete  compensators  revealed  that  the 
transformation  generated  the  proper  poles  but  had  mapped  the 
zeros  to  the  same  location  as  the  poles.  Graphically  this  can 
be  seen  by  comparing  Figures  3-3  and  3-4.  Figure  3-3  is  the 


continuous  frequency  response  of  Equation  (3-4),  which  is 

G^(s)  for  the  lateral  mode  of  the  KC-135.  The  discrete 

frequency  response  of  0^(z)  is  shown  in  Fiqure  3-4.  As  can  be 

seen,  the  frequency  response  is  essentially  flat  in  the 

frequency  range  defined  by  the  primary  strip,  Wg/2,  The  Z 

transform  of  Equation  (3-5)  is  given  by  Equation  (3-7). 

Similar  results  are  obtained  on  all  the  other  compensators 

transformed  using  the  CAD  package  (14)  and  the  standard  Z 

transform.  The  desirable  gain  and  phase  characteristics  of 

the  compensators  have  therefore  been  lost  in  the  transforma¬ 
tion  . 


^  1081075Z(Z-.1306)(Z-.7798  +  j. 01498) 
^  (Z-.08208) (Z-. 1256) (Z-. 7788)^ 


(Z-.04187+j.2072) 

(Z-.09285+j.2029) 


Another  discrepancy  in  the  CAD  package  arose  using  the 
standard  Z  transform.  Equation  (3-8)  shows  Gj^(z)  for  the 
longitudinal  mode  of  the  KC-135  resulting  from  the  CAD 
program. 

.  -654.75Z(Z-1.115)  (3_qj 

^  (Z-.3824+j.3937) 

The  corresponding  continuous  domain  compensator  is  shown  in 
Appendix  E.  The  significant  point  here  is  that  the  CAD 
package  generated  a  compensator  with  a  zero  outside  the  unit 
circle  from  one  that  contained  only  left  half  plane  poles  and 
zeros.  Performing  the  standard  Z  transform  by  hand  produced 


Thus  it  is  apparent  that  the  CAD  package  is  not  producing  the 
correct  results.  In  addition,  examination  of  the  discrete 
frequency  response  of  Equation  (3-9)  exhibited  the  same  loss 
of  desirable  magnitude  and  phase  characteristics  as  compared 
with  the  analog  frequency  response.  In  light  of  the  above 
discussion,  the  standard  Z  transform  is  not  a  useful  tool  in 
extending  the  continuous  domain  compensators  to  discrete 
domain  and  is  not  discussed  further  in  this  study. 

III-4-2  Tustin  Transformation.  The  Tustin  transfor¬ 

mation  is  an  approximation  to  the  standard  Z  transform 
(8:194;  9;237;  13:56).  For  a  properly  chosen  sample  period, 

T,  the  Tustin  transformation  has  shown  to  provide  good 
correspondence  between  the  continuous  and  discrete  domain 
frequency  responses  and  thus  is  quite  popular  as  a  design 
tool.  The  definition  of  the  Tustin  transformation  is 

s  =  i.  LSLJl  (3-10) 

T  z  +  1 

The  transformation  is  applied  by  substituting  the  above 
relation  into  the  S  domain  compensator  transfer  functions 
with  a  specified  sample  period.  In  some  cases,  if  good 
discrete  response  characteristics  are  not  achieved,  it  is 
necessary  to  "prewarp"  the  S  domain  frequencies  prior  to 
using  Equation  (3-10)  (9:246).  The  Tustin  transformation  is 
applied,  at  least  initially,  without  prewarping. 

As  mentioned  previously,  the  accuracy  of  the  mapping  of 
the  Tustin  transformation  is  dependent  on  the  relationship 
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KC-135  Tustin  Mapping  Regions,  =  40  Hz 
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7.  F-16  Tustin  Mapping  Regions,  f  =  60  Hz 
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that  good  mapping  characteristics  can  be  expected  for  the 
majority  of  the  poles  and  zeros.  There  are,  however,  some 


exceptions.  Each  S  domain  QFT  compensator  has  a  high  frequen 
cy  complex  pole  pair  for  the  purpose  of  noise  rejection.  All 
of  the  pole  pairs  (both  aircraft)  except  for  G^(s)  and  G^Cs) 
of  the  KC-135  longitudinal  mode  are  warped  significantly. 

This  warping  has  no  real  effect  since  these  poles  were 
generally  at  or  well  above  the  sampling  frequency.  G^(s)  for 
the  KC-135  lateral  mode  had  several  real  poles  that  were 
warped  onto  the  left-half  Z  plane  real  axis  within  the  unit 
circle  thereby  severally  limiting  there  intended  lag  effect 
(on  the  loop  transmission).  As  is  shown  in  the  next  chapter, 
this  warping  had  no  effect  on  the  robustness  or  performance 
because  all  the  bounds  are  still  met.  In  other  words,  the 
loop  transmission  moved  to  the  right  due  to  the  loss  of  the 
lag  effect  but  still  met  all  the  bounds.  The  F-16  compensa¬ 
tors  each  had  a  zero  that  is  warped  significantly  and  the 
effect  of  the  warping  is  discussed  in  the  next  section. 

The  result  of  performing  the  Tustin  transformation  on 
the  KC-135  lateral  mode  compensators  is  shown  in  Equations 
(3-11)  and  (3-12). 


2490.75Z(Z+1) (Z-. 904762) ^(Z-. 649485) 
( Z-. 777778 ) ^( Z+. 018405 ) (Z+. 111111) 

(Z-. 230769) _ 

(Z+.138462+j. 492308) 

667. 252Z ( Z+1 ) (Z-. 839080) _ 

( Z-. 230769 ) (Z+. 405405+ j. 432432) 


An  example  of  the  good  mapping  characteristics  achieved  by 


the  Tustin  transformation  can  be  seen  in  Figure  3-8  where 
both  the  discrete  frequency  response  for  the  KC-135  lateral 
mode  Gj^(z)  and  the  corresponding  continuous  frequency  response 
are  plotted.  The  Tustin  transformed  discrete  controllers  for 
the  KC-135  longitudinal  mode  and  the  F-16  are  shown  in 
Appendix  E. 

It  should  be  noted  that  the  Tustin  transformation 
algorithm  in  the  CAD  program  (14)  performed  the  transformation 
correctly  for  the  relatively  low  order  compensators  associated 
with  the  KC-135.  The  CAD  program  was  inaccurate  for  the  high 
order  F-16  compensators.  The  algorithm  produced  poles  outside 
the  unit  circle  from  transfer  functions  that  contained  all 
left-half  S  plane  poles  and  also  generated  complex  pole  pairs 
from  real  S  plane  poles.  It  is  assumed  that  these  inaccura¬ 
cies  are  the  result  of  numerical  problems  within  the  compu- 


ter.  Due  to  this  numerical  problem,  the  Tustin  transformation 
for  the  F-16  compensators  were  obtained  by  hand  calculations. 

III-5  F-16  Pole-Zero  Mapping  Problems 

Simulations  for  the  F-16  sampled-data  system  showed 
that  for  F.C.'s  2,  3,  and  4,  the  system  is  unstable  when 
using  the  discrete  compensators  transformed  via  Tustin. 
Simulations  involving  F.C.  1  (all  CSC  modes)  are  stable  and 
met  all  specifications.  As  mentioned  in  the  previous  section, 
one  real  zero  in  each  of  the  F-16  compensators  is  warped  onto 
the  left-half  real  axis  within  the  unit  circle  at  sampling 
frequencies  of  40  and  50  Hz.  With  f^  =  60  Hz  the  zero 
associated  with  G2(z)  is  mapped  into  the  origin  of  the  Z  plane 
while  the  zero  associated  with  Gj^(z)  is  mapped  into  the  left- 
half  Z  plane  real  axis  within  the  unit  circle.  Table  3-1 
shows  the  migration  of  the  poles  and  zeros  of  Gj^(z)  as  the 
sampling  frequency  is  increased  from  40  Hz  to  200  Hz.  Table 
3-2  presents  similar  data  for  G2(z).  Zeros  number  8  and  7  are 
the  zeros  in  question  for  Gj^(z)  and  G2(z),  respectively. 

The  instability  is  produced  by  a  net  increase  in  phase 
lag  introduced  by  the  Tustin  transformation.  In  order  to 
alleviate  this  problem,  some  phase  lead  is  added  to  G2(z) 
directly  in  the  form  of  pole-zero  pair.  Simulations  are 
performed  on  F.C.  2  with  several  pole-zero  combinations.  With 
a  pole-zero  combination  of 

Z  -  .5 


added  to  ^2^^^  overall  compensator  gain  set  equal  to 

-2.0,  a  stable  system  is  achieved  which  met  specifications 
for  F.C.'s  1  and  2  although  F.C.'s  3  and  4  are  still  unstable. 

Poles  or  zeros  located  on  the  negative  real  axis  within 
the  Z  plane  unit  circle  have  only  a  small  effect  on  the 
system  response  (13:36).  Inspection  of  Table  3-2  reveals 
that  even  for  very  high  sample  rates  (200  Hz),  zero  number  8 
of  G,(z)  remains  in  the  left-half  Z  plane.  Essentially,  the 
phase  lead  associated  with  this  zero  in  the  S  plane  (S=-430) 
has  been  lost  through  the  transformation  and  the  UHFB  has 
been  crossed.  Even  though  this  zero  is  located  far  out  in  the 
left-half  S  plane,  it  contributes  a  phase  lead  of  six  degrees 
one  decade  before  the  corner  frequency  and  increases  from 
there.  Due  to  the  method  of  presentation  of  the  template  data 
in  reference  2,  it  is  not  possible  to  determine  which  point  on 
the  templates  is  generating  the  instabilities  for  F.C.'s  3  and 
4. 

It  is  hypothesized  that  the  addition  of  phase  lead  to 
G^(z)  in  the  form  of  a  pole-zero  combination  will  produce  a 
stable  system  and  satisfy  the  specifications,  A  satisfactory 
pole-zero  combination  will  of  course  require  some  amount  of 
"cut  and  try"  and  may  produce  stable  systems  for  F.C.'s  3  and 
4  but  at  the  same  time  make  F.C.'s  1  and/or  2  unstable. 

Instead  of  adding  phase  lead  to  the  discrete  compensators, 
better  results  from  the  Tustin  transformation  might  be 
obtained  by  first  prewarping  the  poles  and  zeros.  Of  course 
only  the  poles  and  zeros  which  show  significant  warping 


need  be  ptewarped.  Another  possible  way  to  obtain  better 
mapping  characteristics  is  to  use  the  matched  Z  transform 
where  both  the  poles  and  zeros  are  mapped  to  the  z  plane  using 
Equation  (3-6).  This  transform  never  maps  stable  real  poles 
or  zeros  to  the  negative  Z  axis  as  occurs  with  the  Tustin 
transformation.  The  application  of  the  systematic  approach  of 
QFT  design  technique  is  not  possible  in  the  Z  domain  without 
mapping  the  plant  templates  and  bounds  into  the  Z  plane  (or 
equivalent).  Due  to  time  limitations  and  the  above  considera¬ 
tions,  only  the  KC-135  aircraft  is  discussed  in  the  Implemen¬ 
tation  Phase  of  this  thesis.  The  time  responses  for  the  KC- 
135  and  F-16  continuous  and  sampled-data  systems  are  presented 
in  Chapter  IV. 


Table 
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Summary 

This  chapter  discusses  a  method  for  reducing  the  lag 
effect  of  the  ZOH  in  the  discrete  system.  This  method  is 
based  on  the  PCT  approximate  model  for  the  sampled-data 
system.  Considerations  in  the  selection  of  a  sample  rate 
are  outlined  and  the  sample  rates  for  the  two  aircraft  are 
chosen.  The  standard  Z  and  Tustin  transformations  are 
discussed  as  they  apply  to  the  continuous  QFT  designs.  It 
is  found  that  the  standard  Z  transform  is  not  a  practical 
tool  for  discretizing  the  QFT  analog  compensators.  The 
results  of  the  Tustin  transformation  are  presented.  The 
mapping  characteristics  of  the  Tustin  transformation,  in 
relation  to  the  QFT  continuous  domain  designs  are  also 
discussed.  Finally,  a  possible  explanation  for  the  poor 
performance  of  the  F-16  sampled-data  system  is  presented. 


Theoretical  Phase  Simulation  Results 


IV- 1  Introduction 

This  chapter  contains  the  KC-135  and  F-16  simulated 
responses  obtained  in  the  "theoretical  phase"  of  this  study. 
The  general  simulation  methodology  is  discussed  in  the  first 
section.  The  KC-135  lateral  mode  responses  are  then  present¬ 
ed.  The  KC-135  longitudinal  mode  simulation  results  for  the 
rigid  and  elastic  aircraft  are  presented  together  in  the  third 
section.  Finally,  the  simulation  results  for  the  F-16  are 
discussed . 

IV-2  Simulation  Methodology 

As  discussed  previously,  the  theoretical  phase  of  this 
thesis  considers  only  the  discrete  time  or  sampling  aspects 
of  the  digital  system.  Therefore,  the  simulations  of  this 
chapter  are  performed  using  the  same  computational  precision 
for  both  the  controllers  and  the  plants.  The  simulations  are 
performed  on  a  Digital  Equipment  Corporation  VAX  11/780 
digital  computer  using  double  precision  floating  point  (DPFP) 
arithmetic.  For  the  purposes  of  this  study,  DPFP  constitutes 
a  very  large  or  "infinite"  wordlength.  Chapters  v  and  VI 
examines  the  finite  word  length  case. 

The  original  continuous  domain  designs  of  references  1 
and  2  are  simulated  first.  This  is  intended  to  be  an  empiri¬ 
cal  validation/verification  step  and  is  done  for  several 
reasons.  First,  a  considerable  amount  of  manipulation  is 
performed  on  the  KC-135  plants,  as  discussed  in  Chapter  II,  to 


obtain  the  form  used  in  this  thesis.  A  good  comparison 
between  the  continuous  domain  results  is  considered  an 
implicit  validation  of  the  equivalent  plants.  The  setup  of 
the  simulations  in  this  thesis  is  essentially  along  the  lines 
of  the  system  block  diagrams  shown  in  Chapter  II  and  is 
slightly  different  than  the  simulation  structures  used  in 
references  1  and  2  (1:262;  2:163).  Also,  various  discrepan¬ 
cies  are  reported  in  reference  2  with  respect  to  simulation 
results  (2:58,60,73).  They  are  attributed  to  either  the  use 
of  compensator  approximations  or  the  need  for  an  adjustment  in 
the  prefilters.  Due  to  a  state  limitation  in  the  computer 
program  used  in  reference  2  (F-16),  the  orders  of  Gj^(s)  and 
G2(s)  were  reduced  from  10th  to  5th  and  from  9th  to  5th, 
respectively.  The  simulation  program  used  for  this  thesis  has 
no  state  limitations,  therefore  the  full  order  compensators 
are  used  in  the  simulations.  Several  of  the  discrepancies 
referred  to  above  did  not  appear  in  the  simulations  involving 
the  full  order  compensators,  as  can  be  seen  in  Appendix  F. 

A  computer  program  called  MATRIX^^  (7)  is  used  in  this 
study  to  perform  the  continuous  and  "hybrid"  simulations. 

This  program  has  a  wide  variety  of  MIMO  system  modeling  and 
simulation  capabilities.  The  actual  simulations  are  perform¬ 
ed  using  a  subset  set  of  the  MATRIX^^  software  called  System 
Build.  The  System  Build  capability  in  MATRIX^^  provides  an 
interactive,  menu-driven  graphical  environment  for  building, 
modifying,  and  editing  computer  simulation  models.  Any 
combination  of  linear,  non-linear,  continuous-time,  discrete- 


time,  or  multi-rate  models  that  describe  a  system  can  be 
constructed  from  a  library  of  basic  building  blocks  (7:SB 
P-1).  This  capability  replaces  the  FORTRAN  programing  and 
batch  operation  that  was  used  in  references  1  and  2. 

Both  the  continuous  and  sampled-data  simulations  are 
performed  using  the  system  block  diagram  structures  of  Chapter 
II.  For  the  continuous  time  simulations,  refer  to  Figures 
2-7,  2-11,  and  2-13,  respectively,  for  the  KC-135  lateral 
mode,  KC-135  longitudinal  mode  and  the  F-16.  Similarly,  refer 
to  Figures  2-8,  2-12  and  2-14  for  the  sampled-data  simula¬ 
tions.  Note  that  each  simulation  plot  in  this  thesis  contains 
all  three  F.C.'s  for  the  KC-135  or  all  four  F.C.'s  and  all  six 


CSC  modes  for  the  F-16  except  where  stated  otherwise.  This 
form  of  presentation  is  chosen  because  it  effectively  illus¬ 
trates  the  robust  characteristics  of  the  OFT  design  technique. 

Computational  time  delay  is  not  included  in  the  simula¬ 
tions  of  this  thesis  in  order  to  simplify  somewhat  the 
system  model.  Although  at  first  this  seems  to  be  a  gross 
simplification,  it  must  be  remembered  that  in  the  limited 
scope  of  a  thesis  study,  inclusion  of  delay  usually  addresses 
the  delay  associated  /ith  the  computation  of  the  difference 
equation  only.  This  does  not  include  the  delay  attributed  to 
the  "overhead*  tasks  of  executing  extensive  software  algo¬ 
rithms  for  system  error  checking  or  utility  functions.  In 
addition,  it  is  reasonable  to  expect  a  decrease  in  computa¬ 
tional  time  for  these  overhead  tasks  anyway,  since  robust  QFT 
controllers  might  relieve  the  need  for  gain  and/or  controller 
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scheduling.  As  microprocessors  (digital  avionics  computers) 
continually  improve  in  speed  and  computational  efficiency,  the 
computational  delay  associated  only  with  the  calculation  of 
the  difference  equations  becomes  less  significant. 

As  an  example,  consider  a  recent  study  done  (15:75) 
using  the  Texas  Instruments  TMS32010  general  digital  control¬ 
ler/digital  signal  processing  computer  (available  since  1983). 
The  TMS32010  can  execute  an  add  or  multiply  in  200  nanoseconds 
and  can  perform  a  combination  add,  multiply,  and  shift  in  400 
nanoseconds.  These  are  the  basic  operations  required  to 
solve  the  difference  equation.  It  was  shown  that,  using  a 
sample  frequency  of  40  Hz,  the  calculation  of  a  1000th  order 
digital  filter  (difference  eqn)  would  require  only  5%  of  the 
sample  period.  Thus  the  effect  would  be  fairly  negligible  for 
a  10th  order  controller.  Another  way  to  see  the  effect  is  to 
consider  the  definition  of  transport  lag 

e-sf  ^  g-jwf 

I 

where  t  is  the  computational  delay  time  associated  with  the 
difference  equation.  The  computational  time  for  a  second 
order  filter  requires  2.5  microseconds  (15:75),  therefore  for 

I 

a  10th  order  filter,  t  =  12.5  microseconds.  Inserting  this 

I 

t  into  Equation  (4-1)  and  considering  the  frequencies 
contained  in  the  system  with  a  sampling  frequency  of  40  Hz,  it 
is  clear  that  the  lag  effect  is  quite  small. 

The  standard  figures  of  merit  (10:92)  are  used  to 
assess  the  performance  of  the  continuous  and  sample-data 


systems.  The  specifications  for  each  system  are  shown  in 
Appendix  D.  Although  the  responses  are  plotted  for  10 
seconds,  they  are  examined  beyond  that  point  to  verify  the 
onset  of  steady-state  conditions.  it  should  be  noted  that  in 
some  cases  for  the  KC-135,  final  values  are  obtained  at 
times  beyond  10  seconds,  as  was  done  in  reference  1. 

In  the  simulation  plots  shown  in  this  thesis,  all 
F.C.  responses  due  to  the  same  input  command  are  shown  on  the 
same  plot.  In  addition,  for  the  longitudinal  mode  of  the 
KC-135  rigid  and  elastic  body  responses  are  also  shown  on  the 
same  plot.  Responses  at  all  F.C.s  and  for  all  CSC  modes  for 
the  F-16  are  shown  on  the  same  plot.  Exceptions  to  the  above 
statement  are  specifically  noted  where  applicable.  The 
individual  F.C.  (and/or  CSC  mode)  responses  are  not  labeled  on 
the  plots.  The  figures  of  merit  data  given  allow  identifica¬ 
tion  of  specific  responses  wherever  the  responses  are  signifi¬ 
cantly  different.  The  figure  label  for  each  sampled-data 
system  simulation  plot  (Chapter  IV  and  V)  contains  a  parenthe¬ 
tical  reference  to  the  corresponding  continuous  system  plot 
for  comparison  purposes. 

IV-3  KC-135  Lateral  Mode  Responses 

The  input  for  the  KC-135  lateral  mode  is  a  30  degree 
bank  angle  command  for  the  first  maneuver  and  a  5  degree 
sideslip  command  for  the  second  maneuver. 


Bank  Angle  Command  Response  -  Continuous  System 

The  bank  angle  command  responses  for  the  continuous 


system  are  shown  in  Figures  4-1  and  4-2.  As  can  be  seen  from 
Figure  4-1,  the  bank  angle  responses  have  very  desirable 
characteristics.  Note  that  there  is  no  visible  difference 


between  the  responses  for  the  three  F.C.s  plotted  in  Figure 
4-1,  thus  showing  that  robust  performance  is  achieved  by  the 
continuous  system.  The  responses  start  out  slowly  and  then 
increases  more  rapidly  as  time  increases.  This  is  desired  in 
a  heavy  transport  since  it  cannot  react  to  rapid  changes  as 
can  a  fighter  aircraft.  The  sideslip  responses  to  a  bank 
angle  command  are  shown  in  Figure  4-2.  It  can  be  seen  that 
the  magnitude  of  the  sideslip  responses  are  quite  small,  as 
desired.  The  characteristics  of  these  responses  are  given  in 
Table  4-1.  Comparison  of  these  responses  with  those  of 
reference  1  shows  that  they  are  identical  for  all  practical 
purposes.  In  addition  all  specifications  are  met  as  can  be 
seen  by  comparing  the  results  shown  in  Table  4-1  with  the 
specifications  given  in  Appendix  D. 

Sideslip  Command  Response  -  Continuous  System 

The  sideslip  command  responses  for  the  continuous 
system  are  shown  in  Figures  4-3  and  4-4.  The  sideslip 
responses  in  Figure  4-3  also  show  the  characteristics  desired 
of  a  large  aircraft.  The  magnitude  of  the  bank  angle 
responses  in  Figure  4-4  are  small  as  desired.  The  figures  of 
merit  for  the  continuous  system  sideslip  responses  are  shown 
in  Table  4-2.  Comparison  of  Table  4-2  with  the  corresponding 
specifications  given  in  Appendix  D  show  that  the  sideslip 
command  response  specifications  are  met.  Comparison  with  the 


Table  4-1 


Bank  Angle  Command  -  Continuous  System 


Response 

F.C. 

Peak 

Value 

Final 

Value 

Rise 

Time  (sec) 

Settling 
Time  (sec) 

1 

30.0 

30.0 

3.53 

6.17 

0  (deg) 

2 

30.0 

30.0 

3.52 

6.17 

3 

30.0 

30.0 

3.51 

6.14 

1 

0.0025 

0.0015 

— 

— 

(deg) 

2 

-0.0019 

0.0015 

— 

— 

3 

0.0361 

0.0178 

— 

— 

responses 

of  reference  1 

indicates  the  systems  are  equivalent 

Bank  Angle  Command  Response  -  Hybrid 

System 

The 

hybrid  system  i 

responses  for 

a  bank  angle  command 

are  shown 

in  Figures  4-5 

and  4-6.  The  bank  angle 

responses 

are  shown 

in  Figure  4-5. 

Figure  4-6 

shows  the  sideslip 

responses  for  the  hybrid  system.  A  comparison  of  the  hybrid 
system  responses  with  the  continuous  system  shows  that 
essentially  equivalent  performance  has  been  achieved  with  the 
hybrid  system.  Table  4-3  gives  the  characteristics  of  the 
bank  angle  responses. 

Sideslip  Command  Response  -  Hybrid  System 

The  responses  for  the  sideslip  command  are  given  in 
Figures  4-7  and  4-8  for  the  hybrid  system.  Figure  4-7  shows 
the  sideslip  respon-^es.  The  bank  angle  responses  are  shown  in 
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System  Sideslip  Respon 


SB 


Response  F.C. 


^  (deg)  2 
3 


0  (deg)  2 
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Table  4-3 

Bank  Angle  Command  -  Hybrid  System 


Response  F.C, 


Peak 

Pinal 

Rise 

Settling 

Value 

Value 

Time  (sec) 

Time  (sec) 

3 


0.0361 


0.0178 


TIME  (SEC) 


Figure  4-8.  The  sideslip  response  figures  of  merit  for  the 
hybrid  system  are  shown  in  Table  4-4.  Comparison  of  both  the 
response  plots  and  figures  of  merit  shows  that  the  hybrid 
system  performance  is  essentially  equivalent  to  the  continuous 
system. 

Table  4-4 

Sideslip  Command  -  Hybrid  System 


Response  F.C. 


fi  (deg) 


Peak 

Final 

Rise 

Settling 

Value 

Value 

Time  (sec) 

Time  (sec) 

4.9 

4.9 

4.95 

8.39 

4.9 

4.9 

4.95 

8.39 

4.9 

4.9 

4.92 

8.36 

-0.319 

-0.319 

— 

— 

-0.337 

-0.337 

— 

— 

-0.179 

-0.179 

— 

0  (degl 


IV-4  KC-135  Longitudinal  Mode  Responses 

A  single  maneuver  is  performed  for  the  longitudinal 
mode  rigid  and  elastic  body  aircraft.  Two  different  inputs 
were  used  in  reference  1  for  the  rigid  body  aircraft  simula¬ 
tions.  The  input  for  F.C.  3  was  a  4  degree  pitch  angle 
command  which  ramps  to  4  degrees  in  1  second  and  remains  at 
4  degrees  thereafter.  For  F.C.'s  1  and  2,  a  similar  input  was 
used  but  the  command  only  ramps  to  1  degree.  The  reason  for 

the  two  different  inputs  is  explained  in  reference  1  (1:131). 
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To  allow  a  direct  comparison  to  be  made  with  the  results  of 
reference  1,  the  simulations  in  this  thesis  investigation  use 
the  same  inputs  as  those  in  reference  1.  For  the  elastic 
body  simulations,  the  input  for  all  three  F.C.'s  is  a  1 
degree  pitch  command  which  rises  to  1  degree  in  1  second  (the 
same  as  reference  1).  Also  note  that  the  rigid  and  elastic 
body  responses  are  shown  on  the  same  plots.  As  reported  in 
reference  1,  there  was  very  little  difference  between  the 
rigid  and  elastic  body  responses  for  the  continuous  system  due 
to  the  small  magnitudes  of  the  first  and  second  body  bending 
modes  (1:97).  As  seen  in  the  responses  on  the  following 
pages,  there  is  very  little  difference  between  the  rigid  and 
elastic  body  responses  for  the  hybrid  system  as  well. 

The  continuous  system  responses  and  figures  of  merit  for 
the  KC-135  longitudinal  mode  are  shown  in  Appendix  F.  A 
comparison  between  these  responses  and  the  those  contained  in 
reference  1  show  that  the  specifications  are  met  and  that  the 
two  continuous  systems  are  equivalent. 

The  responses  to  the  pitch  angle  command  for  the  hybrid 
system  are  shown  in  Figures  4-9  through  4-11.  The  pitch 
angle  response  is  shown  in  Figure  4-9.  The  pitch  angle 
response  is  similar  to  the  tracking  responses  of  the  lateral 
mode  in  that  there  is  a  slow  initial  rise  followed  by  a  more 
rapid  increase  in  the  output.  The  perturbation  altitude  and 
perturbation  velocity  responses  are  shown  in  Figures  4-10  and 
4-11,  respectively.  Note  that,  as  in  the  continuous  system, 
the  magnitudes  of  the  altitude  and  velocity  responses  are 


TIME  (SEC) 


TIME  (SEC) 


Fig.  4-10.  Hybrid  System  Altitude  Responses 
(compare  to  Fig.  F-2) 


TIME  (SEC) 


Fig.  4-11.  Hybrid  System  Velocity  Responses 
(compare  to  Fig.  F-3) 

small.  Table  4-5  shows  the  time  response  characteristics  for 
the  rigid  body  pitch  angle  command  responses  while  Table  4-6 
shows  the  corresponding  results  for  the  elastic  body. 

IV-5  F-16  Simulation  Responses 

Stable  responses  are  obtained  for  F.C.'s  1  and  2  only 
for  the  F-16  as  discussed  in  Chapter  3.  Therefore  the  plots 
for  the  F-16  contain  responses  for  F.C.'s  1  and  2  only  and 
for  each  of  these  two  F.C.'s  there  are  6  CSC  modes  plotted. 

The  input  for  the  first  maneuver  is  a  pitch  rate  pulse  command 
of  10  degrees  per  second  and  the  input  for  the  second  maneuver 
is  a  roll  rate  pulse  command  of  50  degrees  per  second.  In 
both  cases,  the  input  is  removed  at  5  seconds  into  the 


simulation . 


Table  4-5 


Pitch  Angle  Command  -  Hybrid  System  (Rigid  Body) 


Response 

F.C. 

Peak 

Value 

Final 

Value 

Rise 

Time  (sec) 

Settling 
Time  (sec) 

1 

1.0 

1.0 

4.35 

8.28 

0  (deg) 

2 

1.0 

1.0 

4.35 

8.27 

3 

3.9 

3.9 

4.38 

8.32 

1 

-0.040 

0,026 

— 

— 

h  (ft) 

2 

0.222 

0,222 

— 

— 

3 

-0.384 

0.093 

— 

— 

1 

-0,088 

-0.007 

— 

— 

u  (ft/sec) 

2 

-0.137 

-0.010 

— 

— 

3 

0.069 

0.015 

— 

— 

The  continuous  domain  designs  in  reference  2  used  a 
response  model,  as  shown  in  Appendix  D  (Tables  D-13  and  D-16), 
to  develop  disturbance  rejection  bounds  needed  for  the  QFT 
design  process.  These  models  were  expressed  in  terms  of  angle 
rate  variables,  similar  to  the  tracking  response  models.  To 
evaluate  the  disturbance  rejection  performance  of  the  designs, 
a  total  angle  was  used  rather  than  angle  rate.  For  example, 
if  pitch  rate  is  commanded,  then  pitch  rate  is  the  tracking 
response  and  roll  rate  is  the  disturbance  rejection  response. 
In  reference  2,  the  performance  of  the  disturbance  rejection 
response  was  evaluated  by  comparing  the  roll  angle  to  the 


Table  4-6 


Pitch  Angle  Command  -  Hybrid  System  (Elastic  Body) 


Response 

F.C. 

Peak 

Value 

Final 

Value 

Rise 

Time  (sec) 

Settling 
Time  (sec) 

1 

1.0 

1.0 

4.36 

8.28 

e  (deg) 

2 

1.0 

1.0 

4.36 

8.27 

3 

1.0 

1.0 

4.38 

8.34 

1 

0.026 

— 

— 

h  (ft) 

2 

0.222 

0.222 

— 

— 

3 

— 

— 

1 

— 

— 

u  (ft/sec) 

2 

-0.137 

— 

— 

3 

0.017 

0.004 

— 

— 

roll  angle  specification.  Because  of  the  state  limitation  of 
the  computer  program  used  in  reference  2,  the  roll  angle  was 
estimated  by  graphically  calculating  the  area  under  the  roll 
rate  response  curves  resulting  from  cross-coupling  (2:60). 

As  shown  in  the  above  referenced  tables,  it  is  desired  that 
roll  angle  be  less  than  five  degrees  when  pitch  rate  is 
commanded.  Similarly,  when  roll  rate  is  commanded,  it  is 
desired  that  the  pitch  angle  be  less  than  three  degrees. 

The  computer  program  used  for  this  thesis  investigation 
has  no  state  limitations,  therefore  angle  (versus  angle  rate) 
disturbance  rejection  responses  are  obtained  directly  by 


integrating  the  angle  rate  responses.  Thus  only  the  angle 
response  curves  are  shown. 

As  in  the  case  of  the  KC-135,  continuous  domain  simula¬ 
tions  are  performed  for  the  F-16.  These  responses  and  their 
associated  figures  of  merit  are  shown  in  Appendix  F.  A 
comparison  of  these  continuous  system  responses  with  those  of 
reference  2  shows  that  the  systems  performed  similarly 
although  several  of  the  performance  discrepancies  reported  in 
reference  2  (referenced  earlier  in  this  chapter)  are  not 
experienced.  Again,  this  is  most  probably  due  to  the  use  of 
the  full  order  compensators  in  this  investigation. 


Pitch  Rate  Command  Response  -  Hybrid  System 


The  pitch  rate  command  responses  are  shown  in  Figures 
4-12  and  4-13.  Figure  4-12  shows  the  pitch  rate  response. 

The  figure  shows  a  relatively  quick  response,  as  desired  for 
a  fighter  aircraft.  The  roll  angle  resulting  from  the  pitch 
rate  command  is  shown  in  Figure  4-13.  From  the  figure  it  can 
be  seen  that  the  maximum  roll  angle  is  approximately  3.2 
degrees,  well  below  the  5  degree  specification.  The  figures 
of  merit  for  the  pitch  rate  responses  are  shown  in  Table  4-6. 


Roll  Rate  Command  Response  -  Hybrid  System 


The  responses  to  the  roll  rate  command  are  shown  in 
Figures  4-14  and  4-15.  The  roll  rate  response  is  shown  in 
Figure  4-14  where  the  response  is  similar  to  that  for  pitch 
rate.  The  pitch  angle  response  is  shown  in  Figure  4-15. 

The  maximum  roll  angle  is  approximately  -1.9  degrees  which 


PITCH 


Table  4-8 

Roll  Rate  Response  Characteristics  -  Hybrid  System 


CSC  Mode 


Peak  Value 
( deg/sec ) 


Rise 

Time  (sec) 


Settling 
Time  (sec) 


meets  the  3  degree  specification.  The  roll  rate  figures  of 
merit  are  shown  in  Table  4-7.  Examination  of  the  response 
data  shows  that  the  specifications  are  met  for  both  F.C.'s  1 
and  2,  demonstrating  partial  robustness. 


IV-6  Summary 

The  responses  for  the  KC-135  illustrate  that  all 
performance  objectives  are  met  for  the  hybrid  system.  The 
response  data  also  shows  that  there  is  very  little  variation 
over  all  F.C.'s  thereby  demonstrating  the  robustness  of  the 
digital  controllers.  The  performance  of  the  KC-135  hybrid 
"infinite"  wordlength  system  matches  that  of  the  continuous 
system  well.  The  close  agreement  between  the  continuous  and 
hybrid  system  responses  suggests  that  a  slower  sampling  rate 
might  produce  acceptable  results  as  well.  The  F-16  responses 
show  that  for  F.C.'s  1  and  2  the  desired  responses  are 
achieved.  The  variation  over  the  two  F.C.'s  is  very  small, 
thereby  achieving  limited  robustness. 
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V.  Implementation  Phase 


V-1  Introduction 

This  chapter  discusses  the  digital  implementation  of  the 
KC-135  lateral  and  longitudinal  mode  control  systems.  In  the 
implementation  phase  both  the  sampling  and  finite  wordlength 
effects  encountered  in  a  digital  control  system  are  consi¬ 
dered.  Some  of  the  effects  of  finite  wordlength  are  briefly 
discussed.  The  primary  means  by  which  these  effects  are 
considered  is  the  use  of  a  computer  program  which  emulates  the 
finite  wordlength  arithmetic  of  microprocessor.  The  develop¬ 
ment  of  and  the  theory  behind  this  computer  program  are 
presented.  The  last  portion  of  this  chapter  is  devoted  to  the 
results  obtained  from  simulations  performed  using  the  finite 
wordlength  simulation  program.  Included  in  this  discussion  is 
the  choice  of  a  microprocessor  wordlength  which  provides  the 
required  control  system  performance. 

V-2  Finite  Wordlength 

In  a  digital  control  systei  th  control  compensator  is 
implemented  on  a  microprocessor  (digital  computer).  The 
calculation  of  the  discrete  compensation  signal  which  drives 
the  analog  plant  is  performed  by  the  digital  computer  and 
because  of  the  finite  wordlength  of  the  digital  computer, 
errors  are  generated  during  this  calculation  (9:319).  In  this 
sense  "error*  means  that  the  discrete  compensation  signal  is 
different  than  the  corresponding  analog  compensation  signal 
(which  is  considered  to  be  exact).  Most  of  these  errors  can 


generally  be  attributed  to  three  different  sources:  quantiza¬ 
tion  of  the  analog  input  to  the  digital  controller  (9:320), 
quantization  of  the  compensator  coefficients  (9:336)  and  the 
quantization  of  the  results  of  arithmetic  operations  (9:330). 
In  all  three  cases  the  errors  occur  because  a  "number"  with 
high  accuracy  (e.g.  "infinitely"  accurate  analog  signal  or  a 
number  with  many  decimal  digits)  must  be  represented  with  a 
relatively  small  number  (finite)  of  binary  bits.  Analytical 
methods  exist  which  allow  the  generation  of  models  for  each  of 
the  error  sources  mentioned  above.  Usually  these  error  models 
are  described  by  nonlinear  functions  or  stochastic  processes. 
These  models  can  be  used  to  determine  the  effect  of  the  error 
sources.  Rather  than  use  error  models,  this  thesis  investiga¬ 
tion  uses  computer  simulations  to  examine  control  system 
performance  where  the  effects  of  microprocessor  errors  are 
included  in  the  simulation. 

The  implementation  structure  of  the  digital  controller 
can  affect  the  performance  of  the  digital  control  system.  A 
order  digital  controller  is  shown  in  Equation  (5-1)  v'here 
some  of  the  numerator  coefficients  may  be  zero. 

Y(z)  a„  +  a,z”^  +  ...  +  a.z"*^ 

G(z)  = -  - ^ (5-1) 

R(z)  1+b^z  +...+b^z 

Equation  (5-1)  may  be  rewritten  as  in  Equation  (5-2)  where 
this  form  is  generally  called  the  k  order  difference 
(recurrence)  equation.  Implementation  of  a  direct  structure 
would  involve  programing  the  microprocessor  such  that  the 


coefficients  in  Equation  (5-2),  and  ,  appear  as  multipli¬ 
ers  in  the  computer  code  where  i  =  l,2,3,...k.  Different 
direct  structures  exist  but  the  simplest  form  is  coding  the 
microprocessor  exactly  as  is  seen  in  Equation  (5-2)  (called 
the  direct  form  1  or  ID) .  The  effects  of  the  last  two  error 

Y(z)  =  aQR(z)  +  a^z"^R(z)  +  ...  +  a|^z~'^R(z) 

-  b^2~^Y(z)  -  b2z“^Y(z)  -  ...  -  b|^z“'^Y(z)  (5-2) 

sources  mentioned  previously  can  be  decreased  by  implementing 
the  order  compensator  as  second-order  sections  connected 

either  in  a  cascade  or  parallel  structure.  Many  different 
cascade  and  parallel  structures  exist  with  each  having  various 
advantages  and  disadvantage  (8:360,471).  in  f-is  study  the 
digital  controllers  are  coded  using  the  ID  structure  discussed 
above . 


V-3  Microprocessor  Emulation 

The  implementation  of  the  difference  equation  on  a 
microprocessor,  and  also  in  digital  controller  simulation, 
involves  the  relatively  simple  computer  operations  of  addi¬ 
tion,  multiplication  and  memory  storage  and  retrieval.  To 
show  more  explicitly  what  must  be  stored  and  retrieved  in  the 
computer  memory,  Equation  (5-2)  is  modified  to 

y(n)  =  aQr(n)  +  aj^r(n-l)  +  ...  +  a,^r(n-k) 

-b^y(n-l)  -  b2y(n-2)  -  ...  -  b|^y(n-k)  (5-3) 

where  y(n)  is  the  difference  equation  output  at  the  current 
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sampling  period  and  r(n)  is  the  difference  equation  input  at 

the  current  sampling  period.  Note  that  the  products  of  the 

-1  -2 

delay  operators  (z  ,  z  ,  ect.)  and  the  input  and  output  of 
Equation  (5-2)  are  now  reolaced  bv  y(n-i)  and  r(n-i)  in 
Equation  (5-3)  which  represent  the  difference  equation  output 
and  input  at  the  i^*^  preceding  sample  period.  Thus  the 
previous  inputs  and  outputs  of  the  difference  equation  must 
be  stored  and  subsequently  retrieved  from  memory  in  order 
to  compute  the  difference  equation  at  the  n*"*^  (present) 
sample  period. 

A  microprocessor  represents  numbers  inherently  as 
integers.  Despite  this  fact  the  programmer  of  a  microproces¬ 
sor  has  the  freedom  to  assign  an  implied  radix  point  anywhere 
in  the  binary  word.  This  leads  to  fixed-point  arithmetic  as 
a  natural  means  of  computing  the  difference  equation. 
Fixed-point  arithmetic  is  generally  different  from  that  of 
mainframe  computers  that  are  used  to  perform  dynamic  system 
simulations.  The  mainframes  performing  these  simulations  are 
usually  configured  to  use  floating-point  arithmetic.  Thus  if 
a  microprocessor  based  digital  controller  is  to  be  accurately 
emulated  in  a  simulation  program  on  a  mainframe,  the  computa¬ 
tions  must  be  done  in  such  a  way  as  to  generate  the  sane 
results  as  generated  by  the  microprocessor. 

A  microprocessor  simulation  tool  was  developed  in 
reference  23  which  applied  numerical  analysis  techniques  to 
software  to  obtain  a  desired  accuracy  with  the  constraints  of 
a  given  fixed  wordlength  and  number  representation.  Although 


this  tool  (a  computer  program)  can  perform  much  of  the 
analysis  needed  for  this  study,  it  was  designed  primarily  as 
a  stand-alone  tool  and  could  not  be  integrated  with  other 
simulation  software  without  considerable  time  and  effort. 

Thus,  a  new  microprocessor  simulation  tool  is  developed  to  aid 
in  this  study.  An  algorithm  is  suggested  in  reference  16 
to  perform  analysis  similar  to  that  provided  in  reference  23. 
This  algorithm  is  used  to  write  the  microprocessor  emulation 
program  (actually  a  FORTRAN  subroutine)  utilized  in  this 
thesis  investigation.  The  following  discussion  of  the  theory 
of  the  algorithm  paraphrases  that  of  reference  16. 

The  emulation  of  fixed-point  microprocessor  calculations 
can  be  achieved  by  performing  all  calculations  with  floating¬ 
point  arithmetic  and  then  shifting  the  numbers  binarily  and 
taking  the  integer  part.  This  is  essentially  the  same  as 
performing  truncated  calculations  using  integer  arithmetic 
with  an  implied  decimal  point,  as  is  the  normal  operating 
mode  of  a  microprocessor.  The  radix  point  can  be  placed  in 
any  location  desired.  Rounding,  rather  than  truncation,  can 
be  emulated  by  adding  +.5  prior  to  taking  the  integer  part. 

The  above  discussion  can  be  expressed  in  equation  form  as 

N(x)  =  IP{[x]2P  +  [r]sgn(x)}  (5-4) 

where  x  is  a  decimal  number,  p  is  the  number  of  binary  places 
desired  to  the  right  of  the  radix  point,  r  is  0  for  truncation 
and  0.5  for  rounding,  IP{  }  is  the  integer  part  function  and 
sqn(  )  is  the  sign  function  (16:247).  Note  that  any  portion 


of  the  decimal  number  smaller  than  the  least  significant  bit 
(LSB)  of  the  microprocessor  word  is  truncated  or  rounded  off 
(depending  on  the  value  of  r). 

By  manipulating  all  numbers  involved  in  the  calculation 
of  the  difference  equation  with  Equation  (5-4),  all  computa¬ 
tions  are  performed  with  integers  that  represent  fixed-point 
numbers.  Thus  the  arithmetic  is  the  same  as  that  performed  in 
the  microprocessor  which  is  being  emulated.  Specifically, 
Equation  (5-4)  must  be  applied  to  the  difference  equation 
input  and  the  coefficients.  If  overflows  are  managed  correct¬ 
ly,  then  the  result  of  adding  integers  is  the  same  as  would  be 
obtained  from  the  microprocessor.  If  the  computer  wordlength 
is  L+1  (one  bit  is  reserved  for  the  sign  bit),  the  magnitude 
of  the  sum  of  anv  two  integers  cannot  exceed  2^-1.  If  a 
positive  overflow  occurs,  then  the  magnitude  of  the  number  is 
set  to  2^-1.  If  a  negative  overflow  occurs,  then  the  magni¬ 
tude  must  be  set  to  2^.  Note  that  it  is  assumed  that  the 
microprocessor  manages  overflows  as  described  above.  In  some 
microprocessors  this  saturation  capability  is  ayailable  as  a 
inherent  function  (such  as  the  Texas  Instruments  TMS32010)  and 
on  others  it  must  be  implemented  in  the  microprocessor  code. 

Additional  massaging  must  be  performed  on  the  results 
of  multiplications.  When  two  fixed-point  numbers  are  multi¬ 
plied,  a  fixed-point  number  results  which  has  twice  the 
number  of  binary  digits  to  the  right  of  the  radix.  To 
properly  massage  this  number  using  truncation,  the  programmer 
of  the  microprocessor  shifts  the  results  binarily  to  the 


right  by  the  number  of  binary  digits  assumed  to  the  right  of 
the  radix.  The  same  manipulation  is  performed  when  rounding  is 
desired  except  that  prior  to  performing  the  shifting  opera¬ 
tion,  the  programmer  adds  1  to  the  resulting  product.  The 
process  described  above  for  manipulating  products  can  be 
implemented  on  a  mainframe  computer  using  the  following 
function : 

T(x)  =  IP{[x]2'P  +  Ir]sgn(x)}  (5-5) 

where  x  is  now  the  result  of  performing  multiplication  using 
floating-point  arithmetic  on  two  numbers  which  have  already 
been  massage  by  Equation  (5-4)  (16:247).  Note  that  IP{  }, 
sgn(  ),  r  and  p  are  as  defined  previously.  Floating-point 
arithmetic  must  be  used  to  perform  all  computations  inside  the 
braces  of  Equation  (5-5).  Multiplications  can  also  produce 
overflows.  If  an  overflow  does  occur,  the  number  resulting 
from  the  application  of  Equation  (5-5)  must  be  manipulated  in 
the  same  fashion  as  discussed  previously  for  the  addition 
operation.  Examples  illustrating  the  application  of  Equations 
(5-4)  and  (5-5)  to  specific  numbers  are  given  in  reference 
16, 

A  FORTRAN  subroutine  is  written  (24)  which  implements 
the  algorithm  discussed  above.  MATRIX^^  does  not  call  this 
subroutine  directly,  if  MATRIX^^  is  "told"  that  a  FORTRAN 
subroutine  is  to  be  used  in  a  simulation,  it  calls  a  generic 
utility  subroutine.  This  utility  subroutine  then  calls  all 
user  written  subroutines.  The  source  code  for  the  utility 


subroutine  is  provided  in  skeleton  form  with  the  MATRIX^ 
executable  (compiled)  code.  After  compilinq  the  two  subrou¬ 
tines  they  must  be  linked  to  an  image  of  the  MATRIX^^  execu¬ 
table  code.  A  listing  of  both  subroutines  and  a  brief 
discussion  of  the  microprocessor  emulation  subroutine  is 
contained  in  Appendix  G.  Pseudocode  is  included  in  the 
subroutine  as  well. 

The  capabilities  of  the  subroutine  can  be  summarized  as 
follows:  Any  microprocessor  wordlength  can  be  specified  up 

to  approximately  60  bits  (Quad  precision  functions  must  be 
used  for  higher  wordlenqths ) .  Either  truncation  or  rounding 
can  be  chosen  for  quantizations.  The  subroutine  accepts  a 
floating-point  input  which  is  considered  to  be  an  "infinite" 
precision  number  equivalent  to  the  analog  error  signal  of  the 
control  system  at  the  present  sample  period  (E*(s)  in  Figure 
2-3),  The  subroutine  then  generates  the  present  output  of  the 
difference  equation  as  a  floating-point  number  (compensation 
signal).  It  is  assumed  that  the  wordlength  of  the 
analoq-to-digital  (A/D)  converter  is  the  same  as  the  micro¬ 
processor  wordlength.  This  is  done  for  simplicity  but 
the  subroutine  can  be  modified  rather  easily  to  remove  this 
limitation  in  future  studies.  Thus  MATRIX^^  performs  the 
sampling  operation,  the  subroutine  performs  t'-.e  "A/D  conver¬ 
sion"  (quantizes  the  input),  and  MATRIX^^  models  the  ZOH  at  the 
digital  controller  output  using  Equation  (2-1).  Therefore 
the  diagrams  of  the  system  models  contained  in  Chapter  II  are 
still  applicable  except  that  the  microprocessor  emulation 


subroutine  is  perforninq  the  comoensator  calculations  usinq 
finite  wordlength  fixed-point  arithmetic. 

V-4  Implementation  Phase  Results 

This  section  contains  the  responses  for  both  the  lateral 
and  longitudinal  mode  of  the  KC-135  with  the  effects  of  finite 
wordlength  included  in  the  simulations.  The  results  are 
obtained  usinq  the  microprocessor  emulation  subroutine 
discussed  in  the  previous  section.  The  inputs  to  both  modes 
of  the  aircraft  are  the  same  as  that  described  in  Sections 
IV-3  and  IV-4  (i.e.  the  same  maneuvers  are  performed). 

V-4-1  KC-135  Lateral  Mode  Responses  The  lateral 

mode  responses  are  obtained  using  a  microprocessor  wordlength 
of  32  bits  with  17  bits  to  the  right  of  the  radix  point 
(17  BRP).  Rounding  rather  than  truncation  is  used  in  the 
simulations  generating  these  plots.  All  performance  specifi¬ 
cations  are  met  usinq  this  configuration  as  is  shown  later  in 
the  section.  Thirty-two  bits  is  a  relatively  standard 
wordlength.  The  number  of  BRP  (17  bits)  is  chosen  because 
this  number  provides  similar  performance  for  the  longitudinal 
mode  simulations  and  thus  a  common  configuration  is  achieved 
for  the  entire  aircraft.  Simulations  are  run  using  several 
different  combinations  of  rounding/truncation  and  numbers  of 
BRRs  with  the  overall  wordlength  kept  at  32  bits.  The 
"qualitative  results"  of  these  simulations  as  a  function  of 
BRR  and  truncation/roundinq  is  shown  in  Table  5-1. 


Wordlenqth 


Bits  Riqht 
of  Radix 
Point  (BRR) 


24 

unacceptable 

unacceptable 

20 

unacceptable 

unacceptable 

18 

unacceptable 

unacceptable 

17 

acceptable 

acceptable 

10 

acceptable 

acceptable 

6 

acceptable 

acceptable 

4 

acceptable 

unacceptable 

3 

unacceptable 

unacceptable 

In  Table  5-1,  "acceptable"  system  performance  means  that 
all  specifications  are  met  while  "unacceptable"  means  that 
either  some  or  all  of  the  specifications  are  not  met  or  that 
the  system  is  unstable.  From  the  table  it  can  be  seen  that 
the  selection  of  truncation  or  rounding  for  the  lateral  mode 
makes  little  difference  in  the  overall  performance.  Also 
from  the  table  it  is  seen  that  the  fewest  number  of  BRR  for 
acceptable  performance  is  4.  This  implies  that  at  least  4 
bits  ate  needed  to  maintain  the  necessary  precision.  Examina¬ 
tion  of  the  table  also  reveals  that  the  fewest  number  of  bits 
to  the  left  of  the  radix  point  that  results  in  acceptable 
performance  is  15  (32  -  17).  This  implies  that  at  least  15 
bits  are  needed  to  provide  the  necessary  magnitude  (maximum 
integer).  Thus  it  appears  that  a  19  bit  microprocessor 
probably  provides  acceptable  performance.  Examination  of 
Equations  (3-11)  and  (3-12)  show  that  the  compensator  gains 


for  the  KC-135  lateral  mofie  are  relatively  hiah.  If  the 
majority  of  this  qain  is  inplenented  outside  the  digital 
controller  (e.q.  analog  amplifier)  it  is  quite  probable  that 
a  16  bit  (common  wordlength)  microprocessor  achieves  accepta¬ 
ble  performance  since  the  required  number  of  "magnitude 
bits"  is  reduced.  The  above  analysis  is  accomplished  using 
the  rule  of  thumb  that  approximately  3.5  bits  are  required  to 
represent  each  decimal  digit. 

Bank  Angle  Command  Response 

The  figures  of  merit  for  the  bank  angle  command  are 
given  in  Table  5-2.  The  bank  angle  command  responses  are 
shown  in  Figures  5-1  and  5-2.  The  bank  angle  response  is 
shown  in  Figure  5-1  and  the  sideslip  response  is  shown  in 
Figure  5-2.  Examination  of  Table  5-2  shows  that  all  perfor¬ 
mance  specifications  are  met.  Also,  the  responses  shown  in 
Figure  5-1  have  the  shape  desired  for  the  tracking  response 
as  discussed  in  ChaPter  IV.  Comparison  with  the  continuous 
domain  response  data  from  Chapter  IV,  repeated  here  as  Table 
5-3,  shows  that  the  results  are  very  similar. 

Sideslip  Command  Response 

The  figures  of  merit  for  the  sideslip  command  are  given 
in  Table  5-4.  As  with  the  bank  angle  command  responses,  all 
specifications  are  met  for  the  sideslip  command  responses. 

The  sideslip  command  responses  are  presented  in  Figures  5-3 
and  5-4.  The  sideslip  response  to  a  sideslip  command  is 
shown  in  Figure  5-3.  Figure  5-4  shows  the  bank  angle  response 


Table  5-2 


Bank  Angle  Command  -  Finite  Wordlength 


Response 


96  (deg) 


F.C, 


Peak 

Value 


Final 

Value 


Rise  Settling 

Time  (sec)  Time  (sec) 


30.0 

30.0 

30.0 


31 

3(1 

3r 


P  (deg: 


-0.0030 

-0.0019 

0.0361 


0.0015 

0,0015 

0,0178 


Table  5-3 

Bank  Angle  Command  -  Continuous  System 


Response 

mm 

1 - 

Peak 

Value 

Final 

Value 

Rise 

Time  (sec) 

Settling 
Time  (sec) 

1 

30.0 

30.0 

3.53 

6.17 

0  (deg) 

2 

30.0 

30.0 

3.52 

6.17 

3 

30.0 

30.0 

3.51 

6.14 

Hj 

0.0025 

0.0015 

— 

— 

P  (deg) 

-0.0019 

0.0015 

— 

— 

0.0361 

0.0178 

-- 

— 

Response 


B  ( deq  ) 


(8  (deq) 


Table  5-4 

Sideslip  Connand  -  Finite  Wordlenqth 


Peak 

Final 

Rise 

Settlinq 

Value 

Value 

Time  (sec) 

Time  (sec) 

1.4 

4 . 9 

4.9 

4.95 

8 .38 

4,9 

4.9 

4.95 

8 , 39 

4.9 

4.9 

4.95 

8 ,35 

1.  ,  * 

-0.319 

-0.319 

-- 

— 

-0.337 

-0.337 

-- 

-- 

-0.179 

-0.179 

— 

— 

Table  5-5 

Sideslip  Command  -  Continuous  System 


Response  |  F.C, 


p  (deq) 


Peak 

Value 


Final  Rise  Settlinq 

Value  Time  (sec)  Time  (sec) 


0  (deq) 


-0.319 

-0,337 

-0.179 


-0.319 

-0,337 

-0.179 


maintain  the  needed  precision.  This  can  be  partly  attributed 
to  the  fact  that  the  lonqituriinal  mode  plant  equations  contain 
unstable  (right-half  plane)  poles  and  demand  more  precision 
to  control  them  than  the  lateral  mode.  The  table  also  shows 
more  variation  in  the  results  between  the  use  of  truncation  or 
rounding.  Rounding  provides  a  laraer  range  of  acceptable 
performance.  Performing  a  similar  analysis  as  is  done  for 
the  lateral  mode  reveals  that  a  minimum  of  26  bits  can  produce 
the  same  results  shown  in  the  responses  on  the  following 
pages.  It  is  possible,  but  less  likely,  that  a  16  bit 
microprocessor  can  achieve  acceptable  results  for  the  longitu¬ 
dinal  mode.  Implementing  a  majority  of  the  compensator 
gain  outside  the  digital  controller  might  allow  a  16  bit 
microprocessor  especially  considering  that  the  magnitude  of 
the  gain  for  one  of  the  longitudinal  compensators  is  higher 
than  the  maximum  lateral  compensator  gain. 

The  figures  of  merit  for  the  rigid  body  longitudinal 

mode  responses  are  shown  in  Table  5-7.  The  elastic  body 

figures  of  merit  are  shown  in  Table  5-8.  The  pitch  command 

responses  for  the  longitudinal  mode  are  shown  in  Figures  5-4 

through  5-6.  The  pitch  angle  response  is  shown  in  Figure  5-4. 

Figure  5-5  shows  the  perturbation  altitude  response  and  Figure 

5-6  shows  the  perturbation  velocity  response.  Note  that  as 
before  both  the  rigid  and  elastic  body  responses  for  all  three 

F.C.'s  are  shown  on  the  same  figure.  The  response  character¬ 
istics  given  in  both  tables  indicates  the  all  performance 
specifications  are  met.  Examination  of  the  tracking  response 


Table  5-8 

Pitch  Angle  Command  (Elastic  Body) 


Response  | F . C , 


u  ( ft/sec)  I  2 
3 


Peak 

value 

Final 

Value 

Rise 

Time  (sec) 

Settling 
Time  (sec) 

1.0 

1.0 

4.36 

8.27 

1.0 

1.0 

4.36 

8.28 

3.9 

3.9 

4.33 

7.78 

-0.040 

0.025 

— 

— 

0.221 

0.221 

— 

— 

-0.101 

0.024 

— 

— 

-0.088 

-0.043 

— 

— 

-0.138 

-0.064 

— 

— 

0.019 

0.014 

— 

— 

TIME  (SEC) 


Fig.  5-7.  Velocity  Responses 
(compare  to  Fig.  F-3) 


V-5  Summary 

The  response  data  presented  in  this  chapter  shows  that 
the  sampled-data  finite  wordlength  KC-135  control  system 
achieves  all  performance  specifications  and  performs  similarly 
to  the  continuous  system.  As  with  the  continuous  system,  the 
elastic  body  responses  show  that  the  addition  of  the  first 
and  second  body  bending  modes  makes  very  little  difference  as 
compared  to  the  rigid  body  aircraft.  This  performance  has 
been  achieved  using  a  fairly  standard  microprocessor  word- 
length  of  32  bits.  An  empirical  analysis  indicates  that  a 
smaller  wordlength  (possibly  16  bits)  will  achieve 
"acceptable"  performance. 


VI .  Conclusions  and  Recomnendations 


VI -1  Discussion 

The  primary  purpose  of  this  thesis  investigation  is  to 
extend  two  quantitative  feedback  theory  analog  multivariable 
control  system  designs  to  the  discrete  domain.  In  addition  to 
extending  the  analog  designs  to  digital  designs,  another 
objective  is  to  determine  what  design  guidelines,  if  any,  can 
be  established  which  might  predict  the  "success"  or  "failure" 
of  extending  an  existing  analog  QFT  design  to  the  discrete 
domain . 

The  preceding  chapters  indicated  that  in  one  case  (the 
KC-135)  the  analog  system  can  be  extended  to  a  digital  system 
which  meets  all  specifications  and  achieves  essentially  the 
same  performance  as  the  continuous  system.  For  the  c  her  case 
(AFTI/F-16),  specifications  are  met  for  only  one  out  of  four 
flight  conditions.  Simulations  show  that  the  system  is 
unstable  at  the  other  three  flight  conditions  for  reasons  to 
be  discussed  later.  Performing  "cut  and  try"  (i.e.  trial  and 
error)  modifications  to  one  of  the  F-16  digital  controllers 
improved  the  performance  so  that  a  stable  system  is  achieved 
and  specifications  are  met  for  two  out  of  four  flight  condi¬ 
tions  . 

In  both  cases  the  Tustin  transformation,  in  conjunction 
with  a  PCT  system  approximation,  is  used  to  extend  the  analog 
systems  to  digital  systems.  Depending  upon  the  locations  of 
the  compensator  s-plane  poles  and  zeros,  the  frequency  domain 


characteristics  of  some  real  poles  and  zeros  can  be  lost  or 
severely  warped  by  the  Tustin  transformation.  For  the 
purposes  of  extending  QFT  analog  designs  to  the  digital 
domain,  it  is  surmised  that  the  criterion  for  this  loss  of 
effect  occurs  for  poles  and  zeros  which  are  located  at 
S  ^  -2/T.  In  addition,  when  this  loss  of  effect  occurs  for 
compensator  zeros ,  a  poor  performing  or  unstable  system  will 
likely  result  because  of  a  net  increase  in  phase  lag.  It 
should  be  stressed  that  the  criterion  for  a  successful 
extension  of  an  analog  QFT  design  to  a  digital  one  is  based 
upon  experience  from  extending  only  two  QFT  analog  designs. 
Thus  it  is  difficult  to  determine  at  this  point  whether  the 
above  criterion  will  affect  the  extension  of  analog  QFT 
designs  in  most  cases  or  in  very  few  cases.  It  is  clear 
however  that  the  criterion  will  most  likely  come  into  play 
when  there  is  a  "relatively  high  degree"  of  uncertainty  in 
the  controlled  plant  (e.g.  healthy  KC-135  versus  a  reconfigur- 
able  F-16  with  a  larger  flight  envelope). 

VI -2  Conclusions 

Based  upon  the  results  and  discussions  of  the  previous 
chapters  the  following  conclusions  can  be  stated: 

1.  Use  of  the  Tustin  transformation  and  PCT  approxima¬ 
tion  to  extend  existing  QFT  analog  designs  to  the  discrete 
domain  provide  acceptable  system  performance  when  all  real 
compensator  zeros  are  located  at  S  ^  -2/T.  This  criterion  may 
apply  to  real  compensator  poles  as  well  depending  upon  the 
location  of  the  Nichols  chart  bounds  (other  than  the  UHFB). 


that  if  the  criterion  is  met,  the  Tustin/PCT  apor ox inat ion 
effectively  accounts  for  the  problems  associated  with 
sampling,  finite  wordlength  and  ZOH  delay  inherent  in  a 
sampled-data  control  system. 

2.  If  the  above  criterion  is  not  met  then  the 
Tustin/PCT  approximations  provide  "poor  performance"  (i.e. 
predefined  performance  specifications  are  not  met)  and 
possibly  an  unstable  system.  Improving  the  performance  then 
requires  the  application  of  "cut  and  try"  methods.  Ih  this 
situation,  the  system  performance  guarantees  inherent  in  the 
QFT  design  technique  are  lost  and  acceptable  performance  may 
not  be  achieved  for  all  plants  in  the  parameter  space. 

3.  The  standard  z  transformation  is  not  a  useful  tool 
in  extending  QFT  analog  control  system  designs  to  the  discrete 
domain . 

4.  Using  the  extension  method  described  in  this  study, 
realizable  sample  times  (40  to  60  HZ)  and  microprocessor 
wordlengths  (16  to  32  bits)  provide  acceptable  performance  for 
the  resulting  QFT  digital  control  systems. 

VI -3  Recommendations 

1.  The  primary  recommendation  is  to  implement  the  QFT 
design  technique  directly  in  the  discrete  domain.  It  is  a 
clear  result  of  this  thesis  investigation  that  acceptable 
discrete  system  performance  cannot  be  achieved  systematically 
in  all  cases  by  extending  analog  QFT  designs  to  the  discrete 
domain  using  the  Tustin/PCT  approach.  A  significant  and 


desirable  characteristic  of  the  O^T  technique  is  that  effec¬ 
tive  multivariate  control  systems  can  be  designed  in  a 
systematic  and  deterministic  manner  for  plants  that  have 
uncertain  parameters.  Therefore  to  maintain  these  desirable 
characteristics  the  technique  should  be  applied  directly  in 
the  discrete  domain.  Either  the  z  or  w'  domain  could  be  used 
but  the  w'  domain  is  recommended  since  standard  frequency 
domain  techniques  can  be  used.  Discrete  domain  design 
requires  that  the  plant  equations  be  transformed  to  the 
discrete  domain  which  implies  that  a  sample  rate  must  be 
chosen.  To  avoid  "false  starts"  it  is  suggested  that  the 
maximum  tolerable  sample  rate  be  used  initially.  A  possible 
method  to  get  around  this  restriction  would  be  to 
"parameterize"  the  plant  equations  as  a  function  of  the  sample 
rate.  This  would  provide  a  much  more  general  technique  but  it 
may  not  prove  to  be  very  practical  (i.e.  the  parameterized 
equations  could  prove  difficult  to  work  with  due  their 
complexity) . 

2.  A  secondary  recommendation  is  to  expand  the  scope 
of  the  work  done  in  this  thesis  effort.  Frequency  prewarping 
could  be  applied  to  the  compensators  prior  to  applying  the 
Tustin  transformation.  Instead  of  using  the  Tustin  transfor¬ 
mation  the  matched  z  transformation  (8:339)  could  be  tried. 
This  transformation  is  relatively  "simple"  to  apply  and  it 
will  not  generate  real  poles  or  zeros  on  the  negative  z-plane 
axis  as  does  the  Tustin. 
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Appendix  A;  Single  Input-Single  Output  Theory 
(This  appendix  was  taken  from  Reference  1  with  minor  changes.) 

Introduction 

Appendices  A  and  B  present  an  overview  of  the  Quanti¬ 
tative  Feedback  Technique  used  in  the  design  of  the  multiple 
output  flight  control  systems  in  this  thesis.  Examples  are 
presented  to  aid  in  the  understandinq  of  the  material.  The 
technique  is  valid  for  the  general  nxn  case.  However,  for 
simplicity,  the  examples  below  are  either  single  loop,  2x2 
or  3x3  systems. 

The  flight  control  problem  involves  a  multiole 
input-multiple  output  (HIMO)  plant  requiring  regulation 
and  control  due  to  parameter  uncertainty  and  disturbances. 

The  mathematical  equations  describing  the  motion  of  an 
aircraft  are  linearized  about  a  point  in  the  flight  envelope 
or  flight  condition.  Uncertainty  arises  as  the  linearized 
coefficients  vary  with  airspeed  and  with  altitude. 

The  Quantitative  Feedback  Synthesis  Technique  devel¬ 
oped  by  Dr.  Isaac  Horowitz  uses  feedback  to  achieve  closed- 
loop  system  response  within  performance  tolerances  despite 
plant  uncertainty.  The  range  of  plant  uncertainty  and  the 
output  performance  specifications  are  quantitative  parameters 
in  the  design  process  (17;81).  The  fundamentals  of  the  design 
method  are  presented  in  the  discussion  of  the  single  input- 


single  output  design.  The  multiple  input-multiple  output 
design  procedure  is  described  in  Appendix  B,  using  the 
fundamentals  developed  in  Appendix  A. 


Problem  Definition 

The  general  single  input-single  output  (SISO)  problem 
involves  a  plant  transfer  function,  P,  with  uncertain  parame¬ 
ters  (gain,  poles,  and  zeros)  known  only  to  be  members  of 
finite  sets.  The  design  specifications  dictate  the  desired 
response  of  the  plant  to  inputs  and/or  disturbances.  The 
problem  is  to  obtain  a  controller  forcing  the  plant  output  to 
satisfy  performance  tolerances  over  the  range  of  plant  uncer¬ 


tainty  . 


The  basic  SISO  control  loop  structure  is  shown 


in  Figure  A-1 


Figure  A-1.  Two  Degrees-of-Freedom  Control  Loop 


In  this  figure,  r(t)  is  the  command  input  to  the 
system  and  d(t)  is  a  disturbance  input  to  be  attenuated. 

P  is  the  plant  transfer  function,  the  characteristics  of  which 
are  not  precisely  known.  The  compensator,  G,  and  the  pre- 


filter,  F,  are  to  be  desianeri  to  force  the  systen  output, 
y(t),  to  be  an  element  of  a  set  of  acceptable  responses, 
despite  the  uncertainty  in  P  and  the  disturbance  input  d(t). 
The  olant  input  signal,  x(t)  is  identified  since  it  is 
aenerally  of  interest  because  of  physical  constraints.  The 
signals,  r(t)  and  y(t)  are  assumed  measurable  quantities  and 
the  latter  is  available  for  feedback.  Access  to  both  signals 
allc  ;  the  use  of  the  two  degree-of-f reedom  structure  of 
Figure  A-1  and  provides  the  designer  with  two  independent 
compensator  elements,  F  and  G  (18:13).  It  is  also  assumed 
that  r(t),  y(t),  and  (for  now)  P,  such  that  y(t)  =  Px(t),  are 
all  Laplace  transformable  functions  (18:8). 

There  are  four  transfer  functions  of  interest  in 
Fiaure  A-1  where  the  loop  transmission,  L,  is  defined  as  L  = 
GP.  The  system  output  due  to  the  command  disturbance  inputs, 
respectively,  are: 


Tp  =  Y(s)/R(s) 


F  G  P 
( 1  +  G  P) 


(1  +  L) 


(A-1) 


Tp  =  Y(s)/D(s)  = 


( 1  +  G  P) 


(1  +L) 


(A-2) 


and  the  plant  input  due  to  the  command  and  disturbance 
inputs,  respectively,  are: 


Ip  =  X(s)/R(s)  = 


( 1  +  G  P) 


(1  +L) 


(A-3) 


Ij^  -  X(s)/D(s)  = 


(1  +  G  P ) 


(1  +  L) 


(A-4) 


The  desiqn  specifications  may  impose  constraints  on  any 
or  all  of  these  transfer  functions,  but  for  the  purpose 
of  this  example,  only  the  first  two  are  considered. 

Desiqn  Specif ications 

The  desian  specifications,  or  closed-loop  system 
response  tolerances,  describe  the  upper  and  lower  limits 
for  acceptable  output  response  to  a  desired  input  or  disturb¬ 
ance,  Any  output  response  between  the  two  bounds  is  assumed 
acceptable.  The  response  specifications  must  be  determined 
prior  to  applyina  the  desiqn  method.  Typically,  response 
specifications  are  qiven  in  the  time  domain,  such  as  the 
fiqures  of  merit  M^,  t^,  t^,  and  based  upon  a  step  forcinq 
function  (10:346),  or  as  a  bounded  region  as  shown  in  Figure 
A-2,  Response  to  a  step  input  is  a  good  initial  test  of 
system  response.  Bounds  (T^^)  and  (Ty)  of  the  figure  are  the 
acceptable  lower  and  upper  limits  of  a  system's  tracking 
performance  to  a  step  input.  Desired  system  response  to  a 
step  disturbance  generally  requires  maintaining  the  output 
belov/  a  qiven  value,  thus  only  an  upper  bound  is  necessary  as 
shown  by  curve  (T^)  in  Figure  A-2.  Additional  similar  bounds 
are  needed  if  other  inputs  are  to  be  considered.  The  desiqn 
technique  is  a  frequency  domain  approach?  therefore,  the  time 
domain  specifications  must  be  translated  to  bounds  in  the 
frequency  domain.  Desired  control  ratios,  =  [Y/R]|^P  and 

Tyy  =  [Y/D]j^jy.  are  modeled  to  satisfy  the  performance  specifi¬ 
cations  using  the  pole-zero  placement  method  as  described  in 


Fig.  A-2.  Time  Domain  Step  Response  Specifications 


Section  12-2  of  Reference  10.  For  response  to  a  step  input, 
a  third-order  model  with  one  zero  is  suggested. 

Tj,  =  _ ^  ^1^ _  (A-5) 

(S^  +  (S  +  P3) 

The  pole-zero  pattern  of  Equation  (A-5)  is  shown  in  Figure 
A-3.  The  locations  of  the  roots  are  adjusted  until  the  step 
response  of  the  modeled  control  ratio  matches  the  bound. 

The  frequency  domain  characteristics  are  consid¬ 
ered  during  the  response  modeling.  It  is  desirable  to 
keep  the  magnitude  difference  (as  a  function  of  frequency) 
between  the  upper  and  lower  bound  models  of  Tj^j^(j<j)  as  large 
as  possible  at  all  frequencies.  Choosing  a  lower  bound  model 
with  a  greater  pole  to  zero  ratio  than  the  upper  one  ensures 
that  the  magnitude  difference  aoproaches  infinity  in  the  limit 


Fig.  A-3.  Third-Order  Control  Ratio  Pole-Zero  Pattern 

as  (j  approaches  infinity. 

Errors  made  during  this  modeling  process  manifest 
themselves  in  one  of  two  ways.  First,  if  the  worst  acceptable 
response  model  is  not  really  acceptable,  the  system  may  not 
meet  the  specifications  over  the  assumed  range  of  uncertainty 
in  P.  And,  second,  if  the  entire  range  of  allowable  outputs 
is  not  considered,  the  bandwidth  of  the  compensation  will  be 
higher  than  necessary,  increasing  the  cost  of  the  compensator 
(18:5)  . 

Once  control  ratios  are  obtained  for  each  time 

response  bound,  a  magnitude  plot  of  the  frequency  response 

(Bode  plot)  for  each  Tj^(j<j)  is  made  on  the  same  graph  as 

shown  in  Figure  A-4.  These  plots  are  the  frequency  domain 

representation  of  the  design  specifications  on  T_  and  T_. 

R  D 

These  derived  frequency  domain  specifications  are  used  to 
obtain  the  bounds  on  the  loop  transmission,  L ( jw)  ,  as  de¬ 


scribed  later. 


Fig.  A-4.  Frequency  Domain  Specifications 
Nichols  Chart 

The  primary  tool  used  in  the  design  of  the  compensator 
elements,  G  and  F  is  the  Nichols  chart,  shown  in  Figure  A-5. 

If  the  open  loop  transmission  of  a  unity  feedback  system 
(L=GP,  assuming  F=1  for  now,  in  Figure  A-1)  is  plotted  using 
the  horizontal  and  vertical  scales  on  the  chart,  then  at  any 
given  frequency,  the  magnitude  and  phase  angle  of  T  =L/(1+L) 
can  be  read  directly  from  the  curved  scales.  Conversely,  any 
point  corresponding  to  the  magnitude  and  angle  of  Tj^  on  the 
curved  scales  provides  a  point  corresponding  to  the  magnitude 
and  angle  of  L  on  the  horizontal  and  vertical  scales  (10:332- 
334),  This  correspondence  between  L  and  T^^  on  the  Nichols 
chart  is  very  important. 

Likewise,  the  Nichols  chart  can  be  used  for  the 
disturbance  response.  Recall  that  Tp=l/(1+L).  By  way  of 
Sidi's  transformation,  L=l/m  (19:152-155)  the  system  control 


ratio  due  to  the  disturbance  becomes  Tjj=m/(l+m),  which  is  of 
the  same  form  as  Tj^=l/(1+L).  One  could  design  the  inverse  of 
the  loop  transmission,  m,  directly  on  the  Nichols  chart,  but 
it  is  much  easier  to  realize  that  by  turning  the  Nichols  chart 
upside  down,  reflecting  the  vertical  angle  of  L  lines  about 
the  -180  degree  line  (i.e.,  -190  becomes  -170,  -210  becomes 
-150,  etc)  and  reversing  the  signs  on  all  magnitude  lines,  the 
chart  can  be  used  directly  to  design  L  itself.  The  horizontal 
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and  vertical  lines  still  correspond  to  the  magnitude  and  ancle 
of  L,  and  the  curved  macnitude  lines  correspond  to  the 
magnitude  of  (1+L)  (19:155).  For  design  purposes,  only  the 
magnitude  of  (1+L)  is  reauired.  Therefore,  the  curved  angle 
lines  on  the  chart  can  be  ignored.  In  practice,  Ridi's 
transformation  is  merely  implied  by  turning  the  Nichols  chart 
upside  down  and  modifying  the  scales  as  described  above.  The 
dummy  variable,  m,  need  not  be  considered  further. 

Plant  Templates 

A  plant  template  is  a  plot  on  the  Nichols  chart  of  the 
range  of  uncertainty  in  the  plant  P  at  a  given  frequency 
(5:290).  Consider  the  example  P(s)  =  K/s(s+a)  where  the  gain 
K  is  described  by:  2  <  K  <  8,  and  the  location  of  the  second 
pole  is  given  by;  0.5  <  a  <  2.0.  An  infinite  number  of 
possible  P's  exist  due  to  the  variation  of  the  parameters,  K 
and  a;  however,  each  parameter  is  a  member  of  a  set  with 
finite  boundaries.  Likewise,  the  magnitude  and  phase  angle  of 
all  possible  P's  lie  within  finite  boundaries  when  plotted  at 
a  given  frequency.  The  plant  template  is  obtained  by  plotting 
Lm[P(j(j)]  vs.  Ang[P(j<j)]  for  all  possible  P(j(.>)'s,  that  is 
over  the  range  of  uncertainty  at  a  given  frequency  on  the 
Nichols  chart.  Note,  only  the  outer  edges  of  the  template 
need  be  calculated.  The  plant  transfer  functions  at  the 
boundaries  are  found  by  holding  one  parameter  constant  at  a 
boundary  value,  i.e.,  set  K=2,  and  vary  •a*  in  increments  from 
0.5  to  2.0.  The  frequency  response  at(*>=l  for  the  P's 
obtained  above  provides  a  set  of  points  from  A,  (K=2,  a=0.5), 
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to  D,  (K=2,  a=2),  on  the  Nichols  chart  as  shown  in  Figure 


A-5.  The  process  is  continued  until  the  complete  template 
is  formed.  For  example,  for  a=0,5,  vary  K  from  2  to  R  to 
obtain  the  line  from  A,  (a=0.5,  K=2>,  to  B,  (a=0.5,  K=8). 
Templates  are  needed  for  a  number  of  frequencies  taken  at 
regular  intervals,  such  as  every  octave.  A  set  of  tem¬ 
plates  is  shown  in  the  figure  to  demonstrate  the  change  in 
size  and  location  of  the  range  of  uncertainty  in  P  for 
different  frequencies. 

Nominal  Plant 

To  facilitate  the  shaping  of  the  loop  transmis¬ 
sion,  the  designer  needs  a  reference  or  nominal  plant  transfer 
function.  This  nominal  plant,  P^,  chosen  by  the  designer,  is 
nothing  more  than  a  reference  plant  to  be  used  in  the  defini¬ 
tion  and  shaping  of  the  nominal  loop  transmission  =  GP^. 
There  are  no  rules  or  constraints  on  the  selection  P^.  it 
doesn't  even  have  to  be  from  the  set  of  possible  P's,  but  it 
is  usually  convenient  to  choose  such  that  it  lies  at  a 
recognizable  point  on  the  templates.  It  is  convenient,  as  is 
the  case  with  the  example,  to  select  P^  such  that  it  lies  at 
the  lower,  left  hand  corner  of  the  templates.  This  choice  for 
Pq  keeps  the  bounds  on  L^,  to  be  described  next,  as  near  the 
center  of  the  Nichols  chart  as  possible.  Once  selected,  the 
Pq  point  should  be  marked  on  each  template,  as  in  Figure  A-5. 
For  the  example,  the  plant  described  by  =  2/(s  +  0.5)  is 
chosen  as  the  nominal  plant. 
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Derivation  of  Bounds  on  L 

The  system  step  response  y(t)  is  uniquely  determined 
by  the  transfer  function  T(s).  Likewise,  T(s),  for  a  stable, 
minimum  phase  system  (no  right-half-plane  poles  or  zeros),  is 
completely  specified  by  the  magnitude  of  the  frequency 
response  T(  ja>)  as  described  in  References  18  and  14,  From  the 
design  specifications,  the  frequency  response  of  the  output 
Y(  j<j)  can  vary  frjm  the  value  of  the  bound  (T^)  to  the  value 
on  the  bound  (T,  )  at  a  given  frequency  (see  Figure  A-4),  For 
the  given  example,  at  the  frequency,  a>  =  1 ,  assume  that  Y(jl) 
can  vary  from  0.7db  to  -0.8db.  The  relative  variation  in 
Y(jl)  is  (0,7)db  -  (-0.8)db  or  1.5db.  In  general,  the 
allowable  relative  change  in  Y( j  )  at  a  given  frequency  is 
expressed  as: 

4Lm  [Y(j<j.)]  =  Lm  [Ty(j(j.))  -  Lm  (TL(j(j.)]  (A-6) 

where  T^,(jtu)  and  Tj^(j(j)  are  the  frequency  domain  bounds  on 
Y(  jcj)  . 

The  relative  change  in  the  output  is  related  to 
the  control  ratio  as  follows.  From  Figure  A-1  and  Equation 
(A-1),  LmY  =  LmT  =  Lm(FL/(l+L)]  where  L  =  GP  and  it  is 
assumed  that  no  uncertainty  exists  in  G  and  F.  Then, 

L(  jw) 

4Lm  [Y(j(j.  )]  =  4Lm  [T(j<i;.  )]  =  4Lm  -  (A-7) 

[1  +  L( jw.  ) ] 

Likewise,  the  relative  change  in  L(jw)  is  equal  to  the 
relative  change  in  the  plant. 


4Lm  [L(j(j^)]  =  4Lm  [P(jcj^)l  (A-8) 

The  variation  in  P  arises  due  to  parameter  uncertainty, 
thus  the  problem  is  to  find  an  L  such  that  the  relative 
change  requirements  on  the  closed-loop  response  are  satisfied 
for  the  entire  uncertainty  range  of  P,  The  design  specifica¬ 
tions  state  the  requirements  on  the  closed-loop  response  Y(  jo;) 
and  thus  T(  jtj)  as  given  by  Equation  (A-7).  Constraints  on  the 
loop  transmission  L{  jcj)  are  desired  (5:291;  18:18). 

L ( jcj)  Rounds  on  the  Nichols  Chart 

The  relative  uncertainty  in  L  is  shown  to  be  equal  to 

the  range  of  uncertainty  in  P  by  Equation  (A-8).  As  described 

earlier,  the  plant  template  is  a  plot  on  the  Nichols  chart  of 

the  range  of  uncertainty  in  P  at  a  given  frequency.  Because 

Lm  (L)  =  Lm  (P)  +  Lm  (G)  and  also  Ang  (L)  =  Ang  (G)  +  Ang  (P), 
a  template  may  be  translated  (but  not  rotated)  horizontally  or 
vertically  on  the  Nichols  chart,  where  horizontal  and  vertical 
translations  correspond  to  the  angle  and  magnitude  require¬ 
ments  on  G(  j<j)  respectively  at  a  given  frequency  (  5:290). 
Drawing  a  line  on  each  of  the  templates  parallel  to  the 
horizontal  or  vertical  grid  lines  (see  Figure  A-5)  of  the 
Nichols  chart  is  suggested  to  maintain  correct  template 
orientation. 

With  the  template  corresponding  to  w  =  1  of  Figure 
A-5,  translate  it  to  position  1  shown  in  Figure  A-6 .  Since 
the  template  is  the  range  of  uncertaintv  in  P  and  L  =  GP, 


where  G  is  to  be  precisely  determined,  it  follows  that  the 


area  now  covered  by  the  template  corresponds  to  the  variation 
in  L  and  in  T  due  to  the  uncertainty  of  P,  Recall  the 
correspondence  between  L  and  T  on  the  Nichols  chart.  Using 
the  curved  magnitude  contours,  i.e.,  contours  of  constant 
Lm[T(jtj)],  read  the  maximum  and  minimum  values  of  T  covered  by 
the  template.  If  the  difference  between  the  maximum  and 
minimum  values  is  greater  than  the  allowable  variation  in  T  at 
the  frequency  u  =1,  that  is  {4LmtT(jl)]  as  given  by  Equation 
(A-7)  and  determined  from  Figure  A-4,  shift  the  template 
vertically,  as  shown  in  Figure  A-6,  until  the  difference  is 
equal  to  Lm[T(jl)l  (to  position  2).  Conversely,  if  the 
difference  is  less  than  that  allowed,  move  the  template 
vertically  downward  until  the  equality  is  obtained.  When  the 
position  of  the  template  achieves  the  equality  (position  2  of 
the  example)  ,  mark  the  nominal  point  of  the  template  on  the 
Nichols  chart.  The  point  marked  corresponds  to  a  bound  on  the 
magnitude  and  phase  angle  values  of  L^(jl)  read  from  the 
horizontal  and  vertical  scales  of  the  Nichols  chart,  where  the 
nominal  loop  transmission,  is  given  by; 

LQ(jw^)  =  G(  )Pq{  )  (A-9) 

Repeat  he  process  horizontally  across  the  chart  at  differ¬ 
ent  values  of  Ang(L^) .  The  points  marked  on  the  chart  form 
a  curve,  ),  representing  the  boundary  of  at 

the  given  frequency  of  the  template.  As  long  as  lies 

outside  or  above  the  boundary,  B  ( jw.  )  ,  corresponding  to  u 
=  w.  at  the  frequency  u  =  <tt  ■  ,  the  variation  in  T  due  to  the 


uncertainty  in  P  is  less  than  or  equal  to  the  relative 
change  in  T  allowed  by  the  design  specifications  at  that 
frequency.  Repeat  this  boundary,  Rn{ja»-)f  derivation  for 
various  frequencies,  using  the  corresponding  plant  tem¬ 
plates  to  obtain  a  series  of  bounds  on  (  5  :  291292). 

Likewise,  the  step  disturbance  response  specifi¬ 
cation  (line  Tp  on  Figure  A-4),  is  converted  to  bounds  on 
L^(ja»).  In  order  to  effectively  reject  the  disturbance 
the  following  inequality  must  be  satisfied: 

1/1  [1  +  L(  j(j)  ]  I  <  IC(  jtj)  I  (A-10) 

where  C(  ja>)  is  the  magnitude  of  the  boundary,  (T^),  on 
Figure  A-4.  Converting  the  magnitudes  to  decibels  and 
rearranging  terms,  the  inequality  can  be  expressed  as: 

Lm  [1  +  L(j(j)]  >  -  Lm  [C(jw)]  (A-11) 

Now  a  template  is  placed  on  the  inverted  Nichols  chart 
such  that  its  lowest  point  rests  directly  on  the  contour 
of  constant  Lm  [  1  +  L ( jo/)  ]  equal  to  -Lm  [C(jw)]  at  the 
frequency,  cj^,  for  which  the  template  is  drawn.  The  point, 
P^,  is  marked  and  the  template  slid  along  the  same  contour 
forming  a  bound,  Bjj(ju>j^)  for  L^.  Bounds  are  formed  for  each 
frequency,  in  this  manner  using  each  template.  Using  the 
rectangular  (Lm  L)  grid,  transcribe  the  bounds,  0j)(jWj^)f  on 
L^  onto  the  upright  Nichols  chart  which  already  contains 
the  command  response  bounds,  on  L^  as  shown  in  Figure 

A-7.  For  each  frequency  of  interest,  erase  the  lower  of  the 
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two  L  bounds,  where  the  remaining  bound  is  labeled,  B^(jtu). 

The  point  here  is  that  the  worst  case  (most  restrictive)  bound 
must  be  used  in  shaping  of  L^. 

Universal  Frequency  Round 

The  universal  frequency  (UF)  bound  ensures  the 
loop  transmission,  L,  has  positive  phase  and  gain  margins, 
whose  values  depend  on  the  oval  of  constant  magnitude  chosen 
(see  Figure  A-7 ) .  As  the  frequency,  a»  ,  increases,  the  plant 
templates  become  narrower  and  can  be  considered  vertical  lines 
as  u  approacnes  infinity.  The  allowable  variation  in  T 
increases  with  frequency  also.  The  result  is  the  bounds  of 
tend  to  become  a  very  narrow  region  around  the  0  db, 
-180  degree  ooint  (origin)  of  the  Nichols  chart  at  high 
frequency.  To  avoid  placing  closed-loop  poles  near  the  ja» 
axis  resulting  in  oscillatory  disturbance  response,  a  UF  bound 
is  needed  on  the  Nichols  chart.  With  increasing  w,  the  bounds 
on  approximately  follow  the  ovals  encircling  the  origin. 
Choose  one  of  the  ovals  near  the  origin.  In  Figure  A-7,  the 
contour  of  constant  magnitude  equal  to  3db  is  used  in  this 
example.  From  the  templates  corresponding  to  high  frequency, 
find  the  template  with  the  greatest  vertical  displacement,  4v, 
in  db.  Av  may  be  accurately  determined  by  finding  the 
maximum  change  in  Lm[P(jCi»)]  in  the  limit  as  w  goes  to 
infinity.  Translate  the  lower  half  of  the  3db  oval  down 
the  length  of  the  template,  i.e.,  Av  db,  as  shown,  thus 
obtaining  the  UF  bound  (see  Figure  A-7).  (Note:  Professor 
Horowitz  refers  to  this  bound  as  the  Universal  High  Frequency 


Nominal  L  (  ja».  ) 


(  UHF )  Bound )  (  18:20-22). 


Shaping  of  the  Nominal 
Loop  Transmission 

The  shaping  of  a  nominal  loop  transmission  conforming 

to  the  boundaries  of  is  a  most  crucial  step  in  the  design 

nrocess.  A  minimum  bandwidth  design  has  the  value  of  L  on 

o 

its  corresponding  bound  at  each  frequency.  In  practical 
designs,  the  goal  is  to  have  the  value  of  occurring  above 
the  corresponding  bounds,  but  as  close  as  possible  to  keep  the 
bandwidth  to  a  minimum.  Figure  A-7  shows  a  practical  design 
for  .  Note,  any  right-half-plane  (rhp)  poles  and/or  zeros 
of  must  be  included  in  to  avoid  any  attempt  to  cancel 
them  with  zeros  of  as  a  starting  point  in  the  design  of 
is  suggested,  to  avoid  any  implicit  cancellation  of  toots  in 
determining  G. 

Solving  for  G 

The  compensator,  G,  is  obtained  from  the  relation: 

G  =  If  the  found  above  does  not  contain  the  roots 

of  P^,  then  the  compensator  G  must  cancel  them.  Note,  cancel¬ 
lation  occurs  only  for  purposes  of  design  using  the  nominal 
plant  transfer  function.  In  actual  implementation,  exact 
cancellation  does  not  result  (nor  is  it  necessary)  since  P  can 
vary  over  the  entire  uncertainty  range . 

Provided  the  nominal  loop  transmission,  L^,  is 
shaped  properly,  i.e.,  meets  the  requirement  of  being  on 
or  above  the  bound,  B  (jw),  at  each  corresponding  fre- 


quency,  the  variation  in  T  resulting  from  the  uncertainty 
in  P  is  guaranteed  to  be  less  than  or  equal  to  the  allow¬ 
able  relative  change  in  T  allowed  by  the  design  specifica¬ 
tions  (5:291).  The  design  of  the  pre-filter,  F,  is  the 
final  step  in  the  design  process. 

Design  of  F 

Design  of  a  proper  only  guarantees  the  variation  in 
T( jw)  is  less  than  or  equal  to  that  allowed.  The  purpose  of 
the  pre-filter  is  to  position  Lm[T(j(»»)]  within  the  frequency 
domain  specifications.  For  the  example  given  above,  the 
magnitude  of  the  frequency  response  must  lie  within  the  bounds 
shown  in  Figure  A-4  which  are  redrawn  in  Figure  A-8 .  One 
method  for  determining  the  bounds  on  the  pre-filter,  F,  is  as 
follows.  Place  the  nominal  plant  of  the  u)  =  1  template  on  the 


Nichols  chart  where  L^(jl)  point  occurs.  Record  the  naximum 
and  minimum  values  of  Lm(T),  1.2  and  1.0  in  the  example, 
obtained  from  the  curved  magnitude  contours.  Compare  the 
values  found  above  to  the  maximum  and  minimum  values  allowed 
by  the  frequency  domain  specifications  of  Figure  A-4  at  ci>  =  1 
{0.7db  and  -O.Sdb).  Determine  the  range,  in  db,  Lm(T)  must  be 
raised  or  lowered  to  fit  within  the  bounds  of  the  specifica¬ 
tions.  For  example,  at  oi  =  1,  the  actual  Lm(T)  must  be  within 
[Lm(Ty)  =  0.7db]  >  Lm[T(jl)]  >  tLm(Tj^)  =  -O.Sdb].  But,  from 
the  plot  of  L^,  the  actual  range  of  Lm(T)  is:  1.2db  > 

Lm[T(jl)]  >  l.Odb.  To  lower  Lm[T(jl)]  from  the  actual  range 
to  the  desired  range,  the  pre-filter,  Lm(F)  is  required;  (0.7 
-  1.2  db)  >  Lm[F(jl)]  >  (-0.8  -  l.Odb),  or  -O.Sdb  >  Lm[F(jl)] 

>  -l.Sdb  (see  Figure  A-8).  The  process  is  repeated  for  each 
frequency  corresponding  to  the  templates  used  in  the  design  of 
L^.  Therefore,  in  Figure  A-9  the  difference  between  the  Ty 

and  T„^„  curves  and  the  difference  between  the  T.  and  T  . 
max  L  min 

curves  indicate  the  requirements  for  F  as  a  function  of 
frequency. 

Bounds  of  F,  [Lm(T.,)  -  Lm(T„,„)]  >  Lm(F)  >  [Lm(T,  )  - 
Lm(Tmin)]f  are  plotted  as  a  function  of  frequency  as  shown  in 
Figure  A-9.  By  use  of  the  straight  line  approximation, 
determine  a  transfer  function,  F,  such  that  its  magnitude  lies 
within  these  bounds.  The  transfer  function  obtained  in  this 
manner  is  the  pre-filter,  F  (5;301). 

The  single  loop  design  is  complete  with  the  design  of 
F.  The  system  response  is  guaranteed  to  remain  within  the 


Fig.  A-9.  Frequency  Bounds  on  the  Pre-filter,  F 


bounds  of  the  design  specifications,  provided  the  uncertainty 
in  P  stays  within  the  range  assumed  at  the  beginning  of  the 
design  process  (5;288). 

Summary 

This  appendix  presents  an  overview  of  the  SISO 
design  technique  of  Professor  Horowitz  for  single  loop  systems 
with  uncertain  plants.  The  technique  is  entirely  based  in  the 
frequency  domain,  and  makes  considerable  use  of  the  Nichols 
and  Bode  plots.  Much  of  the  designing  can  be  done  by  graphi¬ 
cal  methods. 

Design  specifications  are  translated  into  the  frequen¬ 
cy  domain  and  constitute  limits  or  boundaries  on  the  frequency 
response  of  the  system  control  ratio  and  the  loop  transmis¬ 
sion.  Two  compensator  elements,  G  and  F,  are  synthesized  to 
control  the  system  response  to  inputs  and  disturbances. 


(This  appendix  was  taken  from  Reference  1  with  minor  changes,) 


Introduct ion 

The  design  approach  for  each  loop  of  the  MIMO  system 
is  identical  to  that  for  the  SISO  system  described  in  Appendix 
A.  But  first  the  MIMO  system  must  be  separated  into  SISO 
loops  which  are  equivalent  to  the  actual  MIMO  model. 

In  general,  an  nxn  MIMO  system  can  be  represent¬ 
ed  in  matrix  notation  as  y  =  £u,  where  y  is  the  vector  of 
plant  outputs,  ^  is  the  vector  of  plant  inputs,  and  P  is 
the  plant  matrix  of  transfer  functions  relating  u  to  y. 

<  ■'  This  P  matrix  is  formed  from  either  the  linear  differential 

equations  describing  the  system  or  directly  from  the  system 
state  space  representation. 

Professor  Horowitz  has  shown,  by  the  use  of  fixed 

point  theory,  that  the  inverse  of  the  P  matrix,  referred  to  as 

'  2 
Q  ,  contains  elements  which  are  the  inverses  of  n  single  loop 

transfer  functions  equivalent  to  the  original  MIMO  plant.  The 

MIMO  problem  is  then  broken  up  into  n  loop  designs  and  n 

pre-filter/disturbance  problems,  which  are  each  handled  as 

described  in  Appendix  A  (20:677). 

The  MIMO  Plant 

Consider  the  multiple  input-multiple  output  plant  of 
Figure  B-1.  The  nxl  input  vector,  u  produces  an  nxl  output 


Fig.  B-1.  MIMO  Plant 


vector,  The  relationship  between  ^  and  u  is  described  by 

the  nxn  plant  matrix,  P,  which  is  known  only  to  be  an  element 
of  a  set  of  possible  P's.  It  is  assumed  that  the  range  of 
uncertainty  in  P  can  be  determined,  probably  in  the  form  of 
empirical  data  relating  ii  to  Note  that  the  input  and 
output  vectors  are  assumed  to  be  of  the  same  dimension. 
Although  this  may  appear  to  be  a  restrictive  assumption,  it 
can  be  shown  that  with  n  inputs,  at  most  n  outputs  can  be 
independently  controlled  (21:530-536).  Thus,  if  the  existing 
model  defines  an  unequal  number  of  inputs  and  outputs,  the 
first  step  is  to  modify  the  model  such  that  the  dimensions  of 
the  input  and  output  vectors  are  equal.  An  example  of  such  a 
modification  is  presented  in  reference  1  (1:33). 

The  plant  matrix  P,  can  be  derived  directly  from  the 
set  of  coupled,  linear,  time-invariant  differential  equations 
describing  the  behavior  of  the  plant  in  response  to  its 
inputs.  Consider  a  general  plant  model  of  the  form: 


(a)y]^  +  (bs  +  c)y2  =  (f)Uj^  +  (g)U2 

(ds)yj^  +  (e)y2  =  (h)Uj^  +  (i)U2 


(B-1) 


where  a  through  i  are  the  constant  coefficients,  Y's  are  the 
outputs,  and  the  u's  are  the  inputs  to  the  plant.  The  system 
of  Equation  (B-1)  can  be  represented  in  matrix  notation  as: 


a 

ds 


bs  +  c 


e 


f 

h 


(B-2) 


Define  the  matrix  multiplying  the  output  vector  as  M  and 
the  matrix  multiplying  the  input  vector  as  N.  The  system 
is  now  described  by: 

^  (B-3) 


The  plant  matrix  needed  is  defined  by; 

y  =  Pu  (B-4) 

Thus  the  plant  matrix,  P  is  simply; 

P  =  M~^N  (B-5) 

The  standard  state  space  representation  for  a  system 
is  described  by  the  equations  (4:93); 

X  =  ^  +  Bu  ( B-6 ) 

y  =  Cx 

The  block  diagram  for  this  system  is  shown  in  Figure  B-2. 

Although  any  number  of  states  may  be  represent¬ 
ed,  it  is  again  assumed  that  the  input  and  output  vectors, 
u  and  respectively,  are  of  the  same  dimension.  Assuming 


Fig.  B-2.  Standard  state  Space  Diagram 


the  system  is  linearized  and  the  A,  B,  and  C  matrices  are  time 
invariant,  the  plant  matrix  is: 

P  =  C[sl  -  AI'^  B  (B-7) 

This  plant  matrix  is  actually  a  representative 
member  of  a  set  of  possible  plant  matrices  due  to  the  uncer¬ 
tainty  in  the  plant  parameters,  in  practice,  a  finite  set  of 
P  matrices  are  formed  representing  the  plant  under  varying 
conditions . 

MIMO  Compensation 

The  compensation  scheme  for  the  MIMO  system  is 
similar  to  that  of  the  SISO  system  of  Appendix  A.  The 
basic  MIMO  control  structure  is  shown  in  Figure  B-3  where 
P  is  the  uncertain  plant  matrix,  G  is  a  diagonal  compensa¬ 
tor  matrix,  and  £  is  a  pre-filter  matrix.  The  design  tech¬ 
nique  can  accommodate  designs  involving  a  non-diagonal  G 
matrix  (4:14).  The  functions  of  G  and  £  are  identical  to 
those  of  G  and  F  of  the  SISO  system  of  Appendix  A.  Figure  B-4 
shows  a  more  detailed  breakdown  of  a  2x2  MIMO  system  where; 


Fig.  B-3.  MIMO  Control  Structure 
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Constraints  on  the  Plant  Matrix 

The  set  of  P  matrices  must  be  tested  to  ensure 
that  two  critical  conditions  are  met  {17:86-90); 

1.  P  must  not  be  singular  for  any  possible  combina¬ 
tion  of  plant  parameters;  i.e.,  p”^  must  exist. 

2  .  As  s  u)  t  I  I  >  IPi2  p2-|  I  for  all  possible 

plants.  This  is  the  requirement  for  a  2x2  plant. 

For  explanation  of  the  constraint  inequality  for  the  3x3 
system  see  the  second  half  of  this  appendix.  For  higher  cases 
see  Reference  17. 

The  first  condition  is  absolutely  necessary  to 
ensure  controllability  of  the  plant.  The  inverse  of  P 


Fig.  B-4.  Two-by-Two  MIMO  System 

produces  the  effective  transfer  functions  used  in  the  design. 
If  condition  2  is  not  satisfied,  it  may  be  possible  to  change 
the  ordering  of  the  input  or  output  vector  which  changes  the 
ordering  ot  the  P  matrix  elements. 

Effective  SI SO  Loops 

Now  define  a  matrix  Q*  =  p”^  which  has  elements, 

2 

The  n  effective  transfer  functions  needed  are; 

=  Reference  17  contains  the  derivation  and  proof 

of  this  equivalence.  The  nxn  MIMO  system  is  now  treated 
as  n^"  SISO  problems.  Figure  B-5  shows  the  four  effective 
SIEO  loops  resulting  from  the  2x2  MIMO  system  (20:682). 

Each  loop  in  Figure  B-5  is  handled  as  an  individual 
SISO  design  problem  in  accordance  with  the  procedures  present¬ 
ed  in  Appendix  A.  The  f's  and  g's  are  the  compensator  elements 
of  F  and  G  described  previously  The  disturbances,  d. . 

•  ^  J 


represent  the  interaction  between  the  loops. 


Fig.  B-5.  Effective  SISO  Loops 


The  in  the  above  equation  is  the  upper  response 
bound,  (Ty  or  in  Figure  A-4),  for  the  respective  input/out¬ 
put  relationship.  These  are  obtained  from  the  design  specifi¬ 
cations  (20:681-684).  Note  that  the  first  digit  of  the 
subscript  of  bj^^  refers  to  the  output  and  the  second  digit  to 
the  input.  Thus,  bj^^  is  a  function  of  the  response  require¬ 
ments  on  the  output,  yj^ ,  due  to  the  input,  r^. 

A  recent  improvement  in  the  design  technique  involves 
modification  of  the  q's  on  the  second  and  subsequent  loops 
based  on  the  g's  already  designed.  This  reduces  the  overde¬ 
sign  inherent  in  the  early  part  of  the  design  process.  During 
the  design  of  the  final  loop  the  exact  equation,  representing 
the  loop  and  the  interactions  of  the  other  loops,  is  used 
(22:977).  The  use  of  this  improvement  is  demonstrated  in  the 
original  continuous  domain  designs  used  in  this  thesis 


investigation  (1:33). 


Basically  Non-interacting  ( BNIC)  Loops 

When  the  response  of  an  output,  ,  due  to  an  input, 
tj,  is  ideally  zero,  the  yj^j  loop  is  called  a  basically 
non-interacting  (BNIC)  loop  (20:679).  Due  to  loop  interaction 
and  plant  uncertainty,  this  ideal  response  is  not  achievable. 
Therefore,  the  performance  specifications  describe  maximum 
responses  and  the  loop  is  handled  exclusively  as  a  disturbance 
rejection  problem. 

3x3  MIMO  System  Design  Theory 

The  remainder  of  this  appendix  is  devoted  to  the  3x3 
MIMO  system  structure  with  extension  to  a  general  nxn  system. 
Several  methods  are  discussed  of  ways  to  reduce  the  complexity 
required  in  the  design  process.  Finally,  a  discussion  of  the 
improved  design  technique  is  highlighted. 

Control  Structure 

The  system  structure  is  similar  to  that  of  the 
2x2  system  discussed  previously.  The  basic  structure  is  shown 
in  Figure  B-6  where  P  is  the  uncertain  plant  matrix,  G  is  the 
compensator  matrix,  and  F  is  the  pre-filter  matrix.  The 
purpose  of  G  and  F  is  identical  to  those  of  the  SISO  system. 
Figure  B-7  shows  the  3x3  MIMO  structure  with  G  being  diagonal. 
The  general  form  of  the  aircraft  equations  of  motion,  for  a 
3x3  MIMO  system,  are  of  the  form: 


Fig.  B-6.  The  3x3  Multiple  Input- 
Multiple  Output  Feedback  Structure 

(as  +  b)Y^  +  (cS  +  d)Y2  +  (eS  +  f)Y3 

=  gu^  +  hU^  +  103  (B-9) 

(jS  +  k)Y^  +  (IS  +  in)Y2  +  (nS  +  o)Y3 

=  PU^  =  gU2  +  rD3  (B-10) 

(ts  +  v)Y^  +  (wS  +  x)Y2  +  (aas  +  bb)Y3 

=  ccUj^  +  dd02  +  eeU3  (B-11) 

where  a  through  ee  are  constant  coefficients.  Equations 
(B-9)  through  (B-11)  can  be  written  in  matrix  notation  as; 

as  +  b  cS  +  d 

jS  +  k  IS  +  m 

tS  +  V  WS  +  x 

Equation  (B-12)  is  of  the  form; 

Mv  =  Nu  (B-13) 


eS 
nS  +  o 
aas  +  bb 


.  fl  ^y'^ 


1  , 


3b 


g  h  i 

i  P  q  t 
cc  dd  ee 


I 


lUij 


(B-12) 


L3 


where  the  output  matrix  M  multiplies  the  output  vector  ^ 
and  the  input  matrix  N  multiplies  the  input  vector  u. 

Equation  (B-13)  is  manipulated  to 

^  =  M~^Nu  (B-14) 

where  £  =  M~^N  (B-15) 

is  defined  as  the  plant  matrix.  Thus  equation  (B-14)  becomes 

^  =  Pu  ( B-16 ) 


Fig.  3-7.  The  Complete  3x3  MIMO 
Structure  with  Diagonal  G 


The  remaining  elements  of  the  3x3  MIMO  system  can  be  repre¬ 
sented  in  matrix  form  as: 
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In  this  thesis  G  is  always  a  diagonal  matrix  but  G  can 
be  a  full  matrix  which  gives  the  designer  more  flexibility 
in  the  design  technique. 

Constraints  on  the  Plant  Matrix 

There  are  two  constraints  on  the  P  matrix  which 
must  be  satisfied.  These  constraints  are: 

1.  The  plant  matrix  P  must  be  invertible.  That  is, 
P"^  ,  must  exist  (i.e.,  P  must  not  be  singular). 

2.  As  55  — ►  CO 

'^llP22P33'  ^  'P11P23P21'  'P12P21P33' 

+  'P12P23P31’  *^13^22^31'  'P13P21P32' 

for  all  possible  plants  (20). 

The  first  constraint  ensures  that  the  plant  is  control¬ 
lable.  If  constraint  2  is  not  satisfied,  then  the  original 
order  of  the  input  and  output  matrices  may  be  changed  to 
ensure  that  this  constraint  is  met.  It  is  to  be  noted 
that  constraint  2  is  not  required  for  the  improved  design 
method.  The  only  requirement  is  that  the  diagonal  terms 


do  not  change  sign  as  S 
where  i= j ) . 


through 


00  (i.s.f  ^22 ' 


Effective  SISO  Loops 


First  define  a  matrix  Q'=p  ^  having  elements  q!j  and 

2 


a  matrix  Q=[qj^j]  where  Thus  q^^  are  the  n 

effective  transfer  functions  that  are  needed  for  the  QFT 


technique.  Using  the  above  transformation  the  3x3  system  is 
2 


now  treated  as  n  (nine)  SISO  problems  (20).  Figure  B-8  shows 
the  resulting  nine  effective  SISO  loops.  Note  that  the  first 


Q  = 


'ii 

and  y. j 

refers 

to  the  ith  output 

and 

the 

to  the  jth  input. 

Since  qij=l/q' j,  Q 

can 

be 

■  tr  ix 

form  as 

'll 

^12 

^13 

'21 

‘^22 

^23 

(B- 

18) 

‘31 

^32 

^533 

Using  Figure  B-8  the  input/output  relationship  is  defined 
as : 


t .  .  =  y .  . /r 

13  ^13''  ; 


(B-19) 


r  .  = 


where 


1 


Fig.  B-8.  The  Final  Equivalent  Single  Loop 
Feedback  Structure  which  Replaces  the  3x3  MIMO  System 


Therefore  the  final  equations  for  the  nine  SISO  effective 
loops  are; 


hi 

= 

^ii^i'^ii 

+ 

‘^iihi 

(B-20) 

1  +  9iqii 

1  +  g^qii 

*^12 

+ 

‘^11^12 

(B-20a) 

1  +  qiqii 

1  +  g^qii 

*^13 

+ 

^ll'^13 

(B-20b) 

1  +  g^qii 

1  +  g^qii 

^21 

^21*^2^22 

+ 

^22^^21 

(B-20c) 

1  +  *32^22 

1  +  g2q22 

*^22 

^22^2^522 

+ 

*^22^^22 

(B-20d) 

^  ■*■  ^2*^22 

1  +  52^522 

*^23 

^23^2'^22 

+ 

‘522'^23 

(B-20e) 

^  ^2^22 

1  +  g2‘?22 

*^31 

^31^3^33 

+ 

^33^^31 

(B-20f ) 

1  +  g3q33 

1  +  g3q33 

*^32 

^32^3'^33 

+ 

^33*^32 

(B-20g) 

^  ■*■  ^3'^33 

1  +  ^3^333 

^^33 

^33^3*^33 

+ 

'?33‘^33 

(B-20h) 

1  +  g3q33 

1  +  ^3^33 

where : 

hi 

= 

[t2l/^12 

+  t 

3l/'^l3^ 

(B-21) 

(B-21a) 


°13  "  ^^23^^12  *^33/^13^ 


‘^21  ■  ^  '^31^^23^ 


^^22  ^^12/'^21  ■*■  *^32^^23^ 


*^23  ^  ■  ^'^13/^21  '*'  *^33/^23^ 

^^31  "  "  ^’^11/^31  ^2l/*^32^ 

^32  "  ■  ^*^12^^31  ^22^^32^ 

^33  ■  ^*^13/^31  ■*■  ^33/^32^ 


(B-21b) 

(B-21C) 

(B-21d) 

(B-21e) 

(B-21f ) 

(B-21g) 

(B-21h) 


Finally,  the  above  equations  can  be  written  in  a  general 
2 

form  for  any  n  equivalent  single  loop  feedback  structure 
(20:680).  Thus  in  general  form: 


^ir 


,  j  rf 

i1^i~^ii  ii^ii 

1  +  9  •  q  •  • 
’1^11 


where:  d^^^  =  “  “k  '^kj^'^ik  k  /  j 

Equation  (B-22)  can  be  expressed  as 
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fj^L. 

1  +  L, 


{B-22) 


(B-23) 


where;  L. 


a  .  a .  . 
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Simplification  of  the 
Single-Loop  Structure 

The  single-loop  structure  in  Figure  B-8  has  two 
components,  one  due  to  the  input  r^  (assume  rj=l)  and  one 
due  to  the  disturbance  input  Therefore,  the  control 

ratio  t^j  is  the  sum  of  control  ratios  each  involving  the 
jth  input  and  can  be  written: 


"^i  j 


f.  .L. 

1  +  L. 


d  .  .q  •  . 


1  + 
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(B-24) 

(B-25) 


For  a  fixed  £  which  is  an  element  of  the  universal  set  P, 
'‘^ij'max 

'^ij^max  “  k/i  ^  *'k  *  max 


Since  the  relative  phases  of  are  not  known,  the 

design  process  must  use  the  extreme  cases; 
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k;^i 


max 
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(B-27) 


The  term  b^^  is  the  control  ratio  that  represents  the  maximum 
allowable  magnitude  of  the  output  due  to  the  disturbance  input 
which  is  one  of  the  required  design  specifications. 

Performance  Tolerances 


The  performance  tolerances  are  divided  into  two 
separate  portions.  One  set  of  tolerances  for  the  inter- 


acting  loop  (tracking  loop)  and  a  second  tolerance  for  the 


Basically  Non-Interacting  (BNIC)  loop  (Appendix  A).  Thus, 
for  the  BNIC  loop  it  is  necessary  that: 


rij(j 


)  + 


■d  .  .  < 
11 


max  — 


I  <  It 


1 3  max 


+  It, 


1 3  max 


=  It,  I  <  b.  .(u»)  (B-28) 

ij  max 

The  above  equation  can  be  justified  since  the  relative 
phase  of  the  two  terms  (T^.(ju))  and  t,  ( ju))  )  are  not 
known.  Then  T^j(ja))  can  be  forced  to  zero  by  simply  selecting 
fj^j,  for  i  =  j,  to  be  zero.  b^j(ju>)  are  the  upper  bounds  for 
the  design  specification  on  the  disturbance  rejection 
(20:681).  The  performance  tolerance  for  the  interacting  loop 
is : 


li^(w)  <  It.  .(jcj)  +  T,  (j{*>)l  <  b.  . 
J  J  ij  j 


(u>) 


(B-29) 


where  a^^  and  b^^j  are  the  upper  and  lower  bounds  of  the 
design  specifications  for  the  interacting  loop.  Applying 
the  above  tolerances  to  the  BNIC  loop  results  in  the  following 
inequality: 


T.  .  I  = 
11 


‘11 


Id 
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1  +  L.| 


<  "^ij 


(B-30) 


where  for  the  3x3  system 


Idi-I  =  lb.  ./q..  +  b.  ./q..  I  k/i 


Rearranging  Equation  (B-30)  yields  the  following  inequal¬ 
ity: 


1  +  L.|  ^ 


^ki/^il 


Iqii 


(B-31) 


Equation  (B-31)  is  used  in  the  design  of  all  BNIC  loops, 
which  generates  the  required  bounds  on  the  loop  transmis¬ 
sion  for  a  disturbance  input. 

The  interacting  loop  tolerance  generates  the  following 
conditions  which  must  be  satisfied  in  the  design: 


•rij  =  iijhi  •  (B-32) 

1+L.  ^ij  1+L. 

where  the  bounds  on  are  found  exactly  as  outlined 

i  j 

above  and  the  bounds  on  are  determined  using  the  plant 
uncertainty  (templates)  and  the  given  tolerances  and 
b^j  on  the  interacting  loop.  It  must  be  noted  that  the 
tolerances  given  by  and  b^^  are  for  the  bounds  on  the 
entire  interacting  loop.  Therefore,  the  tolerances  must 
be  divided  between  the  interacting  and  BNIC  cases.  If  one 
of  the  equations  of  Equation  (B-32)  is  much  more  dominant 
than  the  other,  then  the  entire  bounds  can  be  dedicated  to 
that  particular  portion.  The  designer  is  cautioned  that 
if  one  of  or  dominates,  it  must  be  verified  over 

the  entire  frequency  range  of  interest. 


Improved  Design  Technique 

The  improved  technique  reduces  the  inherent  overde¬ 


sign.  The  unimproved  technique  does  not  exploit  the  correla- 


tion  between  the  t^^  of  the  system.  With  the  improved 
technique  this  correlation  between  the  is  taken  into 

account  for  the  second  and  subsequent  loop  design.  This 
technique  is  highlighted  using  a  3x3  MIMO  system  as  an 
example . 

Assume  that  loop  one  is  the  first  loop  to  be  designed. 
Therefore  the  general  equation  for  these  three  SISO  problems 
is  obtained  from  Equation  (B-22)  (22:977-988)  and  is  given  as; 


Li  =  giq 


1^11 


3 

=  L 
k/l 


tij/q] 


(B-33) 


Using  Equation  (B-33)  the  elements  and  fj^^  are  designed 
such  that  are  stable  and  meet  the  desired  tolerances 

or  specifications.  With  loop  one  design  completed,  the 
second  loop  design  can  now  be  accomplished.  The  general 
equations  for  loop  two  are  obtained  using  Equation  (B-22). 
The  resulting  equation  is; 


•2j^2  ^  ^2i^22 


1  +  L, 


(B-34) 


d,.  =  L  ^2j/^2k  j  =  1,2,3 

By  solving  for  t^^^  in  Equation  (B-33)  and  substituting 
these  into  Equation  (B-34)  it  can  be  shown  (22)  that 


I 


a 


I 

I 


I 


i 

i 


V  ■  ;  *  T  «»  ■'  ^  V  ■  W  *  i  ■  v^'  V  "^  I  ^  »~T  ’  * 


'2  j 


^2j^2e  *^26 


1  +  T: 


26 


where  L,  =  ^2*^22 _ 

y,,  =  ^11^22 

^12^21 


‘22e 


2j 


_  q--  (1  +  L,  ) 

-  ,i  .  4.  - 


=  ^2e 


‘11 


y‘^21‘^13^^  ^1^ 


q2i(l  +  L^) 


(B-35) 


(B-35a) 


(B-35b) 


(B-35c) 


(B-35d) 


Using  equations  (B-35a)  through  (B-35d)  the  elements  L2g 

and  f_ .  can  be  designed  to  meet  the  desired  tolerances 
^  J 

since  Lj^  and  are  known  from  the  first  loop  design. 

The  final  loop  is  accomplished  in  a  similar  manner  as 
the  second  loop  design.  The  resulting  equations  are; 


A  =  ^23(1  +  L^)  +  >'13(1  +  L2) 


13 


'3i 


^’'l2‘^l  ^13^2’ 

(B-36e) 

=  (1  +  L^)(l 

■^^2^  ~  ^12 

(B-36f ) 

=  ^21^33  ; 

^23^^31 

u.,  =  ‘531^22 
‘^32*^21 

(B-36q) 

=  ^11^22 
‘^12‘J21 

=  ^^22^33 
‘^23^32 

(B-36h) 

=  ^11^33 
<^13^31 

(B-36i) 

=  ^li’"l‘533"l 

^2i^2‘^33'^2 

(B-36j) 

£ 

-  A 

=  ^22 
'^2l‘^31 

(1  +  L2)/q3i 

(B-36k) 

= 


‘11 


^12^31 


-  (1  +  Lj^)/q 


32 


(B-361) 


In  Equations  (B-36a)  through  (B-361)  L^,  L2f  ^2j' 

all  known;  therefore,  the  only  unknowns  are  and  g^. 
These  now  constitute  single-loop  uncertainty  problems  and 
can  be  designed  using  the  methods  outlined  in  Appendix  A. 
The  required  elements  are  designed  to  the  desired  toler¬ 
ances  or  specifications.  This  completes  the  improved  3x3 
MIMO  design  technique.  If  all  loops  are  designed  to  meet 
the  given  tolerances,  then  the  desired  response  from  the 
MIMO  system  is  guaranteed. 
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Equilibrium  exists  when  it  is  impossible  to  reduce  the 


burden  on  any  L^,  without  increasing  it  on  some  other 
(20:683).  This  simply  states  that  some  of  the  design  toler¬ 
ances  may  be  decreased  or  increased  for  certain  loops  without 
making  it  more  difficult  on  other  loops.  This  results  in  only 
one  column  of  the  equivalent  SISO  systems  being  dominant. 
However,  after  equilibrium  is  reached,  it  may  be  desirable  to 
sacrifice  one  loop  for  the  sake  of  another.  This  involves 
"tradeoffs"  between  the  different  loops.  it  should  be  noted 
that  these  tradeoffs,  when  used,  always  make  it  harder  to  meet 
specifications  on  one  loop  when  reduction  is  accomplished  in 
another . 

Summary 

This  Appendix  describes  the  design  technique  used  in 
the  design  of  single  input-single  output  and  multiple  input- 
multiple  output  systems.  First,  the  SISO  system  technique  is 
outlined  and,  second,  the  technique  is  expanded  to  the  MIMO 
system. 

The  appendix  describes  the  multiple  input-multiple 
output  plant  and  the  corresponding  plant  matrix.  Guidelines 
are  presented  for  finding  the  P  matrix,  which  relates  the 
input  vector  to  the  output  vector. 

The  division  of  the  MIMO  system  into  separate  SISO 
loops  is  presented  using  the  inverse  of  the  P  matrix. 

After  the  equivalent  SISO  loops  are  determined,  each  is 
designed  in  accordance  with  the  SISO  design  theory  pre- 


Appendix  C:  Plant  Transfer  Functions 


Introduction 

Appendix  C  contains  the  input-output  matrices  for  the 
KC-135  F.C.'s  2  and  3  for  both  the  lateral  and  longitudinal 
modes.  The  individual  plant  transfer  functions  for  F.C.'s  2 
and  3  of  the  KC-135  lateral  mode  are  also  listed.  The 
individual  plant  transfer  functions  for  all  three  F.C.'s  of 
the  KC-135  longitudinal  mode  are  presented.  The  plant 
transfer  funcitons  for  all  four  of  the  F-16  flight  conditions 
are  listed  in  this  appendix  as  well. 


KC-135  Lateral  Mode  Input-Output  Matrices 
Flight  Condition  ^ 


(S^+.920643S-. 016597) 


{ .3974S+6. 25223) 


(  .  03994  58S^+.06167S+.00821)  -( S^+ . 27271 9S+1 . 88298  : 


.46223 

1.111 

( .0007655S+. 02546) 

-( .0273552S+1.3286) 

«r 

_ 

Flight  Condition  2 


(S^+1.09665S-. 09548) 


(  .  69382S+2.4876) 


(  .  0  316455^+. 228 35S+. 0  3271)  -( S^+ . 368629S+ . 609055 ! 


.37757 

(  .0024207S+. 01926) 


.36189 

-( .03819S+. 46618) 


’  w 


0 

/S 


KC-135  Longitudinal  Mode  Input-Output  Matrices 


Rigid  Body 

Flight  Condition  2 


L041S 

( .56013-1.113) 

(3+. 0029) 

(-. 07433^-. 03555S) 

( .00633+. 479) 

.0049 

1^.17893 

(S^+.8149S+2.4608) 

-.0091 

.010498 

-.0798  .034 

6 

e 

= 

-.0186 

.0273  0 

-1.758 

.218  0 

^  T 

Flight  Condition 

3 

ri031S 

( .166S-.9848) 

(S+.0234) 

(-.24383^-. 18993) 

( .03353+. 7784) 

.0592 

-.26193 

(3^+1.0923+1.074) 

-.0318 

h 

9 

u 


h 

0 

u 


.0065 

-.0631 

-.0375 

.0646 

-.9218 

.0954 

.1268 

0 

0 


6 

e 

«sb 

1 

o. 

►9 

Elastic  Body 

The  general  form  of  the  elastic  body  input-ouput 


(C-ll) 


C-12) 


matrices  is  given  in  Equation  (c-13).  The  individual  elements 
(a  through  f)  of  the  input-output  matrices  are  given  in 
Equations  (C-14)  through  (C-16). 


Flight  Condition  3^ 

a  =  -.24385^  -  (.00122/A  +  .00059/B)S  (C-16a) 
b  =  (1  +  .00401/A  -  .00517/B)S  +  .00499/A  +  .00244/B  (C-16b) 
c  =  (-.00366/A  +  .00178/B)S  (C-16c) 
d  =  +  (.00121/A  -  .01548/B)S  +  .0015/A  +  .00731/B  (C-16d) 
e  =  .000961/A  +  .00125/B  (C-16e) 
f  =  .00029/A  +  .00374/B  (C-16f) 
A  =  +  1.893S  +  131.014  =  (S  +  .9465  +jll.4069)  (C-16g) 
B  =  +  1.229S  +  547.77  =  (S  +  .6145  +j23.3964)  (C-16h) 


KC-135  Longitudinal  Mode  Plant  Transfer  Functions 


Rigid  Body 


Flight  Condition  1 


p  _  . 277230(S+. 0136689) (S-6. 15013) (S+6. 35439) 
^11  ~  -  - 


S (S-. 0568942 ) (S+. 0599089 ) (S+. 647974+ jl. 66670) 


p  _  ^  416125(5-. 0731283+j. 0712002) (S+. 877711) 

S ( S-. 0568942 ) (S+. 0599089 ) (S+. 6 47974+ jl. 66670) 


13 


.  00 315 410 ( S+. 384 350+ jl. 95325) 


S(S-.0568942) (S+.0599089) (S+. 647974+ jl. 66670) 


21 


_  -1. 64900(5+. 00813148 ) (S+. 500895) 


( S-. 0568 942 ) (5+. 0599089 ) (5+. 647974+ jl. 66670) 


22 


.  1733 ( 5-. 06 33967 ) (5+. 08715  44) 


( 5-. 0568942 ) (5+. 0599089 ) (S+. 647974+ jl. 66670) 


23 


.00051975(5+1.84759) 


( 5-. 0568942 ) (5+. 0599089 ) (5+. 647974+ jl. 66670) 


(C-17) 


(C-18) 


(C-19) 


(C-20) 


(C-21) 


(C-22) 


.011619(S+.294472) (S-1.75751) (S+77.6744) 

(  S-. 0568942)  (S+. 0599089)  (S+. 64797  4+ jl. 6667  0~) 


{C-23) 


0827 (S-. 0389887) ( S+ . 486976 )( S  +  1 . 76180 ) 
(S-.0568942) (S+. 0599089 ) (S+. 6 4797 4+ jl. 66670) 


{C-24) 


.0495(5-. 00 163560) ( S+ . 648795+ jl . 66 3 52 ) 

(  S-. 0568942  ) (S+. 0599089 ) (S+. 647974  + jl. 66670) 


{C-25) 


FI iqlit  Condition  2 


.  250336(5+. 01 50531 ) (S-6. 41871) (5+6.62877) 
S(5-.0499167) (S+. 0625482 ) (S+. 6 41481+ j 1.55626) 


(C-26) 


-.  367429(5+. 00725990 ) (S-. 828748) (5+1.603  30) 

5 ( 5-. 0499167 ) (5+. 0625482 ) (5+. 641481+ j 1.55626) 


(C-27) 


.  00224226 ( 5+. 413 300  + jl. 78313) 

5 (5-.  0499167 ) (S+. 0625482 ) (5+. 641481  + jl.  55626) 


(C-28) 


-1. 758 (S+. 00890 135 ) (S+. 446262) 

(  5-. 0499 167 ) (5+. 0625482 ) (5+. 6 41481+ j’  .55626) 


(C-29) 


.  218 (5-. 0110571 ) (5+. 186891) 

(5-.0499167) (5+. 0625482 ) (S+. 641481+ jl. 55626) 


(C-30) 


.0003094(5+1.77430) 

(  5-. 0499167 ) (5+. 0625482 ) (S+. 6 41481+ jl. 55626) 


(C-31) 


.010498(5+. 256483) (5-1 .73052) (5+95.5837) 

( 5-. 0499167 ) (5+. 0625482 ) (S+. 641481+ jl. 55626) 


(C-32) 


-. 0798 (5-. 308874 ) (5+. 315440) (5+2.62745) 

( 5-. 0499167 ) (S+. 0625482 ) (S+. 641481+ j 1.55626) 


■J  1  \ 


(C-34) 


. 03 4 ( S+. 00793446 ) (S+. 64 2379+ jl. 55311) 
(S-.0499167) (S+. 0625482 ) (S+. 641481+ jl. 55626) 


Flight  Condition  ^ 

.  153815 (S+.  105056  ) (S-4. 12352 ) (S+4. 3 206 5) 
S(S-. 136952) (S+. 192094) ( S+ . 9 19587+ jl . 0 2075 ) 

264971  ( S-. 148416 ) (S+. 208363 ) (S+1. 06381) 
i( S-. 136952) (S+. 192094) ( S+ . 919587+ jl . 0 2075 ) 

.  030790 (S+. 554997+ jl. 08817) 

S (S-. 136952 ) (S+. 192094 ) (S+, 919587+ jl. 02075) 

-. 9218 (S+. 0604067 ) (S+. 697985) 
rs-.136952) (S+. 192094 ) (S+. 919587+ j 1.0 2075) 

. 0954 (S-. 127425 ) {S+. 181287) 

( S-. 136952 ) (S+. 192094 ) (S+. 919587+ jl. 02075) 

.00403224(S+2. 77876) 

(S-.136952) (S+. 192094 ) (S+. 919587+ jl. 02075) 

.0065(S+.446838) (S-4.80063) (S+27.3263) 

( S-. 136952 ) {S+. 192094 ) (S+. 919587+ jl. 02075) 

06 31(S-. 164207 ) (S+. 156521) (5+1.69664) 

( S-. 136952 ) (S+. 192094 ) (S+. 919587+ jl. 02075) 

.1268(5+. 000194722) ( S+ . 935361+ jl . 00664 ) 

(~5-.  136952 )  (5+.  192094  )  (5+.  919587+ jl.  02075) 


(C-35) 

(C-36) 

(C-37) 

(C-38) 

(C-39) 

(C-40) 

(C-41) 

(C-42) 

(C-43) 


Elastic  Body 

Flight  Condition  1^ 


ell 


-.  017  2302 ( S+. 00001 2  507 ) (S-3. 11866) 

S (S-. 0006 29852 ) (S+. 000309291+ j. 000362912) 

(S-.690542+jl7.3342) 

( S+. 659000+ jll. 5603 ) (S+. 556000+ j23. 3392) 


(C-44) 


;av.-- 

i.  t  m.t 

t 

*.  n’.  ■v'’ . 


V  _  ■*, 

f 


*•*  *“■."  % 


f. 

» 


e21 


007590(S+. 0000095) 


( S-. 0006 29852 ) (S+. 000309291+3.000362912) 


(S+.651531+jl2. 783921) 


cll 


(S+.659000+jll.5603) (S+. 556000+ j23. 3392) 


.  277230 (S+. 0000038 ) (S+. 0146521 ) (S+6. 35534) 
S (S+.O 599089 ) (S-. 000629852 ) (S-. 0568942) 


( S-6. 15050 ) (S-. 000531746+ j. 00384070) 


(C-45) 


c21 


( S+. 000309291+ j. 000362912 ) (S+. 647974+ j 1.66670) 


-1. 64900 (S+. 0000043 ) (S+. 0102544 ) (S+. 500952) 
(S+.0599089) (S-. 000629852 ) (S-. 0568942) 


«  (C-46) 


( S-. 00110205+ j. 00453225) 


(C-47) 


(  S+.  000309291+ j. 000362912 ) {S+. 647974+ j 1.66670) 


Flight  Condition  2 


ell 


-.0262853(5+1.13826) (S-. 0000215) 

S ( S+. 0000087 ) (S+. 924000+ j9. 60782) 

(S-3. 80642+j9. 41610) 


( S+. 598000+ j22. 7344 ) (S-. 00031041+ j. 0154348) 


(C-48) 


t  .  .-  4 


w.  ,.  J 


iri 


t  T 
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•  •  •  "ji 


'  •  o‘ 

P 

;.S/y 


.0162(5+4.90139) (S-. 00024978) 

(  S+. 0000087 ) (S+. 924000+ j9. 60782) 


(S-2. 88813) 

(S+. 598000+ j 22. 7344) (S-. 00031041+ j. 01543 48) 


(C-49) 


. 250336(5+. 0149432) (5+6.62950) (S-6. 41928) 

5 ( S+. 06 25482 ) (S-. 0499167 ) (5+. 641481+ jl. 5 5626) 

(5-. 00024721+j. 0153769) 

(S-.00031041+j .0154348) 

-1. 75800 ( 5+. 446247 ) (S+. 00872367) 

( 5+. 06 25482) (5-. 0499167 ) (S+. 641481+ jl. 55626) 

(5-. 00021654+j. 0153884) 

(5-. 00031041+j. 0154348) 


Flight  Condition  2 

-.00906 89(5+. 0000283) (5-. 948309) 

5(5+.0000425) (5+. 00873567 ) (S-. 00875316) 

(S+.364981+jl7.5884) _ 

( 5+, 946 518+ jll. 4069 ) (5+. 614516+ j23. 3964) 

.00403(5-. 0000651) _ 

(5+.0000425) (S+. 00873567) (5-. 00875316) 

(5+.922475+jl2.6551) 

( 5+. 946518+ jll. 4069 ) (5+. 614516+ j23. 3964) 

.  153815(5+. 00004381 ) (S+. 00832179 ) (S+. 10510 5) _ 

5 (5+, 00004252 ) (5+. 00873567 ) (5+. 192094 ) (S-. 00875316) 


(5+4.32110) (5-. 00841569) (5-4.12404) 
( 5-. 136952 ) (5+. 919587+ jl. 02075) 


(C-54) 


921800  (S-t-.  00006170  )  (S+.  00791259  )  (S+.  06060  49) 

( S+. 00004252 ) (S+. 00873 567) (S+. 192094 ) (S-. 00875316) 

(S+. 697999) (S-. 00817314)  ( 

( S-. 136952 ) (S+. 919587+ jl. 02075) 


F-16  Plant  Transfer  Functions 
Note  that  due  to  symmetry, 

Flight  Condition  ^ 


^  ^  -1, 118S ( S+. 0 1822 ) (S+. 4568) 

(S-.  3633  )  (S+1. 3)  (S+.  07683+ j.  2065) 

q_  ^  . 1209S ( S+. 06537 ) (S+. 2589) 

6^  (S-.3633) (S+1. 3) (S+. 07683+ j. 2065) 


p_  ^  -2.239S(S+.205+j.853) _ 

d®  (S+.1041) (S+. 6835) (S+. 2741+ jl. 909) 


p_  ^  -2.142S(S+.3017+jl.562) _ 

i®  (S+.1041) (S+. 6835) (S+. 2741+ jl. 909) 


Flight  Condition  2 

pa  ^  -2. 931S ( S+. 01004 )(S+. 5502) 

(S-1.167) (S+2. 028) (S+. 006472+ j. 07803) 

pa  ^  -.1059S(S+. 006697)  (S+1. 861) 

(S-1.167) (S+2. 028) (S+. 00647 2+ j. 07803) 


P 


a 

22  “ 


(C-56) 

(C-57) 

(C-58) 

(C-59) 

(C-60) 

(C-61) 


-8.723S(S+.219+jl.607) _ 

(S+.07795) (S+. 8265) (S+. 211+ jl. 953) 


(C-62) 


-6 .792S(S+.2442+j2. 101 ) 
(S+.07795) (S+. 8265) (S+. 211+ jl.  953  ) 


(C-63) 


Flight  Condition  2 

-12.03S(S+.01262) (S+1.51) 

(S-.9645) (S+3. 223) (S+. 00755 3+ j, 05384) 

-3. 236s (S+. 01254 ) (S+1. 646) 

(S-.9645) (S+3. 223) (S+. 007553+ j. 05384) 

-25.53S(S+.3541+j2.927) 

( S+. 02719 ) (S+2. 697) (S+. 391+ j2. 962) 

-25.36S(S+.3749+j3.578) 

(S+.02719) (S+2.697) {S+.391+j2.962) 


Flight  Condition  £ 

-16. 45S(S+. 02996) (S+1. 097) 
(S+.01516+j .02343) (S+. 8012+ j6. 592) 

-2. 925S ( S+. 03459 ) (S+. 6861) 
(S+.01516+j.02343) (S+. 8012+ j6. 592) 

-7.084S(S+.408  3+j4.916) _ 

(S+.03448) (S+2.171) (S+.4996+j3.129) 

-23.3S(S+.3774+j3.848) 
(S+.03448) (S+2.171) (S+. 4996+ j3. 129) 


(C-68) 


(C-69) 


(C-70) 


(C-71) 


Sideslip  Command 

The  upper  bound  for  the  sideslip  response  to  a  5  degree 
sideslip  command  is  shown  in  Table  0-4. 


TABLE  D-4 

b-_ — Time  Domain  Specifications 
KC-135  Lateral  Mode 


Rise  time  =  4.91  seconds 
Settling  time  =  8.90  seconds 
Peak  value  =  5.0  degrees 
Final  value  =  5.0  degrees 


Table  D-5  outlines  the  worst  acceptable  response  or 
lower  bound  for  sideslip  response  with  a  5  degree  sideslip 
command . 


TABLE  D-5 

^22 — Time  Domain  Specifications 
KC-135  Lateral  Mode 


Rise  time  =  7,87  seconds 
Settling  time  =  14.4  seconds 
Peak  value  =  5.0 


Final  value  =  5.0 


Table  D-6  shows  the  maximum  allowable  time  response 
characteristics  for  bank  angle  response  to  a  sideslip  command 
of  5  degrees. 


TABLE  D-6 

b_. — Time  Domain  Specifications 
KC-135  Lateral  Mode 


Peak  time  =  .681  seconds 
Settling  time  =  20  seconds 
Peak  value  =  2.0 
Final  value  =  0.0 


KC-135  Longitudinal  Mode  Response  Models 
Pitch  Angle  Command 

Table  D-7  shows  the  upper  bound  pitch  angle  response 
for  a  4  degree  pitch  angle  command  input. 

TABLE  D-7 

b-- — Time  Domain  Soecif ications 
^KC-135  Longitudinal  Mode 


Rise  time  =  2.93  seconds 
Settling  time  =  6.26  seconds 
Peak  value  =  4.08 
Pinal  value  =  4.00 


s.  .  -v  *. 
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The  lower  bound  pitch  angle  response  for  a  4  degree 


^itch  angle  command  is  shown  in  Table  D-8. 


TABLE  0-8 


a__ — Time  Domain  Specifications 
■^KC-llS  Longitudinal  Mode 


Rise  time  =  5.08  seconds 
Settling  time  =  10.1  seconds 
Peak  value  =  4.00 
Final  value  =  4.00 


Table  d-9  shows  the  maximum  allowable  response  for 
perturbation  altitude,  h,  due  to  a  pitch  angle  command  of  4 
degrees . 


TABLE  D-9 

b, 2 — Time  Domain  Specifications 
KC-135  Longitudinal  Mode 


Rise  time  =  2.11  seconds 
Settling  time  =  9.92  seconds 
Peak  value  =  0.5 
Final  value  =  0 
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The  maximum  response  for  perturbation  velocity,  u,  due 


to  a  pitch  angle  command  is  outlined  in  Table  D-10. 


TABLE  D-10 

b22--Time  Domain  Specifications 
KC-135  Longitudinal  Mode 


Peak  time  =  1.94  seconds 
Settling  time  =  9.76  seconds 
Peak  value  =  1.00 
Final  value  =  0 


F-16  Response  Models 

Pitch  Rate  Command 

The  upper  bound  pitch  rate  response  to  a  pitch  rate 
command  of  10  degrees  per  second  is  contained  in  Table-11, 


TABLE  D-11 

b^2^--Time  Domain  Specifications,  F-16 


Rise  time  =  0.65  seconds 
Settling  time  =  3.01  seconds 
Peak  value  =  55.00  deg/sec 


Rise  time  =  0.65  seconds 
Settling  time  =  4.30  seconds 
Peak  value  =  1.4  deg/sec 
Total  pitch  angle  _<  3.0  deg 


Appendix  E:  Continuous  and  Discrete  Compensators 
Introduction 

This  appendix  contains  the  continuous  and  discrete 
compensators  for  the  longitudinal  mode  of  the  KC-135  and  the 
F-16  aircraft.  The  continuous  domain  compensators  are 
modified  to  counter  the  effects  of  the  ZOH,  as  discussed  in 
Chapter  III.  The  KC-135  longitudinal  mode  compensators  are 
used  for  both  the  rigid  and  elastic  body  models.  The  sampling 
rates  for  the  KC-135  and  F-16  are  selected  as  40  Hz  and  60  Hz 


respectively.  Also  included  are  the  S  domain  prefilters  for 
the  F-16  and  both  modes  of  the  kc-135  as  originally  designed 
in  references  1  and  2. 


Continuous  Domain  Prefilters 


KC-135  Lateral  Mode 


.9 


(S+.9) (S+1) 
.512 


(S+. 512) (S+1) 


KC-135  Longitudinal  Mode 


F22(s)  - 


.3506 


(S+.5) (S+.54+j.64) 


F-16 


Fii(s)  = 


S+2 
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(E-1) 


(E-2) 


(E-3) 


(E-4) 


F22(s) 


(S+1.334+jl. 00022) 


Continuous  Domain  Compensators 
KC-135  Longitudinal  Mode 

^  .  .  -18.1875(3+4) {S+80) 

o,  i  S  i  “  ■  .  -  II 

(S+24+j32) 

=  5.62(1o’^)(S+1)(S+80) 

^  (S+30) (S+390+j520) 

_  ,  ,  2680(3+2) (S+. 01846) (S+80) 

b^tS)  =  - 

S(S+.5) (3+10) (3+18+ j24) 


F-16 


(E-6) 


(E-7) 


(E-8) 


-19916.667(3+. 3) (3+1.1) (3+1.167) (3+2.028) (3+20) 
S^(S+.0107) (S+.4845) (S+.5)^(S+35) 


(3+120) (3+430) (S+. 006472+ j. 07803) 
(3+87) (S+636+j848) 

-18750(S+.03448) (S+. 3) (3+2.171) (3+12.9) (5+120)^ 
S^(S+.2) (3+3.46) (3+48) (S+. 62+ j6. 07582) 

(S+.4996+j3. 12896) 

(S+450+j600) 


Digital  Controllers 

KC-135  Longitudinal  Mode 


-20. 6453Z(Z-. 904762) 
(Z-.405405+j. 432432) 


(E-11) 


13477. 3Z(Z-.975309)(Z-H) _ 

(Z-.454545) (Z+. 846937+3,169347) 

. 476 851Z ( Z-. 999539 ) (Z-. 95 1220) 
(Z-1) (Z-. 987578 ) (Z-. 777778) 

(Z+1)^ 

(Z-.540275+j. 377210) 


(E-12) 


(E-13) 


9. 1 198 3Z ( Z+1 ) (Z-. 995012 ) (Z-. 981833) (Z-. 980737) 

( Z-1 )2(Z-. 999822 ) (Z-.991957) (Z-.991701)^ 

( Z-. 966762 ) (Z-. / 14286 ) (Z+. 563636) 

(Z-. 548387) (Z-. 159420) 

(Z-.999891+j. 001300)  ^ 

(Z+.859419+j. 157689) 

10. 258 5Z^ ( Z+1 ) (Z-. 999425 )(Z-. 995012 ) (Z-. 964460) 
(Z-1 )^(Z-. 996672 ) (Z-. 943949 ) (Z-.428571) 


(  Z-. 805869 ) (Z-. 990366+ j. 051683) 
(Z-.984684+j.099972) (Z+. 80026 3+ j. 210250) 


(E-15) 


Appendix  F;  Continuous  System  Simulation  Responses 


Introduction 

This  appendix  contains  the  simulation  responses  for  the 
KC-135  longitudinal  mode  and  the  F-16  continuous  systems. 

For  the  KC-135  longitudinal  mode/  continuous  system  simula¬ 
tions  are  performed  only  for  the  rigid  body.  The  KC-135 
lateral  mode  continuous  system  responses  are  shown  in  Chapter 


The  continuous  system  response  data  in  this  appendix  and 
in  Chapter  IV  is  compared  with  the  original  continuous  system 
performance  in  references  1  and  2  as  an  empirical  validation 
of  the  plant  and  compensator  data  as  well  as  the  general 
simulation  setup  used  in  this  study.  The  continuous  response 
data  is  also  compared  with  the  hybrid  system  responses  to 
assess  the  hybrid  system  performance. 


KC-135  Longitudinal  Mode  Responses  -  Continuous  System 


The  pitch  angle  command  responses  are  shown  in  Figures 
F-1  through  F-3.  The  pitch  rate  response  is  shown  in  Figure 
F-1.  Figure  F-2  shows  the  perturbation  altitude  response  and 
Figure  F-3  shows  the  perturbation  velocity  response.  The 
rigid  body  response  characteristics  for  the  pitch  angle 
command  are  shown  in  Table  F-1.  As  noted  in  Section  IV-4, 
comparison  of  these  responses  with  the  specifications  shows 
that  the  desired  robust  performance  is  achieved. 


Continuous  System  Altitude  Responses 


F-16  Pitch  Rate  Response 


Continuous  System 


The  pitch  cate  command  response  for  the  F-16  is  shown  in 
Figures  F-4  and  F-5.  The  pitch  rate  response  is  shown  in 
Figure  F-4  and  the  roll  angle  response  is  shown  in  Figure 
F-5.  Note  that  all  four  F.C.'s  are  included  in  these  response 
plots.  From  Figure  F-5,  the  maximum  roll  angle  for  the  pitch 
rate  command  is  approximately  2.8  degrees  which  meets  the  5 
degree  specification.  The  figures  of  merit  for  the  pitch 
rate  command  responses  are  shown  in  Table  F-2.  Examination 
of  the  continuous  system  response  data  shows  that  in  general, 
the  specifications  are  met.  In  some  instances  the  settling 
time  specification  is  not  met,  but  as  discussed  in  reference 
2  the  specifications  can  be  met  by  fine  tuning  the  prefilter, 
F^j^.  Note  that  the  maximum  delta  in  the  settling  time  is 
1.33  seconds  which  is  less  than  the  delta  in  the  response 
models  (1.5  secs).  The  pitch  rate  response  data  in  reference 
2  showed  a  2.14  second  delta  in  the  settling  time  which 
theoretically  should  not  occur.  It  is  highly  probable  that 
the  use  of  reduced  order  compensators  in  reference  2  caused 
the  discrepancy. 


PITCH 


Table  F-2 


Pitch  Rate  Response  Characteristics  -  Continuous  System 


Peak  Value  Rise 

F.C.  CSC  Mode  (deg/sec)  Time  (sec) 


Settling 
Time  (sec) 


1 

10.02 

1.00 

1.80 

2 

10.12 

0.77 

1.25 

3 

10.02 

0.90 

1.80 

4 

10.12 

0.71 

1.25 

5 

10.30 

0.85 

2.50 

6 

10.02 

0.90 

1.75 

1 

10.00 

1.02 

1.75 

2 

10.02 

0.90 

1.50 

3 

10.00 

1.02 

1.79 

4 

10.13 

0.82 

1.34 

5 

10.03 

0.86 

1.45 

6 

10.00 

1.02 

1.75 

1 

10.00 

1.09 

1.92 

2 

9.99 

1.04 

1.88 

3 

10.00 

1.08 

1.94 

A 

** 

10.03 

1.02 

1.60 

5 

10.00 

1.05 

1.92 

6 

10.00 

1.08 

1.99 

1 

10.02 

1.28 

2.17 

2 

10.05 

1.38 

2.38 

3 

10.02 

1.22 

2.17 

4 

10.12 

1.57 

2.58 

5 

10.05 

1.35 

2.40 

6 

10.02 

1.28 

2.17 

F-16  Roll  Rate  Command  Response  -  Continuous  System 

The  roll  rate  command  responses  are  shown  in  Figures 
F-6  and  F-7.  Figure  F-6  shows  the  roll  rate  response  and 
Figure  F-7  shows  the  pitch  angle  response.  The  maximum  pitch 


angle  is  approximately  -1.9  degrees.  The  pitch  angle  specifi 
cation  is  3  degrees,  therefore  this  specification  is  met. 
Table  F-3  shows  the  roll  rate,  response  characteristics.  The 
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ystem  Roll  Rate  Responses 


Introduction 


This  appendix  contains  the  microprocessor  emulation 
subroutine  code  listing.  Also  included  is  the  code  listing 
for  the  utility  "hook"  subroutine  provided  with  the  MATRIX^^ 
executable  code  (7:SB  6-2).  A  brief  discussion  of  each 
subroutine  is  given  below. 

Hook  Subroutine 

A  detailed  explanation  of  this  subroutine  and  its  use 
is  contained  in  reference  7  (7:SB  6-4).  Essentially  the 
entire  code  of  this  subroutine  is  provided  with  the  MATRIX^ 
code.  The  name  of  the  subroutine  is  UPDFTN.  If  more  than 
one  user  written  FORTRAN  subroutine  is  to  be  interfaced  with 
MATRIXj^  then  this  UPDFTN  must  be  modified  slightly.  The 
version  of  UPDFTN  shown  below  has  been  modified  to  allow 
three  different  subroutines  to  be  called  from  MATRlXj^  (actual 
ly  from  System  Build)  where  their  names  can  be  seen  in  the 
calling  statements  in  the  code.  All  parameters  sent  from 
System  Build  to  the  microprocessor  emulation  subroutine  (i.e. 
the  "user"  defined  subroutine)  pass  untouched  through  subrou¬ 
tine  UPDFTN.  It  is  very  important  that  none  of  the  UPDFTN 
code  be  changed  (from  the  template  provided  with  MATRlXj^) 
except  in  very  specific  areas.  These  areas  are  explained  in 
reference  (7).  The  code  for  UPDFTN  is  given  below. 


non  ooo  nnn 


SUBROUTINE  UPDFTN  (  NUMBER,  NAME,  T,  U,  NU,  X,  XDOT,  NX,  Y,  NY, 

+  RP,  IP,  DELX,  DELY,  ERROR  ) 

C 

C - 

C  !  SYSTEM_BUILD  FORTRAN  FUNCTION  EVALUATION  HOOK  ! 

C  !  ! 

C  !  Used  by  MATRIXX,  An  ANSI77  version,  A  subset  of  System_Bui Id .  ! 

C - 

C 

DOUBLE  PRECISION 
DOUBLE  PRECISION 
INTEGER 
CHARACTER*6 
LOGICAL 
DATA 
DATA 
C 

BLOCK  =  IABS( NUMBER) 

Test  whether  or  not  NUMBER  is  in  range  : 

IF(  BLOCK. GT.MAXNUM  )  THEN 
ERROR  =  .TRUE. 

IF(  NUMBER. GT.O  )  WRITE(6,111)  BLOCK 
111  FORHAT(/'  >>UPDFTN>>  Not  able  to  update  Function  ',12) 

RETURN 
ENDIF 

If  NUMBER  is  negative  then  just  pass  back  block  name  ; 

NAME  =  FNAME( BLOCK) 

IF(  NUMBER  .LT.  0  )  RETURN 


DBV  09/05/84 

T,  U(l),  X(l),  Y(l) 

RP(l),  XDOTd),  DELX(l),  DELY(l) 

IP(1),  BLOCK 
NAME,  FNAMEO) 

ERROR 

MAXNUM  /  3  / 

FNAME(l)  ,FNAME(2)  , FNAMEO)  /  'GlZ','G2Z'  ,'G3Z'  / 


Fix  up  the  GO  TO  statement  to  branch  to  various  FORTRAN  blocks  : 


GO  TO  (  01, 

CM 

O 

03 

), 

BLOCK 

01 

CALL  FTNOl 
GO  TO  99 

{  T, 

u, 

NU, 

X, 

DELX, 

NX, 

DELY, 

NY, 

RP, 

IP, 

ERROR 

) 

02 

CALL  FTN02 
GO  TO  99 

(  T, 

u, 

NU, 

X, 

DELX, 

NX, 

DELY, 

NY, 

RP, 

IP, 

ERROR 

) 

03 

CALL  FTN03 

(  T, 

u, 

NU, 

X, 

DELX, 

NX, 

DELY, 

NY, 

RP, 

IP, 

ERROR 

) 

GO  TO  99 


C 


C 

C 

c 

04 

CALL 
GO  TO 

FTN04 

99 

{  T, 

u. 

NU, 

X, 

DELX, 

NX, 

DELY, 

NY, 

RP, 

IP, 

ERROR 

) 

c 

05 

CALL 

FTN05 

(  T, 

u. 

NU, 

X, 

DELX, 

NX, 

DELY, 

NY, 

RP, 

IP, 

ERROR 

) 

C  GO  TO  99 
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Convert  to  implicit  form  : 

99  IF(  NX  .GT.  0  )  THEN 
DO  998  I  =  1,  NX 

998  DELX(I)  =  DELX(I) 

ENDIF 

IF(  NY  .GT.  0  )  THEN 
no  999  I  =  1,  NY 

999  DELY(I)  =  DELY(I) 

ENDIF 

RETURN 

END 


-  XDOT(I) 


-  Y(I) 


Microprocessor  Emulation  Code 


System  Build  has  the  capability  to  pass  all  parameters 
needed  by  the  microprocessor  emulation  proqram.  For  different 
simulations  it  is  desirable  to  change  certain  parameters  for 
the  emulation  subroutine.  For  instance,  the  wordlength  or 
rounding/truncation  specification  might  be  changed  from  one 
simulation  to  the  next.  Unfortunately,  if  one  parameter  is 
changed,  all  parameter  must  be  reentered  into  System  Build. 

The  emulation  program  was  originally  written  with  all  para¬ 
meters  being  passed  from  System  Build  but  because  of  the 
above  mentioned  limitation,  the  code  was  modified  so  that  the 
difference  equation  coefficients  are  coded  into  the  subroutine 
itself.  This  alleviates  the  need  to  continuously  reenter  the 
coefficients.  The  following  parameters  are  specified  by 
the  user  and  passed  to  the  emulation  code  from  System  Build; 
the  rounding/truncation  specification,  the  microprocessor 
wordlength,  the  number  of  bits  to  the  right  of  the  radix 
point,  and  the  order  of  the  difference  equation.  These  items 


are  passed  in  an  integer  'parameter  vector"  called  IP.  in 
addition.  System  Build  passes  the  following  items  to  the 
emulation  subroutine  in  support  of  the  dynamic  discrete 
simulation:  the  current  sample  time  and  the  floating  point 

number  representing  the  error  signal  (i.e.  the  input  to  the 
controller).  Other  items  are  passed  which  are  not  used. 

These  unused  items  simply  assume  their  default  values.  See 
the  definition  of  variable  section  of  the  code  for  specific 
variable  names.  After  the  emulation  subroutine  computes  the 
difference  equation  it  passes  back  the  compensation  signal  to 
System  Build  as  a  floating  point  number. 

It  is  critical  for  the  emulation  subroutine  code  to  be 
written  following  certain  conventions  expected  by  system 
Build.  For  a  detailed  explanation  of  these  conventions  see 
reference  7  (7:SB  6-3).  Pseudocode  and  an  explanation  of  each 
variable  is  given  in  the  code  below.  Note  that  a  different 
subroutine  is  used  for  each  digital  controller.  This  greatly 
simplifies  the  interface  with  System  Build,  so  for  the 
KC-135  longitudinal  mode  there  are  three  subroutines  with 
identical  code  except  for  the  difference  equation  coeffi¬ 
cients.  The  particular  version  shown  below  is  for  G2(z)  of 
the  KC-135  longitudinal  mode.  Note  also  that  double  precision 
arithmetic  is  used  throughout  the  subroutine  as  described  in 
Chapter  V.  The  microprocessor  emulation  code  follows: 


SUBROUTINE  FTN02  (  T,  U,  NU,  X,  XDOT,  MX,  Y,  NY, 

>  RP,  IP,  ERROR  ) 

C 

C - 

Cl  DESCRIPTION 

Cl  This  subroutine  interfaces  with  MATRIX  x  -  SYSTEM  BUILD  to 
Cl  simulate  the  finite  wordlength  fixed  point  binary 
Cl  arithmetic  (addition  and  multiplication)  of  a  discrete  control 
Cl  compensator  implemented  on  a  microprocessor.  Currently,  the 
Cl  subroutine  is  capable  of  simulating  only  equal  order 

Cl  compensators,  but  this  limititation  can  be  removed  easily. 

Cl  The  highest  order  denominator  coefficient  is  assumed  to  be  one. 
Cl  The  coefficients  must  be  coded  in  this  subroutine  and  ordered 
Cl  from  highest  powered  coefficient  to  lowest  powered  coefficient. 

Cl  The  leading  denominator  coefficient  (equal  to  one)  must  be 

Cl  entered  also. 

C - 

c 

c - 

Cl  PSEUDOCODE 

Cl  -  initialize  parameters 

Cl  -  initialize  compensator  coefficients 

Cl  -  massage  compensator  coefficients  (integer  form) 

Cl  -  compute  overflow  values 

C|  -  massage  current  input  to  compensator  (error  signal) 

Cl  -  compute  input  portion  of  difference  equation 

Cl  --  massage  multiplies 

Cl  —  check  for  overflow 

Cl  --  add  terms 

Cl  —  check  for  overflow 

Cl  -  compute  output  portion  of  difference  equation 
Cl  --  same  as  input  portion  except  subtract 

Cl  -  shift  outputs  and  inputs 

Cl  -  compute  'microprocessor'  output  (floating  point  number) 

c - 

c 

c - 

Cl  DEFINITION  OF  VARIABLES 
Cl 

Cl  T  -  current  sample  time 

Cl  U  -  input  vector  (error  signal) 

Cl  NU  -  number  of  inputs  (dimension  of  U;  =1) 

Cl  X  -  state  vector  (not  used) 

Cl  XDOT  -  first  derivative  of  state  vector  (not  used) 

Cl  NX  -  number  of  states  (dimension  of  X;  =  0) 

Cl  Y  -  output  vector  (compensation  signal) 

Cl  NY  -  number  of  outputs  (dimension  of  Y;  =  1) 

Cl  RP  -  general  vector  of  real  values  (not  used) 

Cl  IP  -  general  vector  of  integer  parameters 

Cl  ERROR  -  logical  flaq  (TRUE  means  simulation  error,  not  used) 

Cl 

Cl  ZZA  -  numerator  coefficient  vector 


o  n  n  n  n 


Cl  ZZB  -  denominator  coefficient  vector 
Cl  ZZRT  -  rounding/truncate  value  (0  or  ,5) 

Cl  ZZSIGN  -  sign  of  value  to  be  'massaged' 

Cl  ZZL  -  microprocessor  wordlength  minus  1 
Cl  ZZN  -  numerator  order 

Cl  ZZP  -  number  of  digits  to  right  of  radix  point 
Cl  ZZI  -  generic  counter 
Cl  ZZOVRH  -  positive  overflow  value 
C|  ZZOVRL  -  negative  overflow  value 

Cl  ZZR  -  vector  of  current  &  past  compensator  inputs 
Cl  ZZC  -  vector  of  current  and  past  compensator  outputs 
Cl 

Cl  NOTE;  All  variables  unique  to  this  subroutine  begin  with  the 
Cl  letters  'ZZ'  to  insure  uniqueness  with  MATRIX  X 

Cl  variables.  All  other  variables  are  initialized  in 

Cl  MATRIX  X  and  passed  to  this  subroutine  as  arguments, 

c - 

c 

c - 

Cl  SPECIFICATIONS 

C - 

C 

IMPLICIT  DOUBLE  PRECISION  (A-H,0-Z) 

DOUBLE  PRECISION  T,  U(l),  X(l),  XDOT(l),  Y(l),  RP(1) 

DOUBLE  PRECISION  ZZA(0;10),  ZZB(0:10),  ZZRT,  ZZSIGN 
DOUBLE  PRECISION  ZZR(OzlO),  ZZC(0:10),  ZZS,  ZZT,  ZZOVRH,  ZZOVRL 
INTEGER  IP(1),  ZZL,  ZZN,  ZZP,  ZZI 
LOGICAL  ERROR 
COMMON  /ZZFTN2/  ZZR,  ZZC 


I  INITIALIZATION  OF  PARAMETERS 


IF  (  T  .EQ.  O.OD+00  )  THEN 
DO  10  ZZI  =  0,ZZN 

ZZR(ZZI)  =  O.OD+00 
ZZC(ZZI)  =  O.OD+00 
10  CONTINUE 

ENDIF 
C 

ZZL  =  IP(1)  -  1 
ZZP  =  IP(2) 

IF  (  IPO)  .EQ.  1  )  THEN 
ZZRT  =  0.5D+00 

ELSE 

ZZRT  =  O.OD+OO 

ENDIF 

C 

C - 

Cl  INITIALIZATION  OF  COMPENSATOR  COEFFICIENTS 


c 

DO  30  ZZI  =  0,ZZN 

IF  (  ZZA(ZZI)  .LT.  O.OD+00  )  THEN 

•;  ZZSIGN  =  -l.OD+00 

V  ■■  else 

ZZSIGN  =  l.OD+00 

ENDIF 

ZZA(ZZI)  =  DINT(  ZZA(ZZI)  *  { 2 . OD+00* *Z ZP )  +  ZZRT 
+  *  ZZSIGN  ) 

IF  (  ZZA(ZZI)  .GT.  ZZOVRH  )ZZA(ZZI)  =  ZZOVRH 

IF  (  ZZA(ZZI)  .LT.  ZZOVRL  )ZZA(ZZI)  =  ZZOVRL 

30  CONTINUE 
C 


Cl  MASSAGE  DENOMINATOR  COEFFICIENTS  (INTEGER  FORM) 

C - 

c 

DO  40  ZZI  =  0,ZZN 

IF  (  ZZB(ZZI)  .LT.  O.OD+00  )  THEN 
ZZSIGN  =  -l.OD+00 

ELSE 

ZZSIGN  =  l.OD+00 

ENDIF 

ZZB(ZZI)  =  DINT(  ZZB(ZZI)  *  ( 2 . OD+00* *Z ZP )  +  ZZRT 
+  *  ZZSIGN  ) 

IF  (  ZZB(ZZI)  .GT.  ZZOVRH  )ZZB(ZZI)  =  ZZOVRH 
IF  (  ZZB(ZZI)  .LT.  ZZOVRL  )ZZB(ZZI)  =  ZZOVRL 
40  CONTINUE 


nonnn  nnno 


I  MASSAGE  CURRENT  INPUT 


IF  (  U(l)  .LT.  O.OD+00  )  THEN 
ZZSIGN  =  -l.OD+00 

ELSE 

ZZSIGN  =  l.OD+00 

ENDIF 

ZZR(O)  =  DINT(  U(l)  *  ( 2 .0O+00**ZZP)  +  ZZRT  *  ZZSIGN  ) 


I  COMPUTE  INPUT  SEQUENCE 


C 

C 

C 

C 

C 


C 


ZZS  =  O.OD+00 
DO  50  ZZI  =  0,ZZN 

ZZT  =  ZZA(ZZI)  *  ZZR(ZZI) 

IF  (  ZZT  .LT.  O.OD+00  )  THEN 
ZZSIGN  =  -l.OD+00 

ELSE 

ZZSIGN  =  l.OD+00 

ENDIF 

ZZT  =  DINT(  ZZT  *  ( 2 . 0D+00**-ZZP )  +  ZZRT  *  ZZSIGN  ) 
IF  (  ZZT  .GT.  ZZOVRH  )ZZT  =  ZZOVRH 

IF  (  ZZT  .LT.  ZZOVRL  )ZZT  =  ZZOVRL 

ZZS  =  ZZS  +  ZZT 

IF  (  ZZS  .GT.  ZZOVRH  )ZZS  =  ZZOVRH 

IF  (  ZZS  .LT.  ZZOVRL  )ZZS  =  ZZOVRL 

50  CONTINUE 


COMPUTE  OUTPUT  SEQUENCE 


DO  60  ZZI  =  i,ZZN 

ZZT  =  ZZB(ZZI)  *  ZZC(ZZI) 

IF  (  ZZT  .LT.  O.OD+00  )  THEN 
ZZSIGN  =  -l.OD+OO 

ELSE 

ZZSIGN  =  l.OD+OO 

ENDIF 

ZZT  =  DINT(  ZZT  *  ( 2 . 0D+00**-Z ZP )  +  ZZRT  *  ZZSIGN  ) 


IF 

( 

ZZT  .GT.  ZZOVRH 

)ZZT 

=  ZZOVRH 

IF 

( 

ZZT  ,LT.  ZZOVRL 

)ZZT 

=  ZZOVRL 

ZZS 

=  ZZS  -  ZZT 

IF 

( 

ZZS  ,GT.  ZZOVRH 

)ZZS 

=  ZZOVRH 

IF 

( 

ZZS  .LT.  ZZOVRL 

)ZZS 

=  ZZOVRL 

60  CONTINUE 


Pt 


SHIFT  PAST  INPUTS  AND  OUTPUTS 


ZZC(O)  =  ZZS 

DO  70  ZZI  =  ZZN,1,-1 

ZZC(ZZI)  =  ZZC(ZZI-l) 
ZZR(ZZI)  =  ZZR(ZZI-l) 
70  CONTINUE 


I  COMPUTE  'MICROPROCESSOR'  OUTPUT 


Y(l)  =  ZZC(O)  *  (2.0D+00**-ZZP) 
ERROR=  .FALSE. 


RETURN 
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